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CFD  Validation  for  Propulsion  System  Components 

(AGARD  AR-355) 

Executive  Summary 

The  high  performance  of  aircraft  turbine  engines  -  mainly  oriented  towards  economy  in  civil  engines 
and  towards  performance,  economy  and  off-design  reliability  in  military  applications  -  depends  on  the 
correct  design  of  gas  flows  in  the  engine.  Today,  this  is  mainly  done  by  numerical  computer  simulation 
instead  of  the  much  more  costly  experimentation  with  hardware.  The  working  group  undertook  to 
analyse  the  quality  of  computer  codes  in  use  and  to  show  ways  for  improvement.  Many  discrepancies 
between  codes  were  detected  and  it  was  recognized  that  even  the  researchers  actively  working  in  this 
field  were  unaware  of  the  full  scale  of  this  phenomenon.  The  military  will  gain  from  the  newly  designed 
engines  as  well  as  from  possible  improvement  in  upgrades.  It  should  be  noted  that  this  kind  of  down-to- 
earth  analysis,  despite  its  merits,  has  not  been  undertaken  on  an  international  basis  outside 
AGARD/RTO. 
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La  validation  CFD  des  organes  des  propulseurs 

(AGARD  AR-355) 


Synthese 

La  recherche  de  hautes  performances  des  turbomoteurs,  motivee,  pour  les  moteurs  civils, 
principalement  par  la  recherche  d’ economies,  et  pour  les  moteurs  militaires  par  F  amelioration  des 
performances,  les  economies  et  la  fiabilite  dans  des  conditions  hors  tolerances,  passe  par  une 
connaissance  precise  de  la  conception  des  ecoulements  dans  le  moteur. 

Aujourd’hui,  dans  ce  domaine  il  est  fait  appel  en  general,  a  la  simulation  numerique  par  ordinateur,  de 
preference  a  F experimentation  reelle,  beaucoup  plus  couteuse.  Le  groupe  de  travail  No.  26  s’est  donne 
pour  tache  d’ analyser  la  qualite  des  codes  machine  utilises  et  de  faire  des  recommendations  concemant 
d’eventuelles  ameliorations.  Un  nombre  eleve  de  divergences  entre  les  differents  codes  a  ete  constate  et 
il  a  ete  admis  que  meme  les  chercheurs  impliques  activement  dans  ce  domaine  ignoraient  la  veritable 
ampleur  du  phenomene. 

Les  militaires  pourront  tirer  profit  des  moteurs  de  nouvelle  generation,  ainsi  que  des  ameliorations 
possibles  apportees  par  les  remises  a  niveau.  Il  est  a  noter  que  malgre  son  merite  evident,  ce  type 
d’analyse  pragmatique  n’a  pas  ete  entrepris  au  niveau  international  en  dehors  de  l’AGARD/RTO. 


iv 


Contents 


Page 

Executive  Summary  iii 

Synthese  iv 

Preface  vi 

Abstract  vii 

Recent  publications  of  the  Propulsion  and  Energetics  Panel  viii 

Introduction  1 

1.  Review  of  Turbomachinery  CFD  3 

1.1  Introduction  3 

1.2  Navier-Stokes  equations  for  a  rotating  system  3 

1.3  Mesh  generation  3 

1.4  Numerical  schemes  4 

1.5  Turbulence  modelling  6 

1.6  Code  validation  9 

2.  The  Test  Cases  11 

2.1  NASA  Rotor  37  11 

2.2  DLR  cascade  16 

3.  The  Computations  27 

3.1  Introduction  27 

3.2  The  results  for  Rotor  37  27 

3.3  The  results  for  the  DLR  cascade  49 

3.4  Analysis  of  the  computational  grids  58 

3.5  Analysis  of  the  algorithms  64 

3.6  Analysis  of  the  turbulence  models  64 

3.7  Convergence  criteria  71 

4.  Conclusions  and  recommendations  75 

4.1  Conclusions  75 

4.2  Recommendations  for  research  77 

4.3  Recommendations  to  PEP  77 

Appendix  List  of  Working  Group  members  and  contributors  79 


References 


81 


Preface 


Computer  codes  which  solve  the  Reynolds-averaged  Navier-Stokes  (RANS)  equations  are  widely  used  by  gas 
turbine  and  steam  turbine  manufacturers  to  analyse  the  aerodynamic  performance  of  existing  compressors  and 
turbines  and  to  design  new  ones.  A  wide  range  of  computational  fluid  dynamics  (CFD)  approaches  and 
turbulence  models  have  been  developed,  and  each  code  user  validates  his  code  using  some  relevant  test  cases.  But 
the  degree  of  validity  of  any  code  for  application  to  a  wider  range  of  turbomachinery  configurations  than  those 
initial  test  cases  remains  open  to  question.  Experts  in  the  field  recognise  the  importance  of  using  a  “good” 
computational  grid  and  a  “good”  turbulence  model,  but  there  is  no  consensus  among  them  about  which  grids  and 
which  turbulence  models  are  “good”  enough  to  provide  a  reliable  basis  for  design  decisions. 

The  Propulsion  and  Energetics  Panel  set  up  Working  Group  26  to  help  the  international  research  community  to 
clarify  some  of  these  issues,  by  comparing  predictions  (using  as  wide  a  range  of  methods  as  possible)  of  two 
representative  but  difficult  test  cases.  The  Group  comprised  experts  in  this  field,  both  Panel  members  and  non¬ 
members.  Predictions  were  obtained  from  leading  code  developers  and  users,  including  some  of  the  Group 
members,  and  compared  with  measurements  of  the  two  test  cases  in  some  detail  by  all  of  the  members.  The 
Group  held  meetings  during  1994-6,  and  its  members  made  additional  contributions  to  the  Technical  Editor 
during  1997.  This  is  its  Report. 

The  Panel  greatly  appreciates  the  careful  work  of  all  those  code  developers  and  the  Working  Group  members. 


Dr  J.  Dunham  Dr  G.  Meauze 

Technical  Editor  Chairman 

Propulsion  and  Energetics  Panel 
Working  Group  26 
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Summary 


Computer  codes  which  solve  the  Reynolds-averaged  Navier-Stokes  equations  are  now  used  by  manufacturers  to 
design  turbomachines,  but  there  is  no  consensus  about  which  grids  and  which  turbulence  models  are  good  enough 
to  provide  a  reliable  basis  for  design  decisions.  The  AGARD  Propulsion  and  Energetics  Panel  set  up  Working 
Group  26  to  help  to  clarify  these  issues,  by  analysing  predictions  (using  as  wide  a  range  of  codes  as  possible)  of 
two  representative  but  difficult  single  blade  row  test  cases:  NASA  Rotor  37  and  an  annular  turbine  cascade  tested 
by  DLR.  This  report  presents  the  Group’s  results  and  conclusions. 

The  predicted  performance  of  both  test  cases  fell  short  of  the  accuracy  engine  designers  need.  NASA  discovered 
that  the  comer  stall  observed  at  the  hub  of  Rotor  37  was  affected  by  the  presence  of  a  small  axial  gap  in  the  hub 
annulus  line  just  upstream  of  the  rotor,  not  modelled  by  the  codes.  The  flow  in  the  tip  region  also  proved  too 
difficult  for  most  codes,  which  tended  to  overestimate  the  pressure  losses  there.  The  predicted  pressure  loss  of  the 
DLR  cascade  was  up  to  40%  in  error,  and  some  codes  were  unable  to  predict  correctly  the  highly  three- 
dimensional  secondary  flow. 

Recommendations  are  made  about  the  type  and  density  of  grid,  which  depend  on  many  factors.  Mixing-length 
turbulence  models  are  unsuitable  for  turbomachines  with  their  complex  endwall  flows;  some  kind  of  turbulent 
transport  model  is  essential.  No  turbulence  model  was  found  which  always  gave  good  loss  predictions. 
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INTRODUCTION 


Computational  methods  for  the  design  and  analysis  of 
turbomachine  flows  have  been  developed  and  brought  into  use 
progressively  over  the  last  fifty  years.  The  coming  of  large 
modem  computers  in  recent  years  has  allowed  the  practical 
development  of  codes  which  solve  the  Reynolds-averaged 
Navier-Stokes  (RANS)  equations  in  three  dimensions.  Such 
codes  are  already  used  by  the  large  engine  manufacturers  for 
the  advanced  design  of  some  engine  components  like 
compressors,  turbines,  air  intakes,  nozzles,  and  combustion 
chambers.  Notable  improvements  have  already  been  obtained; 
nevertheless  the  physical  representation  of  real  phenomena  is 
not  yet  completely  satisfactory,  and  only  comparative 
predictions  can  be  seriously  considered. 

Different  computational  fluid  dynamics  (CFD)  approaches  and 
turbulence  models  exist,  and  it  seems  essential  today  to 
establish  their  degree  of  validity  for  application  to  typical 
configurations  in  turbomachinery.  Obviously,  the  different 
CFD  code  users  make  their  own  evaluation,  but  that  remains 
limited.  To  make  a  real  improvement  in  knowledge  of  the 
influence  of  the  different  elements  of  a  CFD  code  (numerical 
algorithm,  type  and  density  of  mesh,  turbulence  model,  ...)  on 
the  results  it  is  necessary  to  compare  computations  from 
different  codes.  This  very  important  work  can  only  be  done  in 
an  international  field;  it  fits  AGARD’s  objectives  perfectly. 

PEP  has  regularly  organised  meetings,  Lecture  Series,  or 
Working  Groups  on  turbomachinery  flow  prediction  methods. 
In  1976,  Lecture  Series  83  was  held  (Serovy,  1976),  and  also  a 
Specialist  Meeting  in  which  predictions  for  a  set  of  test  cases 
were  invited  (AGARD,  1976).  Subsequently,  PEP-WG  12  was 
set  up  to  review  through-flow  calculation  methods,  which 
reported  in  1981  (Hirsch  and  Denton,  1981);  it  also  employed 
test  cases.  In  1985,  a  Lecture  Series  on  the  emerging  three- 
dimensional  flow  computation  methods  applied  to 
turbomachines  was  presented  (G.Meauze,  1985).  Following 
that,  PEP-WG  18  was  set  up  to  assemble  another  new  group  of 
test  cases,  published  in  1990  (Fottner,  1990).  A  Symposium  on 
CFD  techniques  for  propulsion  applications  was  held  by  the 
Propulsion  and  Energetics  Panel  in  Spring  1991  (AGARD, 
1991).  In  his  Technical  Evaluation  Report,  Dr  L.A.Povinelli 
recommended  that  a  new  activity  should  be  started  aimed  at 
code  validation.  Other  similar  recommendations  had  also  been 
made  to  the  Panel,  and  following  discussions,  it  was  agreed  to 
form  Working  Group  26:  “CFD  validation  for  propulsion 
system  components”.  This  was  approved  by  the  National 
Delegates’  Board  in  Fall  1993,  for  starting  in  1994.  The  topic 
falls  within  the  specific  recommendations  of  AGARD 
concerning  improved  engine  design  (flexibility),  improved 
affordability,  and  improved  hardware  and  software  reliability. 

A  Lecture  Series  on  turbomachinery  CFD  was  held  in  1994 
(Hirsch,  1994). 

Meanwhile,  several  experts  had  been  consulted,  and  detailed 
plans  about  the  Working  Group  discussed  informally.  Due  to 
the  scope  of  the  work  foreseen,  it  was  agreed  to  limit  the  range 
of  investigation  to  viscous  three-dimensional  steady  and  non¬ 
reacting  flow  configurations.  Two  representative 
turbomachinery  configurations  were  chosen  as  test  cases  for  the 
principal  studies.  Particular  attention  was  given  to  the  accuracy 
of  the  experimental  results,  hr  1993,  the  Turbonrachinery 
Committee  of  the  International  Gas  Turbine  Institute  of  the 
American  Society  of  Mechanical  Engineers  (ASME)  had  issued 
an  open  invitation  to  predict  the  flow  details  of  an  isolated 


transonic  fan  rotor,  NASA  Rotor  37.  This  proved  a  challenging 
case,  so  the  WG  decided  to  select  it  as  one  of  the  two  cases;  the 
other  was  an  annular  turbine  cascade  tested  by  the  German 
Research  Establishment  DLR. 

The  objectives  of  the  Working  Group  were  to  obtain  CFD 
calculations  of  the  specified  flows,  to  evaluate  the  results  and 
the  methods  used,  for  the  purpose  of  advancing  the  technology. 
It  was  intended  to  study  the  effectiveness  of  the  methods  used 
in  the  various  CFD  codes,  including  grid  geometry  and  type 
and  the  turbulence  model,  and  to  make  recommendations  for 
the  guidance  of  code  developers,  and  for  future  research. 

Most  of  the  test  case  activities  have  been  led  by  experts  from 
Government  agencies  and  Universities.  Codes  are  widely  used 
by  gas  turbine  and  steam  turbine  manufacturers,  but  most  aero 
engine  firms  have  usually  regarded  their  proprietary  methods 
as  too  commercially  sensitive  to  publish  their  predictions  of 
test  cases.  In  selecting  which  code  to  use,  organisations  have  to 
consider  several  factors: 

( 1 )  experience  within  the  organisation  of  using  the  code; 

(2)  how  easy  it  is  to  use  and  how  robust  it  is; 

(3)  how  well-supported  it  is  by  its  originators  or  by  in- 
house  experts; 

(4)  what  computing  facilities  are  needed  and  how  long  it 
takes  to  run;  and 

(5)  how  accurately  it  predicts  the  aerodynamic 
performance  of  the  turbomachine. 

It  is  logical  to  expect  that  a  balance  must  be  found  between 
these  factors,  depending  on  the  particular  circumstances  of  the 
organisation  at  the  time,  hi  this  report,  information  about  the 
computing  requirements  and  time  (factor  4)  has  been  tabulated 
for  the  codes  used,  but  not  analysed;  the  focus  is  entirely  on 
accuracy.  So  it  was  not  intended  to  attempt  to  rank  the  specific 
codes  in  order  of  merit  or  to  recommend  to  manufacturers  a 
“best”  code;  indeed  such  a  recommendation  would  probably  be 
out-of-date  by  the  time  this  report  was  printed. 

The  report  has  been  written  following  detailed  discussions  by 
the  Group.  Each  section  was  then  drafted  by  one  or  two 
members,  listed  in  an  Appendix,  but  the  final  report  represents 
the  work  and  conclusions  of  all  the  members. 
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Chapter  1 

REVIEW  OF  TURBOMACHINERY  CFD 


1.1  INTRODUCTION 

The  development  of  numerical  methods  and  computing 
facilities  has  led  to  the  use  of  Computational  Fluid  Dynamics 
(CFD)  as  a  current  tool  for  designing  components  of  turbine 
engines.  Even  if  the  experimental  process  still  remains  the 
final  way  to  calibrate  and  validate  an  engine,  CFD 
simulations  allow  a  small  number  of  configurations  to  be 
selected  for  all  or  part  of  the  turbomachine  in  the  design 
process.  Different  stages  have  been  reached  in  the  use  of 
CFD:  3D  steady  turbulent  calculations  are  performed  daily, 
even  if  the  physical  models  used  today  cannot  reproduce 
faithfully  all  aspects  of  the  real  flow.  Despite  the  differences 
which  still  exist  between  numerical  simulations  and  reality,  it 
is  possible  to  predict  many  of  the  flow  properties  and  the 
losses  due  to  the  non-isentropic  features  of  the  flow  (shocks, 
viscous  layers,  tip  clearance  effects,  passage  vortices,...). 

This  chapter  is  a  short  presentation  of  the  state-of-the-art  of 
numerical  methods  used  for  solving  compressible  Navier- 
Stokes  equations  in  turbomachinery  applications.  This  review 
addresses  specific  features,  as  seen  from  the  standpoint  of 
turbomachinery  flows,  such  as  robustness  and  efficiency, 
space  accuracy  and  time  accuracy,  grid  density  and  grid 
regularity,  numerical  dissipation,  turbulence  model 
(including  transition).  This  chapter  is  far  from  exhaustive, 
and  the  reader  is  referred  to  Hirsch  (1994)  for  detailed 
presentations  concerning  turbomachinery  design  with  CFD, 
and  to  Couaillier  (1995)  and  to  AGARD  (1991)  for  a 
presentation  of  various  CFD  applications  to  turbomachinery 
flows.  Other  more  recent  references  will  be  indicated  as 
appropriate. 

1.2  NAVIER-STOKES  EQUATIONS  FOR  A 
ROTATING  SYSTEM 

This  section  is  devoted  to  the  formulation  of  the  Reynolds 
Averaged  Navier-Stokes  equations  applied  to  turbomachinery, 
flow  calculations,  the  full  3D  system  written  in  a  rotating 
frame  of  reference  (see  also  Veuillot,  1990). 

In  order  to  formulate  the  equations  in  a  rotating  frame  of 
reference,  the  definition  of  the  relative  time  derivative  is 
introduced.  Considering  a  scalar  function  in  space  and  time 
<D(M,  t) ,  if  the  point  M  is  fixed  in  a  rotating  wheel  with  a 

steady  speed  of  rotation  Q  ,  the  following  relationship  exists 
between  the  absolute  and  relative  time  derivatives  : 

The  relative  time  derivative  (dldt)R  is  equal  to  zero  for  any 
time-steady  scalar  field  in  the  rotating  frame  of  reference. 
Consider  now  the  time  derivative  of  a  vector  expressed  in  the 
relative  frame  of  reference.  Denoting  by  q.  the  coordinates 
of  Q  in  the  orthonormal  rotating  system  of  reference  (en ) , 
the  following  equation  is  obtained : 

(f)  -&)  <2> 

The  relative  time  derivative  (9/  dt)R  representing  the  first 
term  on  the  right-hand-side  of  this  equation  is  equal  to  zero 
for  any  stationary  vector  field  in  the  rotating  system  of 
reference. 


By  using  the  relative  time  derivative  definitions  and  by 
recalling  that  the  absolute  velocity  V  and  the  relative 
velocity  W  are  related  by:  V  =  W  +  Qx.r ,  r  being  the 
radius,  the  Navier-Stokes  equations  can  be  written  in  a 
different  form.  In  a  first  approach  the  governing  equations  for 
mass,  momentum,  and  energy  are  written  by  using  the 
relative  time  derivatives  of  the  absolute  variables  : 


l  dt  J 

1  +1 
R 

(3pE) 

+v[. 

V  8t  J 

a  1 

Another  approach  is  to  write  the  governing  equations  by 
using  the  relative  time  derivatives  of  the  relative  variables  : 

^ -  +v[pfF®fF-F-T(]  =  p(Q2r-2Qx#)  (4) 

(?^r)  +  1  w- ~  l ~  1  ] = P& 

The  stress  tensor  a  =  -pi  +  x  ,  where  x  is  the  shear  stress 
tensor  and  p  the  static  pressure,  is  not  affected  by  the  rotation 
(which  is  a  solid  body  motion).  In  system  (4),  as  opposed  to 
system  (3),  the  fluxes  are  formally  identical  to  those  of  the 
Navier-Stokes  equations  written  in  an  absolute  frame  of 
reference.  Nevertheless  the  discretisation  of  the  source  terms 
in  system  (4)  introduces  more  approximation  errors  than  the 
discretisation  of  the  source  terms  in  system  (3). 

The  expression  of  the  relative  specific  energy  conservation  in 
system  (4),  obtained  from  system  (3)  by  combining  the  energy 
equation  with  the  continuity  equation  and  the  momentum 
equation,  can  also  be  written  in  the  following  conservative 
form : 

ftr)  +  _  p  +  f,).#  _  ^  =  0  (5) 

In  equation  (5),  where  the  quantity  E'  is  equal  to 

Er  - 1/2  x  rj  ,  Er  being  the  relative  energy  (Cv7)),  it  is 

seen  that  the  rothalpy  I  is  constant  along  a  streamline  of  an 
inviscid  rotating  steady  flow. 


1.3  MESH  GENERATION 

The  generation  of  the  computational  mesh  is  certainly  today 
one  of  the  most  important  requirements  to  obtain  reliable 
numerical  solutions  for  turbomachinery  configurations.  It  is 
necessary  to  deal  with  various  geometries:  axial  or  centrifugal 
machines,  blades  with  thick  leading  or  trailing  edges,  highly 
cambered  turbines,  tip  clearance  representation,  rotor/stator 
interaction,  etc.  The  improvements  in  this  area  over  the  last 
ten  years  since  the  beginning  of  3D  calculations  are  very 
important.  A  large  variety  of  mesh  choices  have  been 
proposed,  from  a  simple  mono-block  C-grid  around  the  blade 
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to  a  complex  hybrid  multidomain  mesh.  Nevertheless,  it  has 
to  be  noted  that  generally  the  3D  meshes  are  built  from  a 
stacking  of  2D  blade-to-blade  meshes.  Moreover,  due  to 
complex  geometries  and  also  to  sophisticated  turbulence 
models,  an  increasing  number  of  mesh  strategies  take  into 
account  adaptive  refinement.  This  variety  of  approaches 
corresponds  to  different  levels  of  mesh  generation,  and  also 
to  different  levels  of  numerical  methods  in  terms  of  data 
structure  management.  By  considering  some  typical 
examples,  a  short  review  of  the  main  mesh  generation 
strategies  is  presented  now. 

1.3.1  Monodomain  structured  meshes 

With  a  mono-domain  approach,  three  types  of  grid  can  be 
considered.  An  H-grid  is  well  suited  to  applied  far-field  and 
periodicity  conditions  and  is  generally  easy  to  set  up,  but  is 
often  highly  skewed  near  the  leading  and  trailing  edges  of  the 
blades.  A  C-grid  provides  a  good  resolution  around  the 
leading  edge  and  in  the  wake,  but  becomes  skewed  at  the 
inflow  and  at  the  periodic  boundaries.  An  O-grid  allows  good 
resolution  of  both  leading  and  trailing  edges,  but  induces 
skewness  at  inflow,  outflow  and  periodic  boundaries. 

1.3.2  Multidomain  structured  meshes 

A  first  approach  used  to  remove  the  drawbacks  at  the  periodic 
boundary  consists  of  using  non-periodic  grids  (see  for 
instance  Veuillot  (1985)  and  Amone  (1993)).  A  more  general 
technique  is  based  on  a  structured  multidomain  approach, 
which  keeps  the  advantages  and  removes  the  drawbacks  of 
the  single  grids.  A  large  number  of  solvers  use  such  a 
technique  today,  with  the  capability  to  use  overlapping 
meshes.  However,  care  must  be  taken  to  ensure  that  the  grids 
are  smooth  across  the  block  interfaces. 

The  five-domain  mesh  presented  by  Heider  (1993)  for  a 
highly  cambered  turbine  with  a  tip  clearance  gap  shows  the 
advantages  of  such  an  approach:  for  each  blade-to-blade 
section  an  O-mesh  is  used  around  the  blade  and  H-meshes 
are  used  upstream  and  downstream  of  the  blade,  providing  a 
good  quality  of  mesh  everywhere.  In  the  clearance  an  H-O 
decomposition  allows  a  good  regularity  to  be  achieved  in  the 
radial  direction.  Choi  (1993)  and  Fougcres  (1994)  used  a 
multidomain  technique  with  or  without  overlapping  for 
meshing  the  holes  in  a  film  cooling  configuration.  This  type 
of  approach  allows  a  good  discretisation  of  the  viscous  layer 
and  of  the  hole  too,  in  order  to  capture  accurately  the  flow 
interaction.  Madavan  et  al  (1993)  use  multidomain 
overlapping  grids  for  a  turbine  rotor-stator  interaction. 

1.3.3  Unstructured  meshes 

The  use  of  an  unstructured  tetrahedra  approach  leads  to  a 
greater  flexibility  in  the  mesh  generation  process.  Moreover, 
even  if  a  multidomain  approach  with  overlapping  grids  or/and 
chimera  grids  allows  local  refinement,  unstructured  methods 
seem  to  be  more  appropriate.  Nevertheless,  unstructured 
solvers  are  still  today  less  efficient  in  CPU  time  than 
structured  solvers,  and  the  implementation  of  algebraic 
turbulence  models  leads  to  serious  difficulties.  2D  Euler 
cascade  flows  (Bassi,  1991)  and  3D  Euler  rotor/stator 
interactions  (Trepanier,  1993)  have  been  performed  with 
adaptative  unstructured  meshes  (see  also  Kwon  (1993)  for  3D 
unstructured  Euler  calculations  on  turbine  blades).  Dawes 
(1992,  1993)  presents  complex  3D  turbulent  calculations  for 
different  turbomachinery  configurations  with  adaptative  mesh 
refinement.  In  particular,  a  very  accurate  shock/boundary 
layer  interaction  can  be  represented  after  different  levels  of 
refinement  starting  from  a  coarse  grid.  However,  tetrahedral 


meshes  are  less  tolerant  of  the  high  cell  aspect  ratios  which 
generally  occur  near  solid  boundaries. 

1.3.4  Hybrid  structured/unstructured  meshes 

The  hybrid  approach,  which  combines  structured  and 
unstructured  meshes  (and  solvers)  is  a  way  to  make  the  best 
of  both  approaches.  For  example,  hybrid  calculations  have 
been  used  by  Nakahashi  (1987)  for  2D  and  3D  Navier-Stokes 
turbine  cascade  simulations,  and  also  by  Mathur  (1993)  for 
2D  unsteady  Navier-Stokes  rotor/stator  interaction.  In  both 
cases,  the  authors  use  a  structured  approach  around  the  blade 
to  get  a  regular  and  orthogonal  'mesh  in  the  boundary  layer, 
and  an  unstructured  mesh  in  the  outer  part  which  allows  a 
good  mesh  quality  to  be  generated  in  the  far  field  and  near 
the  periodic  boundary.  Care  must  again  be  taken  at  the 
interfaces. 

1.4  NUMERICAL  SCHEMES 

The  numerical  schemes  used  for  solving  the  Euler  or  the 
Navier-Stokes  equations  for  turbomachinery  flow  simulations 
do  not  differ  from  those  used  for  other  applications.  As  is 
well-known,  a  large  variety  of  numerical  schemes  have  been 
developed  within  the  CFD  community,  and  it  is  very  difficult 
to  list  them  exhaustively.  Therefore,  this  section  gives  a  short 
summary  of  the  most  popular  of  them  for  turbomachinery 
applications.  For  a  detailed  review  of  numerical  schemes 
used  in  fluid  dynamics,  refer  to  Hirsch  (1988,  1990). 

Most  of  the  computations  performed  up  to  now  assume  steady 
flow,  and  most  of  the  methods  used  for  these  calculations  rely 
on  an  unsteady  or  pseudo-unsteady  approach.  In  such 
approaches,  different  combinations  of  space  discretisation 
(centred,  upwind)  and  time  discretisation  (explicit,  implicit) 
exist,  sometimes  based  on  a  coupling  between  space  and  time 
and  sometimes  not,  which  may  use  convergence  acceleration 
techniques.  However,  some  methods  will  also  be  mentioned 
for  solving  the  equations  written  in  steady  formulation. 

1.4.1  Space  discretisation 

The  discretisation  of  the  gradients  is  mainly  performed 
through  finite  volume  methods  or  finite  difference  methods 
on  structured  meshes.  However,  on  structured,  unstructured 
or  hybrid  meshes,  finite  volume  methods  and  finite  element 
methods  have  been  used  (Holmes,  1988  and  Ivanov,  1993). 
Because  of  the  importance  of  mass  flow  conservation  in 
turbomachinery  applications,  finite  volume  methods  are 
generally  preferred. 

1 .4 . 1 . 1  Space  Centred  Schemes 

The  Lax-Wendroff  scheme  (Lax,  1964)  is  based  on  a  Taylor 
expansion  of  second  order  in  time,  where  the  time  derivatives 
are  replaced  by  3-point  space-centred  derivatives  for  the 
Euler  equations.  Several  extensions  of  this  scheme  have  been 
done  in  multi-space-dimensions,  for  the  Euler  and  the  Navier- 
Stokes  equations,  differing  by  their  non-linear  properties  and 
their  grid  dependence.  For  more  details  concerning  these 
properties,  see  the  complete  study  by  Lerat  (1979). 

One  original  formulation  of  the  Lax-Wendroff  scheme 
associated  with  an  efficient  multigrid  method  has  been 
proposed  by  Ni  for  Euler  and  Navier-Stokes  turbomachinery 
simulation  (Ni,  1982,  Davis,  1987),  and  used  by  several 
authors  (see  for  instance  Giles  (1988),  Cambier  (1988),  and 
Heider  (1993)).  This  formulation  is  very  easy  to  implement 
but  does  not  preserve  the  good  numerical  properties  of  the 
original  Lax-Wendroff  scheme  extended  on  a  curvilinear 
mesh  by  Lera  (1979). 
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The  MacCormack  scheme  (MacCormack,  1969),  which  is 
also  second  order  accurate  in  space  and  time,  is  based  on  a 
predictor-corrector  approach  of  the  Lax-Wendroff  type, 
allowing  the  calculation  of  the  Jacobian  matrix  for  tire 
evaluation  of  the  second  order  term  to  be  removed.  It  has 
been  extensively  applied  to  different  external  and  internal 
configurations.  For  turbomachinery  flow  simulations,  see  for 
instance  the  various  single  and  multistage  configurations 
detailed  by  Veuillot  (1985)  and  Fourmaux  (1987),  and  more 
recently  it  has  been  used  for  a  coupled  unsteady 
inviscid/boundary  layer  calculation  by  Tran  (1992). 

It  must  be  noted  that  in  the  above  schemes  coupling  space 
and  time  discretisations  (for  such  types  of  scheme  see  also 
Denton,  1982),  the  numerical  solution  of  a  steady-state 
problem  is  not  independent  of  the  time-step  used  to  reach  this 
steady  state,  which  can  be  seen  as  a  drawback  or  as  an 
advantage  to  build  Lax-Wendroff  schemes  without  need  of 
artificial  viscosity  (Lerat,  1988). 

The  oilier  important  class  of  space  centred  scheme,  used  in 
particular  for  turbomachinery  applications,  is  based  on  the 
scheme  proposed  originally  by  Jameson  et  al  (1981),  built  on 
a  simple  centred  discretisation  of  the  physical  flux  completed 
by  a  non-linear  second  order  dissipation  accompanied  by  a 
linear  fourth  order  dissipation.  Different  versions  of  this 
scheme,  which  is  in  fact  tire  most  popular  scheme  used  in  the 
CFD  community,  have  been  implemented  in  3D  codes  using  a 
node  centred  approach  (Chima,  1991,  Dawes,  1992)  or  a  cell 
centred  approach  (Amone,  1993,  Calvert,  1997,  Dawes, 
1986,  Jennions,  1993,  Kang,  1993,  Liu  1993,  McNulty,  1994, 
Vuillot,  1993,  Shabbir,  1996,  and  Zimmermann,  1992),  with 
various  implementations  of  the  original  numerical  dissipation 
adapted  to  the  type  of  mesh. 

1 .4. 1 .2  Upwind  schemes 

The  centred  discretisations  require  the  addition  of  non-linear 
numerical  dissipation  to  remove  oscillations  near  the  shocks, 
and  can  also  require  the  addition  of  linear  numerical 
dissipation  to  ensure  the  stability  of  the  scheme.  On  the  other 
hand,  the  upwind  approach  uses  directly  the  information 
linked  to  the  hyperbolic  nature  of  the  Euler  equations  which 
is  propagated  along  the  characteristic  lines.  This  upwind 
approach,  which  is  certainly  sounder  theoretically  than  the 
centred  approach,  has  generated  different  classes  of  scheme. 

Godunov's  method  considers  a  piecewise-constant 
approximation  of  the  field  in  each  cell,  associated  with  the 
resolution  of  a  Riemaim  problem  at  each  interface.  Ivanov  et 
al  (1993),  for  turbomachinery  applications,  use  a  piecewise 
cubic  distribution  in  each  cell,  leading  to  a  third  order 
scheme. 

Schemes  based  on  Flux  Vector  Splitting  (Steger,  1981,  Van 
Leer,  1982),  Flux  Difference  Splitting  (Roe,  1981,  Osher, 
1982),  and  Total  Variation  Diminishing  (Harten,  1983,  Yee, 
1987)  have  been  constructed  with  different  space  order 
approximations.  Various  turbomachinery  applications  have 
been  performed  with  such  schemes,  for  instance  by  Allmaras 
(1992),  Benetschik  (1992),  Chen  (1993),  Domey  (1991), 
Engel  (1994),  and  Madavan  (1989,  1993). 

1.4.2  Other  schemes 

Some  authors  use  the  steady  formulation  of  the  governing 
equations  to  solve  steady  problems.  Patankar  and  Spalding 
(1972)  introduced  a  pressure  correction  method  for  solving 
the  incompressible  equations,  extended  by  Hall  (1984,  1992) 
for  the  simulation  of  3D  viscous  compressible  and  transonic 


flows  in  turbomachinery  configurations.  Another  example  of 
such  a  scheme  based  on  a  steady  approach  is  that  described 
by  Moore  and  Moore  (1991)  for  the  calculation  of  the  flow 
through  a  turbine  with  tip  clearance. 

1.4.3  Time  Discretisation 

The  time  discretisation  methods  used  to  solve  the 
Euler/Navier-Stokes  equations  can  be  split  into  two  classes: 
coupled  space-time  methods,  like  the  Lax-Wendroff  or 
MacCormack  schemes,  and  separate  space  and  time  methods. 
For  tire  methods  belonging  to  the  second  class,  the  time 
integration  is  commonly  performed  by  using  a  Runge-Kutta 
approach.  This  Runge-Kutta  time  integration  was  first  used 
for  tlie  solution  of  Euler  equations  by  Jameson  et  al  (1981), 
associated  widi  a  space  centred  discretisation,  but  it  is  also 
used  widi  upwind  space  discretisation.  All  these  schemes 
were  initially  developed  in  an  explicit  approach,  and  for 
turbomachinery  simulations  some  examples  of  2D  or  3D 
calculations  have  been  presented  by  Engel  (1994),  Couaillier 
(1991),  Denton  ( 1 982),  and  Kunz  ( 1 992). 

For  steady  state  problems  based  on  unsteady  approaches, 
implicit  schemes  allowing  the  use  of  large  time  steps  can  be 
very  efficient  for  reaching  the  converged  solution  rapidly.  To 
remove  the  time  step  limitation  of  the  explicit  schemes 
corresponding  to  a  condition  on  die  Courant-Friedrichs-Lewy 
(CFL)  number,  various  studies  have  been  performed  since  the 
contribution  of  Beam  and  Warming  (1978).  For  the  Lax- 
Wendroff  scheme,  Lerat  (1982)  proposed  implicit  methods 
with  ADI  approaches  using  a  block  technique  or  a  spectral 
radius  technique,  leading  to  an  unconditionally  stable  second 
order  scheme  for  a  linear  problem.  This  spectral  radius 
technique,  which  simply  requires  die  resolution  of  a 
tridiagonal  system,  has  been  coupled  widi  Ni's  scheme  (see 
for  instance  Heider,  1993),  and  has  also  been  coupled  by 
Jameson  (1985)  to  die  Runge-Kutta  scheme.  This  approach, 
which  is  well-known  as  Implicit  Residual  Smoothing  (IRS), 
has  been  widely  used  for  ten  years,  for  either  steady  or 
unsteady  problems,  on  structured  meshes  (see  for  instance 
Chima,  1993,  Amone,  1993,  and  Liamis,  1994),  and  also 
adapted  for  unstructured  meshes  widi  an  unfactored  approach 
(Dawes,  1992). 

Rai  et  al.  (1989)  presented  a  diird  order  space  accurate 
upwind  scheme  associated  with  a  fully  implicit  algoridun, 
solved  by  a  subiteration  procedure  for  die  simulation  of 
unsteady  viscous  flow  (see  also  Madavan,  1993).  In  diis 
scheme,  factorisation  and  linearisation  errors  can  be  driven  to 
zero  at  each  time  step.  Such  types  of  subiteration  procedure 
associated  widi  upwind  schemes  are  presented,  for  instance, 
by  Benetschick  (1992)  for  steady  simulations  and  by 
Furukawa  (1992)  for  unsteady  simulations. 

1.4.4  Acceleration  methods 

Different  techniques  can  be  used  to  accelerate  the 
convergence  to  steady  state,  such  as  local  time-stepping  and 
implicit  residual  smootiiing.  Implicit  methods  aim  to  reduce 
die  computational  costs. 

Perhaps  die  most  effective  technique  leading  to  convergence 
acceleration  is  based  on  die  use  of  multigrid  methods.  The 
multigrid  methods  consist  of  transferring  residuals  on  to 
coarse  grids,  solving  die  modified  system  on  diese  grids,  and 
interpolating  die  new  residual  back  to  the  fine  grid.  They 
have  die  property  of  smoothing  die  long  wave  errors  much 
faster  dian  by  using  only  die  fine  grid.  These  mediods  were 
first  developed  at  the  beginning  of  die  1970s  for  the  solution 
of  elliptic  problems,  and  were  later  extended  to  die  solution 
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of  Euler  and  Navier-Stokes  time-dependent  problems  by  Ni 
(1982)  and  Davis  (1987),  associated  with  the  Lax-Wendroff 
scheme,  and  by  Jameson  et  al  (1985)  associated  with  the 
space-centred  Runge-Kutla  scheme.  Jameson's  multigrid 
technique  combined  with  the  IRS  stage  enables  very  efficient 
solvers  to  be  written,  as  for  instance  by  Amone  (1993),  where 
a  3D  Navier-Stokes  calculation  on  a  380,000  mesh  point  grid 
is  presented.  This  calculation  achieves  a  five  order  decrease 
of  the  residuals  after  200  cycles,  corresponding  to  45  minutes 
on  a  CRAY-YMP. 

1.5  TURBULENCE  MODELLING 
In  order  to  close  the  Reynolds-averaged  Navier-Stokes  model 
(RANS)  formulated  with  the  system  of  equations  (3)  and  (4), 
it  is  necessary  to  model  the  Reynolds  stress  T,  and  the 
turbulent  heat  flux  q, .  There  exist  two  broad  classes  of 
turbulence  model,  used  at  the  engineering  level. 

The  first  one,  representing  most  of  the  turbulence  models 
used  in  practice,  is  defined  by  the  eddy  viscosity  concept. 
With  an  eddy-viscosity  concept,  also  called  tire  Boussinesq 
hypothesis,  the  tensor  x,  and  the  flux  q,  are  expressed  by 

relations  analogous  to  those  existing  for  the  laminar  tensor  x 
and  heat  flux  q  : 

f ,  =  -2  /  3 (pk  +  p,V.F)7 +2  p,D 
q,  =  -Cpp,  /Pr,  VI 

Hie  modelling  problem  is  thereby  reduced  to  the 
determination  of  two  scalar  quantities,  the  turbulent  viscosity 
p,  and  the  turbulent  kinetic  energy  k.  Among  the  models 
adopting  the  Boussinesq  hypothesis,  there  are  at  least  two 
distinct  categories : 

Hie  first  approach  consists  of  algebraic  turbulence  models,  in 
which  eddy  viscosity  is  expressed  directly  from  the 
characteristic  properties  of  the  boundary  layer  by  considering 
the  local  equilibrium  of  turbulence  (production=dissipation). 
In  this  approach,  the  turbulent  kinetic  energy  is  generally 
assumed  to  be  zero,  since  it  cannot  be  correctly  estimated. 

The  second  approach  is  transport  equation  models,  where  the 
eddy  viscosity  is  generally  related  to  the  turbulent  kinetic 
energy  and  to  a  characteristic  length  scale  of  turbulence, 
expressed  in  terms  of  the  local  turbulent  quantities. 

Hie  second  class  of  turbulence  model,  based  on  Reynolds 
Stress  Models  and  Algebraic  Stress  Models,  does  not  use  the 
Boussinesq  hypothesis.  Whereas  the  RSM  models  solve  a 
transport  equation  for  each  of  the  stresses,  the  ASM  Models 
use  only  two  transport  equations,  completed  by  algebraic 
relations  for  the  Reynolds  stresses.  Very  few  turbomachinery 
simulations  have  been  performed  up  to  now  with  RSM  or 
ASM  models.  For  example,  see  the  3D  calculations 
performed  by  Hall  (1984)  with  an  ASM  model  associated 
with  wall  functions  for  various  conditions  in  a  turbine  blade 
configuration.  No  model  belonging  to  this  class  has  been  used 
in  the  WG26  calculations. 

Blade  surface  boundary  layers  undergo  transition  from  the 
laminar  to  the  turbulent  state,  in  some  cases  by  “natural 
transition”,  and  in  some  cases  before  a  separated  laminar 
boundary  layer  reattaches  (after  a  “transition  bubble”).  It  is 
also  possible  for  turbulent  boundary  layers  to  revert  to  a 
laminar  state  in  a  strong  favourable  pressure  gradient.  Hie 
prediction  of  these  effects  is  critical  to  the  accurate 
estimation  of  surface  heat  transfer,  and  in  some  cases  it  is 
also  important  to  loss  prediction. 


hi  a  turbomachine,  rapid  transition  has  been  shown  to  be 
forced  in  some  cases  by  the  impact  of  a  shock  wave  or  a  wake 
arriving  from  an  upstream  blade  row,  followed  by  reversion 
to  the  laminar  state  when  the  disturbance  has  passed. 
Transition  is  known  to  depend  on  the  free  stream  turbulence 
and  also  on  the  roughness  of  the  surface.  So  it  is  not 
surprising  that  there  is  as  yet  no  generally  accepted  way  of 
predicting  transition  within  a  RANS  computation,  although 
some  of  the  turbulence  models  used  include  a  way  of 
estimating  where  transition  is  likely  to  occur.  Mayle  (1997) 
has  recently  summarised  transition  prediction  methods  for 
two-dimensional  boundary  layers,  but  three-dimensional 
transition  is  obviously  more  difficult  to  predict. 

Almost  all  the  computations  undertaken  for  the  WG  have 
assumed  the  blade  surface  boundary  layers  to  be  fully 
turbulent  from  the  leading  edge.  Hris  assumption  is  believed 
to  be  adequate  for  the  test  cases  chosen  by  the  WG,  but  in 
general  it  is  likely  to  be  a  significant  limitation. 

1.5.1  Mixing-length  algebraic  models 
a)  Most  of  the  Navier-Stokes  simulations  for  turbomachinery 
are  still  carried  out  with  simple  algebraic  turbulence  models. 
Hie  most  popular  of  them  is  the  Baldwin-Lomax  model 
(Baldwin,  1978)  derived  from  the  Cebeci-Smith  model 
(Cebeci,  1974),  and  a  large  number  of  turbomachinery 
computations  have  been  done  with  it  (see  for  instance  Cliima, 
1993,  Davis,  1987,  Engel,  1994,  Liu,  1988,  Nakahashi,  1987, 
and  Scott,  1986). 

This  model  considers  the  boundary  layer  split  into  two  parts: 
an  inner  region  denoted  by  the  subscript  i  and  the  outer 
region  denoted  by  the  subscript  e.  Hie  turbulent  viscosity  p, 
is  defined  as  follows : 

p,  =  p„  if  d  <,  d0  and  p,  =  ptt  if  d  £  d0 

where  d  is  the  distance  to  the  wall  and  d0  is  determined  as 
the  smallest  value  of  d  for  which  p„  =  p,e . 

hi  the  inner  region  the  mixing  length  is  defined  by:  1  =  Kd , 
where  K  =  0.41  is  the  Von  Karman  constant.  Hie  eddy 
viscosity  p«  is  then  given  by  the  relation:  p„  =  pl2D7  |tn| , 
where  m  is  the  vorticity.  Hie  damping  function  D  is  defined 
as  follows : 

D(d)  =  l-exp(-d+  IA+) 

where  A*  =  26.  The  non-dimensional  distance  ct  is  defined 
using  the  skin  velocity  Kt(xw  =  pwV2)  and  the  density  at  the 
wall  pw  and  molecular  viscosity  p d+  =  p wVtd  /  pH, . 

In  the  outer  region,  the  eddy  viscosity  \i,c  is  given  by  the 
relation : 

P«  =  ^CcppFwakeFKieb{d) 

where  C;  =  0.0168  and  Ccp  =  1.6.  Hie  quantity  Fwau  is 
defined  as  follows : 

Fwake  ~  ^inax  mill(/~Imx  ,  Cwk  Vm  /  7*ntax  ) 
where  Cwk  was  originally  given  the  value  0.25,  dmax  is  the 
value  where  the  function  F(d)  =  d.D(d).w  attains  its 
maximum  value  Fmax,  and  where  Vm  denotes  the  difference 
between  the  maximum  and  the  minimum  velocity  inside  the 
shear  layer.  Hie  Klebanoff  intermittency  factor  Fiaeb  is 
defined  by : 

Fk>'„  =  [l  +  ^Mf-Klebd  /  dnm  )6  ]  ' 
with  Ckm  =  0.3. 
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b)  In  the  model  of  Michel  et  al  (1969),  the  turbulent  viscosity 
p,  is  obtained  as  follows : 

p,  =  p/2F2jur| ,  where  /  =  0.085  8  tanh  (Kd  /  0.0855) . 

where  F  is  the  Van  Driest  viscous  sublayer  damping  function 
given  by : 

F(£)  =  1  -  exp(-^ / 26K) ,  where  £  =  p/2(p  +  p,)/p2|®| . 

Hiese  relations  lead  to  an  implicit  equation  for  p,,  which  is 
solved  by  Newton’s  method.  The  quantities  d  and  8  are 
respectively  the  distance  to  the  wall  and  the  boundary  layer 
thickness. 

Theoretical  and  practical  difficulties  appear  when 
implementing  tliese  models  in  the  presence  of  several  walls, 
where  comer  flows  exist.  In  that  case,  if  the  quantity  d  can  be 
relatively  easily  defined  for  both  models,  the  evaluation  of 
tlie  quantities  dmax  and  Fmax  for  the  Baldwin-Lomax  model 
and  8  for  tire  Michel  model  causes  great  difficulties. 

1.5.2  One-equation  Spalart-Allmaras  model 
Recently  there  lias  been  significant  interest  in  one  equation 
transport  models  as  a  way  of  obtaining  the  advantages  of  a 
transport  model  at  a  minimum  computational  cost.  Perhaps 
the  most  popular  at  present  is  the  Spalart-Allmaras  model 
(Spalart  and  Alhnaras,  1992,  1994)  which  solves  an  equation 
for  the  high  Reynolds  number  eddy  viscosity  v  written  in  the 
following  form : 

^  =  CM[1  -  fnfw  +  -k[v.((v  +  v)Vv)  +  Ci2(W)2 )] 

~  +fnAU2 


where  the  turbulent  kinematic  viscosity  v,  is  given  by 


v,  =  v./v, ,  where  /„,  = 


x3+cvl 


j  with  %  =  - 


S  =  S  +  £2(j2  A i  w'iere  S  is  tlie  vorticity, 


A  2=1- 


1+X/vl’ 


-  l+C^  .  n  /  6  \  _v 

fw  ~  S  6  6  9  g-r+C^f*  ~r)>  r  “  rrr2  i2 

18  +C*3j  SK  d 

fa  =  C/3  exp(-Ct4  %2) 

/,1=C„g,exp(-Ct2-^Tp2+g,2c/,2] 

g,  =  inin(0.1,At/  /  w>,Ax,) 

Ax,  =  grid  spacing  along  the  wall 

at  the  nominated  transition  point 
d,  =  distance  between  the  local  point  and  the  transition  point 
w,  =  vorticity  at  the  wall  at  the  nominated  transition  point 

values  of  constants: 

Cbl  =  0.1355 , 0  =  2  /  3  ,  Cb2  =  0.622  ,  K  =  0.41 
Cwi  =  Q|  /  K  +  (1  +  Ctj)/  ct  ,  Cw2  =  0.3  ,  Cw3  =  2 
Cv i  =  7.1 ,  C(|  =  1 ,  C(2  =  2  ,  Cy  =  1.1 ,  Ct4  =  2 

(hi  Spalart  and  Allmaras  (1994),  C, 3  =  1.2,  C, 4  =  0.5.) 

Tlie  boundary  condition  on  tlie  wall  is  expressed  by  setting 


1.5.3  Two-equation  turbulence  models 
Because  they  provide  a  good  compromise  between  accuracy 
and  computational  cost,  the  use  of  turbulence  models  with 
two  transport  equations  for  turbomachinery  applications  has 
been  increasing  for  several  years.  Examples  of  these  models 
are  tlie  Jones-Launder  k-e  model  (Jones,  1973),  tlie  Wilcox  k- 
co  model,  tlie  Wilcox-Rubesin  k-co2  model  (Wilcox  and 
Rubesin,  1980,  Wilcox,  1991)  and  tlie  Coakley  q-®  model 
(Coakley,  1983).  The  low-Reynolds-number  Jones-Launder 
model  and  Wilcox  model,  which  are  in  fact  popular  two- 
equation  models,  are  presented  hereafter.  For  these  two 
models,  the  low  Reynolds  number  closure  coefficients  do  not 
depend  on  geometrical  quantities,  such  as  the  reduced  wall 
distance  y+  depending  on  the  wall  distance  and  on  the  friction 
velocity  inaccurately  defined  in  separation  zones,  but  instead 
tlie  damping  functions  are  based  on  the  turbulent  Reynolds 
number  Re,,  only  depending  on  local  flow  quantities: 

Re,=ril,A 

\iZ  p® 

(Note  that  this  equation  defines  the  symbol  ®  ,  which  is  not 
the  same  as  tlie  vorticity  m  used  in  Section  1.5.1.)  Hence, 
from  a  practical  point  of  view,  these  two-equation  models  are 
well-suited  for  complex  three-dimensional  configurations. 
However,  there  can  be  numerical  difficulties  associated  with 
solving  two-equation  models  with  explicit  schemes  due  to  the 
stiffness  of  the  equation  system  (Dailey  et  al,  1994). 

a)  The  equations  of  the  k-e  system  can  be  written  as  follows : 

— +  V.(pfrF)  =  V.  +?.:VF-p e-D 

dt  ctkJ  J 

^£  +  V.(pFe)  =  V.  [p  +  — |Ve  +C.f,-T ,:VF 

8f  aj  J  k 

-C2f2~  +  E 

k 

In  these  equations,  D  is  a  wall  term  and  E  is  a  low  Reynolds 
term  introduced  to  solve  tlie  k-e  system  in  the  region  down  to 
tlie  wall  with  low  turbulent  Reynolds  number  (Re,  £  100) . 
Hie  turbulent  viscosity  is  evaluated  by  the  relation: 
p(  =  c^pk2  /  pe .  Some  examples  of  the  terms 
E,f^ ,/,  and  f2 ,  including  empirical  constants  and  damping 

functions,  were  presented  by  Patel  (1984).  Hie  Jones- 
Launder  model  is  based  on  the  following  definitions  : 

D  =  2p(vV*)2,  E  = 

P  l  dn2  J 

ak  =Lae  =  13,  Cj  =  157, C2  =2, 
f  -  I./2  =  [1  -0.3exp(-Re2 )], 

f„  =  expf - — — |,C„  =  0.09 

Al  +  Re,/50j 

Hie  boundary  condition  on  the  wall  is  expressed  by  setting 
k= 0  and  s=0. 


b)  The  equations  of  the  k-®  system  can  be  written  as  follows  : 
+  V.(pFA-)  =  V  j(p  +  ct*p,)va]  +  T ,:  VF  - p'pA® 


■^^-  + V.(pF®)  =  v|(p  +  CTp,)V®j  +  ay-x,:VF  -Pp®2 

The  turbulent  viscosity  is  evaluated  by  tlie  relation: 
p,  =a*pA/ro.  The  parameters  a,a’,p,p”,a,a”  in  the 
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above  equations  are  the  closure  coefficients  of  the  model  and, 
according  to  Wilcox  (1994),  are  given  by  : 


a  =  a, 


1  +  Re,  /  Ra  }  ’ 

p=p„  p*=p; 


.  .  +Re./  Rk 

a  =a,  — - — — - — - 
'  1  +  Re,/F* 

p;+(Re,/^)4 

1  +  (Re(/  Rp)4 


a  =  a’  =1/2, 

a0  =  1/10, aj  =P,  / 3, Pq  =5/18, 
a,  =5/ 9,«;  =  1,P,  =  3/40,p*  =9/100, 
■Rp  =8,R*  =6,/^  =2.7 . 


At  solid  walls,  £=0  and 
u  2 

to  =  -^—SR  where  wt  =  friction  velocity  = 
v 

Sr  is  related  to  the  non-dimensional  surface  roughness 
k  u 

kp  =  -  1  tlirough  the  correlation 
/  \2 

for  kx  <  25 

for  k%  S  25 
where  k£  £5  for  smooth  surfaces. 


The  k-<o  model  presented  here  includes  low  Reynolds  number 
modified  closure  coefficients,  which  aim  to  simulate 
transition.  This  form  of  the  coefficients  is  not  strictly 
necessary  in  a  turbulent  boundary  layer  even  if  .the  equations 
are  integrated  down  to  the  wall. 


1.5.4  Wall  functions 

1. 5.4.1  Wall  function  corrections  for  skin  friction 
The  calculation  of  the  viscous  shear  stresses  at  solid  surfaces 
is  a  critical  factor  in  predicting  boundary  layer  development 
accurately.  A  simple  calculation  is  adequate  if  the  near  wall 
grid  spacing  is  small,  i.e.  xw  =  p(Vp  -Vw)/Yp  ,  where  Vp  is 

the  flow  velocity  at  the  point  nearest  to  tire  wall,  Tw  is  the 
distance  of  the  point  from  the  wall,  and  V„  is  the  velocity  of 
the  wall. 

For  equilibrium  turbulent  boundary  layers,  this  gives  accurate 
results  if  the  point  nearest  the  wall  is  in  the  laminar  sub¬ 
layer,  (i.e  y+  <  5),  but  it  underestimates  tw  for  coarse  grids. 
For  example,  assuming  Spalding’s  expression  for  the  law  of 
the  wall,  the  error  is  just  over  50%  if  the  first  point  is  at  y+  = 
25.  It  is  then  necessary  to  apply  corrections  based  on  the  law 
of  the  wall  to  give  a  closer  estimate  to  tire  velocity  gradient  at 
the  wall.  The  corrections  can  be  correlated  against  p(Vp- 
V„)Ypl\i  in  order  to  avoid  tire  expense  of  solving  the  law  of 
tire  wall  during  the  calculation. 


1. 5.4.2  Wall  functions  for  high  Reynolds  number  turbulence 
models 

hi  order  to  avoid  the  use  of  grid  nodes  very  near  to  the  wall 
(y+<5),  wall  functions  are  often  used.  Let  Ay  be  the  size  of 
the  nearest  cell  to  the  wall,  V  be  the  velocity  V  at  Ay  ,  and  © 
-  (Tav/2-Tj-o  ),  where  T  =  temperature.  Then 


pVA y 


-  1£ 


p  p  2pK’  2  \xcp® 

The  boundary  conditions  for  k  and  e  at  the  wall  are: 


KAy 


hi  these  equations,  Pr  is  the  Prandtl  number,  qp  is  the  wall 
heat  flux,  xp  is  the  wall  shear  stress,  and  u,  is  the  friction 
velocity. 


According  to  Shabbir  et  al  (1997),  the  following  relations  due 
to  Shih  and  Lmnley  (1993)  must  be  used  for  y+>l  1.  They  also 
allow  the  integration  of  the  equations  up  to  the  wall. 

4 

k  =  025u.1,  e  =  0.25— *— 

Vlum 

The  friction  velocity  ux  (or  the  wall  shear  stress  xp)  is 
deduced  from  the  coefficient  ap.  Similarly,  the  wall  heat  flux 
is  deduced  from  the  coefficient  pp..  One  of  two  models  could 
be  used,  or  a  blending  of  tire  two  models  to  form  a  third 
model: 


1st  model 

Assume  a  Couette  flow  with  turbulence  in  equilibrium  near 
the  wall: 

alp  =  l+Ft,  Pip  =  1  +F,G 


where 


2nd  model 

The  turbulence  is  produced  outside  the  wall  boundary  layer, 
and  is  diffused  towards  the  wall: 
a2p  =  l+F2,  p2p  =  l+F2G 

where 


F,  =  expl  — j  + 

2  F|  20  J  3.05  +  ln(l  +  v+) 


-1,  v;=-!- 


3  rd  model 

a  blending  of  the  first  two  models  (Duchene,  1995): 
ot3;,  =  l  +  F3,  P3p  =  1  +  FjG 
F3  =  max(F|,F2) 


1.5.5  Validation  of  turbulence  models 
The  use  of  a  reliable  turbulence  model  being  an  important 
key  to  capturing  the  main  features  of  complex  internal  flows, 
and  then  to  predict  the  losses,  several  studies  comparing 
different  turbulence  models  according  to  different  levels  of 
grid  refinement  have  been  done. 

Cliima  et  al  (1993)  presented  a  modified  Baldwin-Lomax 
model  showing  good  agreement  with  experimental 
measurements  of  kinetic  energy  efficiency  for  an  annular 
turbine  cascade,  and  of  heat  transfer  at  the  endwall  for  a 
linear  turbine  cascade. 

Several  other  authors  have  presented  results  with  a  modified 
Baldwin-Lomax  model  (see  for  instance  Granville,  1987  and 
York,  1985).  The  “2D  ONERA  bump”  is  a  typical 
configuration  used  to  validate  turbulence  models  for 
shock/boundary  layer  interaction.  This  test  case  has  been 
used  in  the  EUROVAL  project  (Haase,  1993)  and  more 
recently  for  an  ETMA  Workshop  organized  by  UMIST 
(ETMA,  1994).  In  these  references,  several  calculations 
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performed  with  a  mixing-length  model,  a  two-transport 
equation  model,  or  an  ASM  or  RSM  model,  using  space 
centred  or  upwind  discretisations,  are  detailed,  hi  particular, 
systematic  comparisons  of  mean  and  turbulent  quantities  are 
presented.  From  these  results  it  appears  that,  for  methods 
using  second  order  space  accuracy,  the  influence  of  the 
turbulence  model  is  greater  than  the  influence  of  the  space 
discretisation  if  sufficient  mesh  refinement  is  used.  The 
Baldwin-Lomax  and  Michel  models  do  not  predict  the 
characteristic  X  shock  interacting  with  the  viscous  layer.  The 
k-s  model  provides  a  better  representation  of  the  shock 
behaviour  but  underpredicts  the  separation  zone,  whereas  a 
multi-scale  model  (Gleize,  1994)  is  as  close  to  the 
experiment  as  an  RSM  simulation  presented  in  (Lien,  1993). 

Biswas  (1994)  gave  some  examples  of  the  different  k-s  model 
versions  and  compared  them  on  a  flat  plate  test  case  (Savill, 
1993)  and  on  a  turbine  blade.  A  detailed  study  of  the  flow  in 
a  transonic  fan,  modelled  by  different  Baldwin-Lomax  and  k- 
s  versions,  was  presented  by  Jennions  and  Turner  (1993).  In 
this  paper,  it  is  shown  how  an  extended  k-e  model,  which 
takes  into  account  multiple  time  scales,  provides  a  better 
prediction  of  the  shock  position. 

A  3D  shock/boundary  layer  interaction  in  a  transonic  channel 
with  a  swept  bump  has  been  measured  and  compared  with 
simulations  performed  using  the  Michel  model  and  using  the 
k-s  model  (Calien,  1993).  The  results  show,  as  for  2D  cases, 
a  better  prediction  using  the  k-s  model  in  the  shock/boundary 
layer  interaction  region,  as  in  the  comer  regions,  which  is  an 
interesting  point  for  3D  turbomachinery  applications.  A 
comparison  of  skin-friction  pattern  lines  between  calculation 
and  experiment,  presented  for  the  wall  bearing  the  bump  and 
for  the  adjacent  wall  with  a  strong  interaction,  shows  the 
interesting  behaviour  of  the  predicted  solution.  Nevertheless, 
even  if  the  main  tendencies  are  well  represented,  the  k-s 
calculation  does  not  predict  correctly  the  levels  of  kinetic 
energy  and  the  length  of  tire  separation  zone  in  the  interaction 
region.  A  comparison  between  a  Baldwin-Lomax  model  and  a 
k-E  model  for  a  3D  subsonic  blade  row  configuration  was 
presented  by  Matsuo  (1991).  As  above,  it  shows  also  that  the 
k-s  model  leads  to  a  better  prediction  of  the  separated  region, 
even  if  the  secondary  flows  are  not  correctly  reproduced. 

Fougeres  (1994)  undertook  numerical  simulations  of  heat 
transfer  with  film  cooling  using  the  Michel  model  for  a  flat 
plate  test  case  and  for  a  nozzle  guide  vane  configuration. 
These  calculations  were  done  by  using  overlapping  grids  to 
refine  the  fluid  injection  holes  accurately.  It  is  demonstrated 
that  a  good  quality  fine  mesh  is  required  to  obtain  well- 
predicted  heat  transfer  coefficients. 

1.6  CODE  VALIDATION 

As  already  indicated,  3D  steady  viscous  turbomachinery 
applications  are  commonly  computed  today  with  solvers 
which  have  been  “calibrated”  on  typical  test  cases  where 
detailed  experimental  data  are  available.  In  order  to  validate 
these  solvers,  not  only  the  wall  pressure  distribution  is 
needed  but  also  finer  experimental  measurements,  like 
boundary  layer  velocity  profiles,  heat  transfer  coefficients, 
wall  stresses  and  velocity  vectors  in  the  field.  In  this  section 
some  examples  are  presented  of  compressor  and  turbine 
blade  flow  calculations  which  are  typical  of  tire  state  of  the 
art  reached  for  these  configurations,  and  which  have  been 
performed  on  grids  probably  finer  than  those  currently  used 
in  tlie  blade  design  process.  See,  for  instance,  the  proceedings 
of  tire  AGARD  Conference  on  CFD  techniques  (AGARD, 
1991)  and  tire  proceedings  of  tire  seminar  organized  by  tire 


ERCOFTAC  Turbomachinery  Special  Interest  Group 
(Gregory-Smith,  1993)  to  find  various  turbomachinery  flow 
calculations  compared  with  experiment. 

1.6.1  Compressor  flows 

Several  authors  have  tested  their  codes  on  tire  NASA  Lewis 
Rotor  67  fan  documented  by  Fottner(1990),  and  some  of 
these  calculations,  including  tip  clearance,  are  reported  in 
AGARD  (1991).  As  regards  tire  detailed  analysis  of  tip 
leakage  phenomena  in  compressor  configurations,  in  tenns  of 
tire  physics  and  also  tire  influence  of  the  mesh  on  numerical 
solutions,  see  for  instance  Adamczyk  (1993),  Chen  (1991), 
and  Kang  (1993).  For  such  a  transonic  compressor,  different 
complex  flow  features  provide  an  interesting  challenge  to 
validate  numerical  methods:  a  bow  shock,  a  separation  zone 
induced  by  tire  passage  shock/boundary  layer  interaction  on 
tire  suction  surface,  tip-leakage,  and  a  strong  radial  flow,  ... 
Two  flow  configurations,  corresponding  to  the  peak  efficiency 
condition  and  the  near  stall  condition,  are  presented  by  Hah 
(1992)  using  a  k-s  model  on  a  250,000  point  H  grid  and  by 
Amone  (1994)  using  a  Baldwin-Lomax  model  on  a 
350,000  point  H  grid.  In  comparison  with  laser 
measurements,  both  calculations  reproduce  the  main  features 
of  the  shock  structure  at  different  sections  of  the  span. 

Couaillier  (1991)  performed  a  calculation  on  the  SNECMA 
TS27  wide  chord  fan  at  tire  peak  efficiency  condition,  using  a 
617,000  point  H-O-H  grid.  No  tip  clearance  was  taken  into 
account  and  a  slip  condition  was  imposed  on  the  casing.  The 
use  of  an  O  Mesh  around  the  blade  led  to  an  accurate 
representation  of  the  leading  edge  bow  shock,  and 
comparison  with  laser  measurement  data  showed  good 
agreement  with  the  shock  location. 

Nozaki  (1993)  simulated  tire  flow  through  a  complete  stage  of 
a  fan  configuration  for  different  operating  points.  The 
calculations  did  not  take  into  account  unsteadiness  because 
the  boundary  conditions  on  tire  rotor/stator  interface  were  set 
to  circumferentially  averaged  values.  Nevertheless,  tire 
overall  time-averaged  performance  showed  a  qualitatively 
good  agreement  with  the  experimental  measurements. 

The  ASME  test  case  exercise  in  1994  using  Rotor  37  has 
already  been  mentioned  as  one  of  the  starting  points  of  the 
present  WG. 

1.6.2  Turbine  flows 

The  use  of  3D  Navier-Stokes  solvers  gives  the  opportunity  to 
analyse  such  typical  flow  features  as  film  cooling,  heat 
transfer,  and  transition  phenomena.  As  an  example  of  a 
turbine  test  case  with  strong  effects  of  transition  phenomena 
and  secondary  flows,  the  Durham  cascade  (Gregory-Smith, 
1992)  has  been  computed  by  several  groups  and  reported  by 
Gregory-Smith  (1994).  It  appears  that  the  calculated  solutions 
are  largely  dependent  on  the  turbulence  models  and  transition 
criteria.  The  analysis  presented  by  Wegener  (1992)  showed 
comparisons  of  the  wake  development  behind  an  annular 
turbine  stator  predicted  by  two  different  solvers  with  the 
experiment.  A  numerical  study  of  heat  transfer  for  a  2D 
turbine  rotor  cascade  and  for  a  3D  turbine  cascade  near  the 
endwall  was  presented  by  Hah  (1989).  The  3D  computation, 
which  was  performed  with  a  k-e  model  on  a  600,000  grid 
point  mesh,  provided  a  numerical  solution  showing  an 
acceptable  accuracy  for  heat  transfer  prediction.  A  Navier- 
Stokes  analysis  of  turbine  blade  heat  transfer  was  also 
presented  by  Domey  (1991). 

Zimmermann  (1992)  used  the  Baldwin-Lomax  model  to 
undertake  a  numerical  analysis  of  tire  total  pressure  losses  in 
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a  transonic  turbine  cascade.  Comparisons  with  experiment 
show  acceptable  predictions,  and  an  analysis  of  the  pressure 
loss  components  leads  to  the  conclusion  that  the  total 
pressure  losses  depend  strongly  on  secondary  flows.  The 
analysis  presented  by  Horton  (1991)  showed  the  influence  of 
the  mesh  quality  on  total  loss  pressure  prediction  with  two 
endwall  geometries  of  a  transonic  turbine. 

Linear  cascades  and  rotor  turbine  configurations  with  tip 
clearance  were  calculated  by  Fourmaux  (1993)  using  the 
Michel  model.  Both  for  the  linear  cascade  and  the  rotor 
configurations,  a  five-domain  mesh  with  1 1  points  in  the  gap 
between  tire  blade  tip  and  tire  casing  was  used.  When 
comparing  results  obtained  with  and  without  tip  clearance 
representation,  large  differences  in  the  flow  field  were  seen 
in  tire  hub  region,  covering  25  %  of  the  blade  in  the  spanwise 
direction. 
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Chapter  2 

THE  TEST  CASES 


The  two  test  case  data  sets  described  in  this  report  were 
selected  by  the  members  of  the  Working  Group  from  a  larger 
group  of  candidate  sets.  They  are 

(1)  an  isolated  axial-flow  compressor  rotor  designed  and 
studied  experimentally  at  the  National  Aeronautics  and 
Space  Administration  Lewis  Research  Center  (US), 
identified  as  NASA  Rotor  37. 

(2)  an  axial-flow  turbine  inlet  stator  blade  row  designed 
and  studied  experimentally  in  the  Institut  fur 
Antriebstechnik  of  the  Deutsche  Forschungsanstalt  filr 
Luft  -  und  Raumfahrt  e.V.  (GE).  This  configuration  is 
identified  here  as  the  DLR  turbine  stator. 

These  data  sets  were  selected  following  a  thorough  review  of 
the  experimental  programmes  and  the  available  data.  Both 
configurations  were  representative  products  of  proven  design 
methods  and  carefully-controlled  manufacturing  processes. 
The  test  programmes  were  carried  out  in  well-developed  test 
installations  operated  by  experienced  research  groups.  In  both 
programmes  a  wide  range  of  data  were  available,  based  on 
proven  instrumentation  and  data  acquisition  methods,  and 
reduced  by  accepted  and  traceable  procedures.  Finally  both 
test  cases  involved  experiments  under  steady-state  entrance 
flow  conditions,  with  internal  and  exit  flows  representing  a 
substantial  challenge  to  current  CFD  codes. 

In  the  following  sections  the  two  test  cases  are  individually 
discussed,  with  reference  to  design  background,  flow  path 
and  blade  geometry,  test  installation  and  methods,  and 
definition  of  the  data  available  and  relevant  to  the  current 
study. 

2.1  NASA  ROTOR  37 


Equivalent  rotor  tip  speed  U, 


=  454.1  m/s 


Equivalent  mass  flow  per  unit  annular  area 
=  200.5  kg/s/m2 


Aa„  Pt\ 


Jref 


where  Pti  =  inlet  total  pressure 

Pref  =  101.33  kN/m2  (sea  level  standard  atmosphere) 
Aan  =  annulus  area 


Rotor  total  pressure  ratio  =  2.106 
Rotor  polytropic  efficiency  =  0.889 
Number  of  rotor  blades  =  36 


2.1.2  Configuration  Geometry  for  WG26  Test  Case 
Subsequent  to  the  tests  of  Moore  and  Reid  (1980)  on  NASA 
Stage  37,  the  rotor  was  retested  as  an  isolated  component. 
This  is  the  geometry  identified  by  NASA  as  Rotor  37.  Fig.2.1 
shows  in  schematic  form  the  annular  flow  path  and  blade 
airfoil  geometries  with  coordinate  reference  definitions  for 
computational  purposes.  These  coordinate  definitions  are 
those  utilized  in  all  data  reported  and  in  all  performance 
computations  reported  in  other  sections  of  this  document.  All 
of  the  values  are  estimated  design  point  operating  values  of 
dimensions  and  angles.  Detailed  numerical  values  of  passage 
and  blade  coordinates  are  not  given  in  this  document  but  are 
available  by  electronic  data  file  transfer  as  outlined  in  Section 
2.1.5. 

Additional  geometry  of  interest  for  CFD  purposes  is  shown  in 
Fig.  2.1  .  The  blade  hub  fillet  radius  was  2.5  mm  and  the 
RMS  blade  surface  roughness  was  0.5  - 1.25  microns. 


2.1.1  Design  Background 

Rotor  37  was  designed  and  initially  tested  as  part  of  a 
research  programme  involving  four  related  axial-flow 
compressor  stages.  These  stages  were  intended  to  cover  a 
range  of  design  parameters  typical  of  aircraft  turbine  engine 
high-pressure  (core)  compressor  inlet  stages.  In  the  case  of 
stage  37,  representative  values  were: 

Rotor  inlet  hub-to-tip  diameter  ratio  0.7 

Rotor  blade  aspect  ratio  1.19 

Rotor  tip  relative  inlet  Mach  number  1 .48 

Rotor  hub  relative  inlet  Mach  number  1.13 

Rotor  tip  solidity  1.29 

Blade  airfoil  sections  Multiple-Circular-Arc  (MCA) 

No  inlet  guide  vanes  were  specified  for  any  of  the  stages. 
Some  design  information  and  overall  stage  performance 
results  were  reported  by  Reid  and  Moore  (1978).  More 
detailed  stage  performance  was  reported  later  by  Moore  and 
Reid  (1980).  It  should  be  noted  that  while  the  designs  and 
stage  tests  were  initiated  during  the  1970's,  geometries  and 
performance  levels  are  similar  to  those  for  current  turbine 
engine  stages. 


2.1.3  Test  Installation 

Figure  2.2  is  a  schematic  diagram  of  the  single-stage  axial- 
flow  compressor  stage  test  installation  at  the  NASA  Lewis 
Research  Center  in  Cleveland,  Ohio  (US).  This  is  the  facility 
used  for  all  experimental  work  on  the  NASA  Stage  37  and 
Rotor  37  configurations.  Atmospheric  air  enters  the  system 
through  an  inlet  and  filters,  followed  by  a  standard  thin-plate 
orifice  for  compressor  mass  flow  rate  measurement. 
Compressor  inlet  pressure  and  mass-flow  rate  are  set  by 
combined  use  of  parallel  butterfly  valves  in  the  inlet  ducting 
and  a  sleeve  throttle  between  the  test  unit  and  the  collector 
ducting.  Air  enters  the  test  compressor  from  an  upstream 
plenum  tank  containing  screens  and  a  honeycomb  grid.  Flow 
enters  the  annular  flow  path  of  Figure  2.3  through  a  bell- 
mouth  inlet  with  a  central  nose  cone  and  three  airfoil-shaped 
support  struts.  The  installation  exhausts  to  the  atmosphere  or 
to  the  Center  altitude  exhaust  system.  Atmospheric  exhaust 
was  utilized  in  all  Rotor  37  experiments. 

The  test  compressor  is  driven  by  a  3000  hp  alternating 
current  motor  through  a  speed-increasing  gear  box.  Motor 
speed  is  varied  by  changing  supply  power  frequency. 


Design  point  values  for  the  rotor  as  estimated  in  the  design 
computations  were: 

Equivalent  rotational  speed 


=  17188.7  rpm  (1800  rad/s) 


where  Tti  =  inlet  total  temperature 

Tref  =  288.15  K  (sea  level  standard  atmosphere) 


2.1.4  Experimental  Programme 

The  data  sets  selected  for  the  PEP  Working  Group  26  test 
case  were  developed  during  1992  and  1993.  The  same  data 
were  the  basis  for  an  extensive  CFD  code  assessment  effort 
sponsored  by  the  ASME  International  Gas  Turbine  Institute's 
Turbomachinery  Committee  during  1993  and  1994.  Wisler 
(1993)  and  Denton  (1996)  discuss  the  organization  and 
results  of  this  work. 
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Figure  2.1  Blade  and  fiowpath  coordinates 


2. 1 .4. 1  Test  Conditions 

All  of  the  data  used  in  the  AGARD  and  ASME  code 
validation  studies  were  measured  with  NASA  Rotor  37 
operating  at  the  design  equivalent  rotational  speed  defined 
above.  To  establish  reference  points  for  detailed  flow  field 
measurements,  overall  performance  was  determined  at 
equivalent  mass  flow  rates  from  the  maximum  attainable, 
referred  to  in  this  document  as  mchokt ,  to  a  minimum  slightly 
above  the  rotor  stall  flow.  This  near  stall  flow  rate  was 
experimentally  determined  to  be  m  /  mchokt  =  0.925.  A  total  of 
13  sets  of  overall  performance  data  were  measured  over  this 
mass  flow  rate  range.  These  points  are  shown  on  Figure  2.4. 

Data  points  tit / mchokt  =  0.98  and  rii/mfhokt~  0.925  were 
designated  by  both  ASME  and  the  Working  Group  as  flow 
rates  for  concentration  of  CFD  effort.  Tire  experimental 
™ choke  as  determined  by  NASA  was  20.93  kg/s.  All  reported 


values  are  equivalent  values  referenced  to  sea-level  standard 
atmosphere  values  of  Prer  and  T  ref . 

2. 1.4.2  Performance  Measurement 

The  performance  data  acquired  for  the  NASA  Rotor  37  test 
case  include  averaged  overall  total  pressure  ratio  and 
adiabatic  efficiency  for  a  range  of  equivalent  mass  flow  rates 
with  the  rotor  operating  at  the  design  value  of  equivalent 
rotational  speed.  Mass  flow  rates  reported  are  based  on  tire 
test  facility  orifice  plate  measurements.  Rotational  speed  was 
measured  using  a  shaft-mounted  magnetic  pick-up  system  and 
an  electronic  counter. 
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Figure  2.2  Test  facility 


In  addition,  flow  passage  data  based  on  radial  surveys  with 
pneumatic  pressure  probes  and  thermocouples  were  recorded 
at  each  equivalent  mass  flow  rate  at  design  equivalent 
rotational  speed.  The  surveys  were  at  Station  1  and  Station  4 
(Fig.  2.3).  The  probes  used  were  a  combination  probe  with  a 
"cobra"  head  for  total  pressure  and  absolute  flow  angle 
information,  with  a  thermocouple  for  total  temperature  data; 
and  a  wedge  probe  for  static  pressure  data.  Wall  static 
pressure  taps  were  also  located  on  hub  and  tip  walls  at 
Stations  1  and  4.  The  circumferential  locations  of  survey 
positions  are  shown  in  Figure  2.5.  Radial  surveys  included  18 
measurement  locations. 

Turbulence  intensity  was  measured  at  Station  1.  Laser 
anemometer  velocity  data  were  acquired  in  four  hub-to-tip 
measuring  planes,  and  also  on  five  blade-to-blade  surfaces  of 
revolution.  Locations  of  these  planes  and  surfaces  are  defined 
schematically  on  Figure  2.3.  The  laser  anemometer  system 
used  and  the  procedures  followed  have  been  described  in 
detail  by  Suder  et  al  (1995),  Suder  and  Celestina  (1996)  and 
Hathaway  et  al  (1993).  Laser  data  were  acquired  only  at  flow 
rates  of  tit  /  riichokt  =  0.98  and  0.925. 

2. 1.4. 3  Data  Reduction 

For  the  ASME  test  case  exercise,  tire  averaging  scheme  used 
to  calculate  tire  performance  from  the  CFD  results  was 
closely  defined  by  Wisler  (1993).  and  these  same  definitions 
were  adopted  by  the  Working  Group. 

For  obtaining  pitchwise  averages,  tire  measured  total 
pressures  and  temperatures,  axial  and  radial  velocities  are 


mass-averaged,  and  the  mean  angle  is  calculated  from  those 
averaged  velocities.  For  comparing  overall  performance  with 
measured  values,  the  preferred  method  is  to  interpolate  the 
pitchwise-mean  computed  values  at  the  radial  locations  of  tire 
measurement  points  given  in  Tables  2.1  and  2.2,  and  then 
mass-average  those  values  radially.  The  flow  area  associated 
with  each  measurement  location  is  also  given  in  those  tables. 
This  ensures  that  the  overall  performance  figures  are  directly 
comparable  despite  the  relatively  sparse  distribution  of  the 
measurement  points. 

The  equations  defining  these  averages  can  also  be  obtained  as 
described  in  the  next  Section. 

2.1.5  Test  Case  Data 

The  data  sets  selected  for  the  WG26  code  evaluation  included 
averaged  overall  performance  values,  rotor  entrance  and  exit 
(Stations  1  and  4)  probe  survey  information  and  selected 
laser  velocimeter  results.  These  data  as  well  as  the 
configuration  geometry  may  be  obtained  by  mailing  a  request 
to  Dr  K.  L.  Suder  at  the  NASA  Lewis  Research  Center, 
Cleveland,  OH  44135,  U.S.A.  Request  letters  should  include 
a  current  electronic  mail  address. 

2. 1 .5. 1  Defined  Values  for  CFD  Code  Input 
Turbulence  intensity  =  3%  at  station  1 . 

P/Pref  and  T/Trer  values  in  Table  2.1  for  Station  1  should  be 
used  in  describing  entrance  conditions  for  CFD  codes, 
including  calculation  of  m  l  mclloke 


Total  Pressure  Ratio 
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Table  2.1  Aerodynamic  Survey  Measurement  Locations  and  Flow  conditions  at  Station  1,  z  =  4.19  cm 


point 

%  span 

radius 

(cm) 

P/Pref 

T/Tref 

AAan  (cm2) 

%  area 

: 

1 

97.0 

25.4203 

0.9435 

1.0008 

59.3723 

0.0537278 

2 

94.0 

25.1765 

0.9762 

1.0004 

45.7621 

0.0414115 

3 

90.0 

24.8412 

0.9959 

1.0004 

57.0518 

0.0516279 

4 

85.0 

24.4450 

1.0020 

1.0002 

63.1837 

0.0571768 

5 

80.0 

24.0182 

1.0041 

1.0000 

62.1197 

0.0562140 

6 

75.0 

23.6220 

1.0048 

1.0002 

58.8042 

0.0532137 

7 

70.0 

23.2258 

1.0041 

1.0006 

60.0314 

0.0543242 

8 

65.0 

22.7990 

1.0034 

1.0008 

69.7699 

0.0631369 

9 

58.0 

22.2504 

1.0041 

1.0004 

78.8012 

0.0713096 

10 

51.0 

21.6713 

1.0048 

0.9996 

78.8530 

0.0713565 

11 

44.0 

21.0922 

1.0054 

0.9994 

74.7585 

0.0676512 

12 

37.0 

20.5435 

1.0054 

0.9990 

72.7605 

0.0658431 

13 

30.0 

19.9644 

1.0054 

0.9988 

61.3109 

0.0554821 

14 

25.0 

19.5682 

1.0054 

0.9988 

50.5746 

0.0457665 

15 

20.0 

19.1414 

1.0054 

0.9988 

49.5106 

0.0448036 

16 

15.0 

18.7452 

1.0054 

0.9987 

46.6640 

0.0422276 

17 

10.0 

18.3490 

1.0040 

0.9990 

47.4221 

0.0429137 

18 

5.0 

17.9222 

0.9864 

1.0004 

68.3072 

0.0618133 

hub  radius  =  17.5259  cm 
tip  radius  =  25.6692  cm 
area  =  1105.06  cm2 


Table  2.2  Aerodynamic  Survey  Measurement  Locations  at  Station  4,  z  =  10.64  cm. 


point 

%  span 

radius 

(cm) 

AAan  (cm2) 

%  area 

1 

97.0 

23.7134 

25.3625 

0.0420720 

2 

94.0 

23.5915 

22.5758 

0.0374494 

3 

90.0 

23.4086 

31.3590 

0.0520192 

4 

85.0 

23.1648 

33.2834 

0.0552114 

5 

80.0 

,22.9514 

30.7671 

0.0510373 

6 

75.0 

22.7381 

32.6485 

0.0541582 

7 

70.0 

22.4942 

32.3201 

0.0536135 

8 

65.0 

22.2809 

36.2286 

0.0600970 

9 

58.0 

21.9761 

42.0868 

0.0698147 

10 

51.0 

21.6713 

43.5641 

0.0722654 

11 

44.0 

21.3360 

42.9205 

0.0711977 

12 

37.0 

21.0312 

40.2774 

0.0668134 

13 

30.0 

20.7264 

35.7477 

0.0592993 

14 

25.0 

20.4826 

29.4308 

0.0488206 

15 

20.0 

20.2692 

27.1778 

0.0450834 

16 

15.0 

20.0558 

28.7959 

0.0477674 

17 

10.0 

19.8120 

28.4675 

0.0472227 

18 

5.0 

19.5986 

39.8218 

0.0660575 

hub  radius  =  19.381  cm 
tip  radius  =  23.823  cm 
area  =  602.835  cm2 


16 


2.2  DLR  CASCADE 


2.2.4  Uncertainty  and  evaluation  methods 


2.2.1  Test  Case 

The  annular  test  cascade  is  a  scaled  version  of  a  subsonic, 
low  aspect  ratio  turbine  stator.  This  stator  is  scaled  by  a 
factor  of  2.348  and  has  constant  hub  and  tip  radii.  The  aim  of 
die  tests  was  to  provide  detailed  information  about  die  flow 
field  inside  die  passage  as  well  as  upstream  and  downstream 
of  die  blade  row.  The  measurements  carried  out  were 
performed  widi  a  5-hole  probe  downstream  of  the  stator  and 
with  die  diree  dimensional  Laser-Two-Focus  (3D-L2F) 
measuring  technique  inside  the  passage  and  downstream. 
Additionally  the  upstream  boundary  layers  were  traversed 
using  a  “fish-mouth”  pitot  probe. 

2.2.2  Facility  and  Instrumentation 

Hie  tests  were  carried  out  in  die  DLR  annular  turbine  test  rig 
with  25  blades  mounted  on  the  stator  hub  (Fig.  2.6).  It  is  an 
open  loop  continuously  operating  facility  supplied  with  cold 
air  (maximum  temperature  315K,  maximum  pressure  2.0  bar) 
and  the  following  main  dimensions: 

hub  diameter,  dH  0.315  m 

tip  diameter,  dr  0.400  m 

The  operating  point  was  determined  by  measuring  the  mass 
flow,  the  total  pressure  and  the  total  temperature  upstream  of 
die  stator. 

The  5-hole  probe  used  to  measure  the  flow  field  downstream 
of  the  blade  row  is  mounted  in  computer  controlled  actuators 
which  allow  spanwise  traversing.  Pitchwise  traversing  was 
achieved  by  turning  the  complete  stator  hub  with  the  blades. 
The  resultant  tip  leakage  is  prevented  by  a  tiny  plastic  tip 
seal.  The  flow  field  was  measured  in  a  plane  40%  axial 
chord  length  behind  the  stator  measuring  plane  3  (MP3). 
There  die  number  of  radial  measuring  positions  was  18  and 
the  number  of  circumferential  positions  was  15,  arranged 
equidistandy. 

The  applied  L2F  measuring  technique  is  described  in  detail 
in  Schodl  (1989).  Because  die  3D-L2F  measurements  are 
very  time  consuming,  the  flow  field  is  highly  resolved  only  in 
the  regions  of  interest.  The  undisturbed  mainflow  is 
described  at  some  characteristic  points.  For  MP1  and  MP2 
the  nearest  measuring  points  to  the  suction  side  are  at  a 
distance  of  0.5  mm  from  the  surface.  The  distribution  of  the 
measuring  points  is  shown  by  the  beginning  of  the  contour 
plot  presentation  for  each  measuring  plane. 

The  axial  position  of  die  measuring  planes  (MP)  of  the  5-hole 
probe  measurements  and  the  Laser-Two-Focus  measurements 
are  shown  in  Figs.2.7  and  2.8. 


2.2.3  Cascade  and  geometry 

The  blade  geometry  at  mid-span  is  shown  in  Fig.2.7.  The 
main  dimensions  are  indicated  below: 


number  of  blades  25 

chord  length  at  hub,  Ch  0.0622  m 

Chord  lengdi  at  mid  span,  Cms  0.0698  m 

chord  length  at  tip,  Ct  0.0768  m 

axial  chord  lengdi  at  mid  span,  Cat,  ms  0.0445m 

aspect  ratio,  1i/Cms  0.61 


outlet  flow  angle  rel.  to  tangl  direction,  a3  20.5° 

The  3D-blade  coordinates  are  available  on  a  floppy  disk, 
which  may  be  obtained  from  C.  Langowsky,  DLR,  Institut  Air 
Antriebstechnik,  D-51 140,  KOln,  Germany. 


2.2.4. 1  5-hole  probe  measurements 

The  measurable  variables  of  the  5-hole  probe  are  evaluated 
diree  dimensionally.  For  this  application  die  probe  had  to  be 
calibrated  for  the  two  flow  directions  a,  p  and  the  Mach 
number.  Polynomial  approximations  were  used  to  relate  die 
calibration  data  to  die  measured  data  (Gallus  and  Bolui, 
1976).  A  correction  mediod  for  die  influence  of  die  head 
geometry  to  the  flow  witii  gradients  on  die  measured  results 
is  not  used.  Thus,  die  original  data  are  presented.  Unsteady 
effects  as  described  for  die  3D-L2F  measurements  could  not 
be  detected  with  the  5-hole  probe. 


The  estimated  uncertainty  of  die  5-hole  probe  data  is  as 
indicated  below: 

flow  angles  a,  p  ±0.5° 

pressures  pt,  pst  ±0. 1  mbar 


For  the  radial-distribution  plots,  two  circumferential 
averaging  techniques  were  used  depending  on  the  specific 
quantity  being  averaged.  These  are  die  area-average  and  the 
mass-average  defined  as  follows: 

Area-average  of  a  quantity  F: 


Mass-average  of  a  quantity  F: 


The  total  pressure  Pt  is  mass-averaged,  while  die  static 
pressure  Pst  is  area-averaged.  The  Mach  number  emerged 
from  the  quantities  of  die  averaged  total  pressure  and  static 

pressure,  M  =  ffpt  ,pst ).  Reported  average  flow  angles  are 
calculated  from  the  area-averaged  velocity  components. 


2.2.4.2  3D  Laser-Two-Focus  (3D-L2F)  measurements 
The  measured  variables  with  the  3D-L2F  are  the  projected 
mean  velocity  v  and  the  flow  angles  a  and  p.  For  die  results 
presented  here,  the  mean  velocity  v  is  converted  to  the  Mach 
number. 

Unsteady  effects  in  the  wake  region  (e.g.  von  Kamian  vortex 
streets)  could  be  resolved  witii  this  measuring  technique. 
The  statistical  evaluation  mediod  describes  an  averaged 
condition  different  from  die  true  alternating  flow  condition. 
These  unsteady  effects  have  an  influence  only  on  the 
measured  flow  angle  a,  which  differ  slightiy  in  the  wake 
region  at  measuring  plane  3  (MP3)  for  die  two  measuring 
techniques,  because  they  are  averaged  in  different  ways. 

For  die  measurements  a  window  of  quartz  glass  with  a  plane 
surface  is  necessary  in  the  casing.  An  unavoidable  tiny  tip 
leakage  results  due  to  the  deviation  to  the  contour  of  die 
casing.  For  that  reason  die  flow  field  is  slightly  disturbed  at 
the  tip.  The  laser  device  had  to  be  inclined  to  the  quartz 
glass  window  for  the  measurements  inside  the  passage, 
increasing  die  possible  measurable  region.  The  accuracy, 
dierefore,  arises  as  following: _ 


. t  ‘ 

s 

mean  velocity  v 

±  1 .5  m/s 

dk  1.5  m/s 

flow  angle  a 

±0.1° 

±0.2° 

flow  angle  p 

±0.4° 

±  1.5° 
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MEASURING  INSTRUMENTATION 


Fig  2.6  DLR  test  facility 


The  flow  variables  M,  a  and  p  are  based  on  the  area- 
averaged  velocity  components  for  the  radial  distribution 
plots. 

2.2.5  Experimental  conditions  and  results 

The  characteristic  data  of  the  operating  point  for  which  the 
measurements  were  carried  out  are  listed  below: 


Mass  flow  rate,  m0  5.490  kg/s 

total  pressure,  pto,  ms  1 .6760  bar 

total  temperature,  Tto  306. 6K 

inlet  flow  angle,  cto  (circumf.)  90° 

inlet  flow  angle,  po  (radial)  0° 

inlet  turbulence  level,  Tuo  4.4% 

inlet  Mach  number,  Mao  _  0. 1 76 

mean  outlet  Mach  number,  Ma3  0.74 
static  pressure  at  hub,  Pst,  3  1 .0750  bar 

Reynolds  number,  Re  1  x  106 


The  results  of  detailed  5-hole  probe  measurements  prove  that 
the  inlet  flow  angle  at  MP=0  is  uniform.  The  incoming 
boundary  layers  at  hub  and  tip  have  a  thickness  of  about 
4.5%  blade  height  and  their  profile  is  shown  in  Fig.2.9. 

Derived  from  the  5-hole  probe  measurements,  Figs.2. 10-1 3 
present  the  contour  plots  of  the  total  pressure  ratio,  the  Mach 
number,  the  circumferential  flow  angle  a  and  the  radial  flow 
angle  P  at  MP  3.  Tire  radial  distribution  of  the 
circumferential  averaged  total  pressure  ratio  po/pt,o,Ms ,  Mach 
number  M3,  static  pressure  ratio  pst,3/pt,o,Ms,  and 
circumferential  flow  angle  a  is  shown  in  Figs.2. 14-1 7.  Tire 
data  for  tire  5-hole  probe  measurements  are  also  available  on 
tire  disk  mentioned  in  Section  2.2.3. 

Tire  distribution  of  the  measuring  points  and  tire  results  of  tire 
3D-L2F  measurements  at  MP3  are  shown  in  Figs.2. 18-21. 
Tire  radial  distribution  of  the  averaged  tliree  flow  variables 


M,  a  and  p  is  presented  in  Figs.2.22-24.  Contour  plots  of  tire 
measurements  inside  the  passage  at  MP2  and  MP1  and  tire 
distribution  of  the  measuring  points  follow  hr  Figs.2.25-32. 
The  results  of  the  3D-L2F  measurements  are  also  stored  on 
tire  disk. 

2.2.6  Nomenclature  for  Section  2.2 
C  [m]  chord  length 

d  [mj  diameter 

h  [nr]  blade  height 

m  [kg/s]  mass  flow  rate 

M  Mach  number 

p  [bar]  pressure 

r  [m]  radius 

v  C 

Re  Reynolds  number  =  — — ■ 

v3 

Tu  [%]  turbulence  level 

v  [nr/s]  velocity 

x,y,z  cartesian  coordinates 

a  [°]  circumferential  flow  angle,  rel  to  tangl.direction 

P  [°]  radial  flow  angle 

5  [nr]  boundary  layer  thickness 

4>  [°]  stator  circumferential  position 

v  [m2/s]  kinematic  viscosity. 


Subscripts  and  Superscripts 

Svnrbols 

0...3 

measuring  plane 

MP 

measuring  plane 

ax 

axial  direction 

PS 

pressure  side 

H 

hub 

SS 

suction  side 

MS 

nrid-span 

TE 

trailing  edge 

st 

static  condition 

t 

total  condition 

T 

tip 

averaged 

5.00  20.00  25.00 


Results  of  the  3D-L2F  measurements  at  MP3 
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Chapter  3 

THE  COMPUTATIONS 


3.1  INTRODUCTION 

The  members  of  the  Working  Group  invited  research  workers 
to  undertake  computations  of  tire  test  cases  using  any  3D 
RANS  codes  they  wished  to  use.  Computations  using  fifteen 
different  codes  were  supplied  to  the  WG,  in  some  cases  with 
a  range  of  grids  and  turbulence  models.  All  the  codes  solved 
tire  steady  flow  RANS  equations.  Most  of  tire  computations 
were  undertaken  specifically  for  the  WG.  Some  solutions 
obtained  earlier  (perhaps  as  a  result  of  ASME’s  use  of  Rotor 
37  as  a  test  case)  were  also  provided.  Not  all  the  codes  in 
common  use  in  industry  were  included,  but  all  the  important 
types  were  represented,  hi  all,  a  wide  range  of  grids  and 
turbulence  models  were  employed. 

Tables  3.1  and  3.2  give  details  of  the  codes  used.  These 
details  are  collected  in  this  way  for  reference  purposes,  but  in 
the  analysis  sections  of  this  chapter  specific  details 
considered  to  be  of  particular  importance  are  repeated. 

As  both  test  cases  involved  only  an  isolated  blade  row,  the 
question  of  how  to  “average”  or  otherwise  to  model  in  a 
steady  code  the  unsteady  aerodynamic  interference  between 
adjacent  blade  rows  did  not  arise  in  the  present  study.  The 
WG  recognises  this  as  a  key  issue  to  be  addressed  at  a  later 
date,  when  unsteady  codes  have  become  established. 

3.2  THE  RESULTS  FOR  ROTOR  37 

3.2.1  Computations  requested 

hi  1993,  The  ASME  Turbomachinery  Committee  invited 
research  workers  to  predict,  using  RANS  codes,  the  overall 
mass  flow  range,  pressure  ratio,  temperature  ratio,  and 
efficiency  of  Rotor  37  at  design  speed,  and  also  to  predict  at 
both  92.5%  and  98%  of  the  choked  flow  the  spanwise 
distribution  of  pitchwise-mean  pressure  ratio,  temperature 
ratio,  and  efficiency.  (Only  the  geometric  data  not  the  flow 
measurements  were  provided  at  that  time,  but  contributors  to 
the  present  study  also  had  the  measurements.)  These  same 
parameters  were  also  specified  to  the  WG;  but  in  addition 
contributors  were  also  asked  to  provide  certain  parameters  in 
the  form  of  contour  plots  to  enable  a  thorough  analysis  to  be 
undertaken. 

The  overall  performance  predictions  and  the  pitchwise-mean 
traverse  predictions  at  98%  of  the  choked  flow  are  shown  in 
Figs  3. 1-3.5  in  comparison  with  the  measurements. 
(Relatively  little  attention  was  given  to  the  stall  point 
results.)  Owing  to  the  enormous  volume  of  the  three- 
dimensional  results,  only  selected  plots  will  be  shown  to 
illustrate  specific  points  in  the  analysis. 

3.2.2  Flow  description 

This  description  of  the  flow  is  based  on  published  results. 
Denton  (1996)  has  given  a  good  global  analysis  of  the  flow  in 
this  compressor.  Chima  (1996b)  and  Suder  and  Celestina 
(1996)  analysed  the  tip  leakage  flow,  and  Hah-Loellbach 
(1997)  and  Shabbir  et  al  (1997)  analysed  the  hub  comer  flow. 

Three  regions  are  considered:  the  mid-span,  the  hub  wall  and 
the  tip  wall  regions. 

(1)  The  mid-span  region  is  dominated  by  a  strong  shock 
attached  at  the  blade  leading  edge.  This  shock  interacts 
strongly  with  tire  suction  side  boundary  layer.  The 


boundary  layer  after  the  shock  may  separate  either  up  to 
the  trailing  edge  according  to  some  authors,  or  reattach 
before  the  trailing  edge.  A  strong  radial  movement  is 
also  observed  in  the  separated  area  from  the  hub  to  the 
tip  wall,  as  for  example  in  Fig  3.6. 

(2)  Near  the  hub  wall,  the  measurements  of  absolute 
stagnation  pressure  P0  at  station  4  show  a  region  of  low 
pressure  at  20%  of  the  blade  span,  for  both  98%  and 
92%  of  the  choked  mass  flow.  Strazisar  has  pointed  out 
that  this  pressure  deficit  is  present  over  a  wide  range  of 
mass  flow  for  the  nominal  speed  of  rotation  (Fig.3.7). 
He  mentions  also  that  the  stagnation  temperature 
distribution  T0,  deduced  from  the  Euler  energy 
equation,  shows  the  same  trend  at  station  3 
immediately  downstream  of  the  rotor  trailing  edge.  This 
indicates  that  the  spanwise  character  of  the  P0 
distribution  is  not  an  artifact  of  the  radial  mixing. 

There  are  strong  indications  that  a  comer  stall  occurs 
near  the  suction  side  comer  (Fig.3.8).  This  comer  stall 
greatly  reduces  the  axial  momentum  in  the  region  of  the 
stagnation  pressure  deficit  near  the  hub  (Fig.3.9).  The 
axial  momentum  is  then  redistributed  all  over  the  blade 
span,  thereby  reducing  the  amount  of  the  blade  work, 
and  therefore  the  overall  pressure  ratio. 

This  comer  stall  seems  mainly  influenced  by  the 
supersonic  Mach  number  ahead  of  the  leading  edge,  in 
the  hub  region.  Both  experiments  and  simulations  do 
not  show  the  strong  P0  deficit  if  they  are  performed 
with  a  lower  subsonic  Mach  number  ahead  of  the  hub 
leading  edge  (Shabbir  et  al.,  1997,  Hah,  Loellbach, 
1997).  This  stall  may  be  partially  due  to  the  glancing 
side  wall  shock  wave  interaction  with  the  hub  boundary 
layer,  as  shown  in  Fig  14  ofPovinelli  (1997). 

Two  other  effects  have  a  very  important  influence: 

•  A  change  of  the  overall  mass  flow  at  the  nominal  speed 
strongly  modifies  the  radial  P0  distribution  over  the 
whole  blade  span.  As  the  mass  flow  is  reduced,  the  To 
deficit  is  always  present  at  the  hub,  while  the  Po  deficit 
is  slightly  smeared  out  for  the  highest  mass  flow  (Fig. 
3.7.) 

•  A  possible  leakage  flow  occurs  through  the  axial  gap 
between  the  non-rotating  upstream  hub  and  the  rotor 
(Fig.3.9).  Shabbir  et  al  (1997)  report  experiments  and 
CFD  studies  on  the  effect  of  modifying  the  upstream 
axial  gap  between  the  non-rotating  hub  and  the  rotor. 
All  these  studies  show  an  increase  of  the  stagnation 
pressure  downstream  near  the  hub,  with  a  reduction  of 
the  upstream  axial  gap.  Their  simulations  show  a 
deficit  of  stagnation  pressure  at  20%  from  the  hub  in 
agreement  with  the  experimental  observations.  More 
important,  the  simulations  shows  this  deficit,  either 
with  a  very  low  leakage  mass  flow  introduced  of  the 
order  of  0.25  to  0.33%  of  the  choking  mass  flow,  or 
even  when  there  is  no  net  mass  flow  through  a  gap 
connected  with  a  blind  cavity.  For  this  last 
configuration,  the  upstream  wave  of  the  leading  edge 
shock  could  drive  a  circumferential  inflow/outflow 
pattern  through  the  upstream  gap.  Most  of  the  CFD 
simulations  have  neglected  this  upstream  gap.  Their 


Table  3.1a  Code  information  for  Rotor  37  case 
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Fig  3  3  Rotor  37  total  pressure  ratio 
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(b)  turbulent  transport  models 


Fig  3  4  Rotor  37:  total  temperature  ratio 


(b)  turbulent  transport  models 
Fig  3.5  Rotor  37  efficiency 
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Fig  3.7  Sensitivity  of  performance  to  mass  flow  (1  is  the  highest  flow  and  7  the  lowest  flow) 


overprediction  of  P0  in  the  hub  region,  and  all  over  the 
blade  span,  is  also  consistent  with  hub  leakage  flow 
influences. 

(3)  The  tip  wall  region  is  dominated  by  a  strong  tip  leakage 
flow.  Chinra  (1996b)  and  Suder  and  Celestina  (1996) 
have  presented  some  detailed  CFD  results  for  the  flow 
in  this  region.  Detailed  experimental  results  and 
analysis  have  also  been  given  by  Suder  and  Celestina. 
The  main  conclusions  of  these  studies  are: 

(a)  The  leakage  flow  issuing  from  the  suction  side, 
upstream  of  the  passage  shock,  behaves  almost  as  an 


inviscid  jet  owing  to  a  strong  supersonic  expansion 
(Fig.  3.10). 

(b)  The  limiting  surface  between  the  tip  leakage  jet  and  the 
main  crossflow  is  a  region  of  strong  shear,  owing  to  the 
very  large  difference  of  direction  between  the  two  flows 
(Fig.3.10).  This  shear  layer  is  attached  to  the  blade 
comer  on  the  suction  side.  Vorticity  may  be  created  at  a 
particularly  high  rate  in  this  layer. 

c)  At  the  exit  from  tire  tip  gap,  the  particle  traces 
presented  by  Chima,  and  Suder  and  Celestina,  exhibit 
clearly  a  negative  axial  velocity  (Fig.3.11).  The  strong 
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Axial  Momentum 


Fig  3.8  Computed  particle  traces  at  99%  choke 
flow  (Hah  &  Loellbach,  1997) 


Fig  3.10  Relative  Mach  no.  contours  on  a  cross- 
section  at  10%  chord,  near  stall  (Chima,  1996b) 

pressure  gradient  across  the  tip  gap  between  the  suction 
and  pressure  sides  dominates  the  gradient  along  the 
blade  chord.  As  a  consequence,  the  leakage  flow 
velocity  vector  points  in  a  direction  roughly  normal  to 
the  blade  suction  surface. 

(d)  The  leakage  flow  and  the  main  crossflow,  both 
supersonic  flows,  merge  along  a  line  connected  with  the 
leading  edge.  The  gradients  of  axial  velocity  component 
are  very  high  across  this  line.  Chima  has  shown  that 
Chen's  model  (1991)  predicts  remarkably  well  the 
global  flow  path  along  this  line.  This  is  thought  to  be  a 
three-dimensional  separation  line,  along  which  the 


Fig  3.9  Effect  of  leakage  flow  on  momentum  mass 
distribution  (  Shabbir  et  al,  1997) 


Fig. 3.1 1  Shock  system  above  95%  span  at  peak 
efficiency,  with  tip  particle  traces  (Chima,  1996b) 

axial  velocity  is  very  small  (Perrin  et  al,  1992).  As  soon 
as  the  leakage  flow  is  ejected  from  the  casing  wall,  it 
wraps  itself  up  in  a  vortical  movement.  The  origin  of 
this  vorticity  is  certainly  linked  to  the  shear  layer 
attached  to  the  suction  side  comer  (see  (b)  above). 

(e)  The  interaction  of  this  3D  separation  line  and  the 
downstream  leg  of  the  passage  shock  occurs  at  mid¬ 
distance  between  the  pressure  and  suction  sides.  A 
strong  area  of  low  velocity  then  appears  downstream  of 
the  shock.  This  area  extends  over  most  of  the  blade 
passage  as  the  mass  flow  decreases.  The  generation  of  a 
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Fig  3.12  Mach  number  contours  at  peak 
efficiency  (Chima,  1996b) 


low  longitudinal  velocity  area  downstream  of  the 
passage  shock  is  generated  by  an  inviscid  phenomenon. 
As  mentioned  by  Schlechtriem  and  LOtzericlr  (1997), 
the  interaction  of  a  longitudinal  clearance  vortex  with 
the  passage  shock  creates  a  transverse  vorticity  in  the 
direction  parallel  to  the  tip  wall,  that  could  lead  to  a 
flow  separation.  Chima  (1996b)  shows  a  separation  in 
this  area  on  the  tip  wall,  which  extends  downstream  of 
the  passage  shock  over  30-40%  of  the  blade  chord  (Fig. 
3.12).  This  part  of  the  leakage  flow  reaches  the  trailing 
edge,  near  the  pressure  side  of  blade. 

(f)  The  leakage  flow,  introduced  downstream  of  the 
passage  shock,  deviates  the  separated  boundary  layer  on 
tire  suction  side  that  migrates  radially  towards  the  tip 
wall.  This  forms  a  second  vortical  movement 
particularly  detected  at  part  speed  (Suder  and 
Celestina,  1996).  This  second  vortex  does  not  merge 
with  the  first  leakage  vortex  according  to  the  simulation 
of  Suder  and  Celestina  (Fig.3.13). 

(g)  Tire  mixing  process  of  the  leakage  flow  with  the 
primary  flow  extends  over  10%  of  the  blade  span, 
which  is  a  distance  equal  to  20  times  the  tip  clearance 
gap.  This  mixing  process  extends  also  probably  far 
downstream. 

The  difficulties  of  analysing  tire  simulated  results  arise  from 
the  lack  of  some  plotted  quantities  in  the  experimental  data, 
and  in  tire  simulated  results  as  well.  For  example,  the  axial 
velocity  component  has  not  been  plotted,  although  it  was 
either  measured  or  computed.  As  a  consequence,  tire  velocity 
triangles  are  difficult  to  analyse,  and  it  is  not  possible  to 
deduce  always  a  firm  conclusion  from  tire  comparisons 
between  tire  experimental  and  CFD  results. 

Two  main  sets  of  data  are  considered: 

•  Blade-to-blade  contour  plots  of  Mrei  are  given  at 
different  percentages  of  tire  blade  height;  pitchwise 
plots  of  Mrei  are  given  after  tire  trailing  edge,  at  stations 
3  and  4,  and  at  20%  chord  inside  tire  rotor  at  station  2. 


•  Radial  plots  of  cicunrferentially  averaged  quantities  are 
given  at  stations  3  and  4.  Most  of  the  authors  have 
plotted  the  absolute  stagnation  pressure  Po/TW,  the 
absolute  stagnation  temperature  To/T(er,  the  absolute 
flow  angle  a  =  tan'fVq/Vz),  and  tire  adiabatic  efficiency 
His- 

Two  supplementary  sets  of  plots  are  sometimes  added: 

•  Contour  plots  are  given  in  different  surfaces  (tire 
meridional  surface,  a  section  nomral  to  the  axis  of 
rotation,  or  on  tire  blade  surface)  for  tire  above 
quantities,  plus  the  turbulent  viscosity  pt/p,  the  entropy, 
and  the  static  pressure. 

•  Particle  traces  show  details  of  tire  tip  leakage  flow,  tire 
hub  flow  or  the  suction  side  flow  behaviour. 

3.2.3  Overall  performance 

The  overall  adiabatic  efficiency  is  generally  predicted  to  be  at 
a  lower  level  than  in  the  experiments  (Fig.3.2).  The  radial 
plots  of  the  same  quantities  at  station  4  (Fig  3.5)  show  an 
underprediction  of  3%  on  efficiency  between  10  and  80%  of 
the  blade  span.  Simultaneously,  there  is  a  strong  deficit  of 
efficiency  for  most  of  the  simulations  from  80%  to  the  tip 
wall.  It  is  concluded  that  the  tip  wall  region  is  an  important 
origin  of  the  global  loss  overprediction  in  tire  simulations. 
Hah's  simulation  produces  a  higher  value  of  the  global 
efficiency,  this  seems  to  be  associated  with  an  underestimate 
of  tire  losses  from  30%  to  85%  of  the  span.  It  produces  also 
an  overestimate  of  the  losses  near  the  side  walls  for  98%  of 
the  choked  mass  flow. 

Most  of  the  simulations  generate  a  higher  pressure  ratio  than 
tire  measured  one  (Fig.  3.1).  A  minority  produces  a  lower 
value  (Kang  and  Hirsch,  CANARI  using  Michel’s  turbulence 
model).  Two  sets  of  results  are  close  to  the  experiments  at 
98%  of  choked  mass  flow  (TRACE-S,  CANARI  using  the  k-e 
turbulence  model),  and  two  others  (CMOTT-CKE, 
TASCflow)  are  close  to  experiments  at  the  near  stall 
condition.  The  highest  value  of  pressure  ratio  is  obtained  by 
Hah;  that  is  consistent  with  the  low  losses  at  mid  span  as 
noted  above. 

The  following  remarks  refer  to  the  CFD  results  compared 
with  tire  experimental  ones  in  the  form  (Qamuiation  -  Qexpoiment). 
Tire  pitchwise  plots  for  98%  of  the  reduced  mass  flow  show 
that,  at  station  2  at  20%  chord  inside  the  rotor,  (Chima, 
Denton,  1996;  Fig.3.14): 

1)  the  computed  Mmi  is  always  lower  downstream  of  the 
shock, 

2)  the  shock  is  located  too  far  upstream,  provided  a  fine 
mesh  is  used. 

Part  of  the  difference  may  be  due  to  shock  smearing  in  tire 
calculation,  though  the  smearing  seems  more  important  in  the 
experimental  results  (Fig  3.14).  A  small  particle  lag  effect 
could  also  have  displaced  downstream  the  experimental 
location  of  the  shock  wave. 

At  station  3,  (Chima,  Denton,  1996,  McNulty,  Fig.3.15): 

3)  tire  wake  is  always  too  deep, 

4)  the  wake  position  at  station  3  is  for  most  of  the 
simulations  in  agreement  with  the  experiments. 

At  station  4, 

5)  tire  wake  deficit  is  strongly  reduced,  although  tire 
experimental  deficit  is  overpredicted  (Chima),  provided 
a  sufficiently  fare  grid  is  used  downstream. 


Second 


Blade-to-blade  view  of  particle  traces  indicating  Perspective  view  of  particle  traces  indicating 

the  formation  of  the  leakage  and  second  vortex.  the  formation  of  the  leakage  and  second  vortex. 


Fig  3.13  Computed  particle  traces  at  60%  speed  (Suder  and  Celestina,  1996) 
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Fig  3.14  Mach  number  variation  across  the  pitch  at  mid-span  (Hildebrandt) 


Relative  Mach  Number  Ma 
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Fig  3.15  Mach  number  variation  across  the  pitch 
at  mid-span  (Hildebrandt) 

Some  hypotheses  are  proposed  for  the  origins  of  the 
differences  between  the  results  of  the  simulations  and  the 
experiment.1 

a)  For  a  fixed  passage  shock  strength,  the  low  post-shock 
Mach  number  Mrei  (observation  no.l)  could  be  induced 
by  too  thin  a  post-shock  boundary  layer  on  the  suction 
side  (Denton,  1996).  This  is  because  the  subsonic  Mach 
number  downstream  of  the  shock  is  reduced  if  the  flow 
blockage  diminishes.  This  implies  a  lower  predicted 
loss  in  the  wake.  However,  this  hypothesis  is  in 
contradiction  with  (3)  the  deeper  simulated  wakes  at 
station  3  (Denton,  1996). 

b)  It  was  also  observed  that  (2)  the  shock  is  slightly 
displaced  upstream,  implying  a  stronger  shock  than 
expected.  This  stronger  shock  is  compatible  with  the 
lower  Mach  number  downstream  of  the  shock.  This  is 
obvious  from  tire  blade-to-blade  plots  of  Chima  or 
Hildebrandt  for  instance.  As  the  shock  is  displaced  on 
the  suction  side  towards  the  leading  edge,  the 
downstream  static  pressure  is  higher  than  expected, 
implying  that  the  mass  flow  used  in  the  simulation  is 
too  low.  Chima's  simulations  have  shown,  by  means  of 
a  downstream  adjustment  of  the  static  pressure,  that  an 
increase  of  0.24%  for  the  mass  flow  (from  98%  to 
98.24%)  is  sufficient  to  fit  the  experimental  blade-to- 
blade  Mach  number  distribution  (Fig.3.16).  Note  that 
while  this  modification  has  a  strong  influence  on  the 
shock  location  and  tire  downstream  Mach  number 
distribution,  it  has  almost  no  effect  on  the  upstream 
Mach  number  value.  This  value  (0.24%)  is  also  of  the 
same  order  as  the  difference  between  the  averaged 
choked  mass  flow  deduced  from  the  simulations  (20.86 
kg/s)  and  the  experimental  value  (20.93  kg/s).  It  may  be 
thought  necessary  to  compensate  for  the  small  error  in 
die  predicted  choked  mass  flow  by  making  a 
comparison  at  a  slightly  higher  mass  flow  than  the 
required  98%  value.  However,  this  small  difference  on 
tire  choked  mass  flow  will  not  explain  alone  the  strong 
discrepancies  between  tire  experimental  and  simulated 


1  The  numbers  in  parentheses  designate  tire  observations 
quoted  at  the  beginning  of  tire  paragraph 


Fig  3.16  Comparison  of  three  solutions  nominally  at 
98%  choked  flow  (Chima) 

results  observed  at  20%  of  tire  blade  height  from  tire 
hub. 

c)  The  observation  (3)  of  deeper  simulated  wakes  in 
station  3  is  a  common  feature  of  all  the  simulations.  It 
may  be  thought  that  some  flow  phenomena  are  not 
simulated.  The  rotor  wake  deficit  (M^  max  -  Mrei,  min) 
/Mrd,  max  has  been  plotted  as  a  function  of  the  axial 
distance  by  McNulty  for  70%  of  tire  blade  span  (Fig. 
3.17).  It  shows  that  the  experimental  deficit  is  1/3 
lower  than  in  the  simulation  near  the  trailing  edge.  But 
the  wake  may  be  unsteady  due  to  periodic  vortex 
shedding.  If  so,  a  laser  velocimeter  will  detect  two 
velocity  minima  at  slightly  different  positions,  and  the 
“average”  wake  will  be  less  deep  than  either  of  the  two 
time-resolved  minima.  So  it  is  not  clear  whether  the 
difference  between  tire  simulations  and  the 
measurements  is  due  to  unsteadiness  in  the  real  flow,  to 
forcing  tire  simulations  to  a  “steady”  solution,  or  to  the 
turbulence  model. 

d)  It  must  be  remembered  that  tire  experimental  flow 
gradients  are  strongly  dependent  on  the  mass  flow.  The 
gradients  computed  at  station  4  are  also  a  consequence 
of  tire  amount  of  loss  generated  in  tire  rotor  and 
downstream  of  tire  rotor.  The  slope  of  the  wake  decay  is 
higher  in  the  simulation  (Chima)  after  a  distance  of 
150%  chord  downstream.  The  higher  wake  dissipation 
observed  in  tire  simulation  between  stations  3  and  4  (5) 
may  be  comiected  with  the  mesh  topology  used  in  the 
downstream  zone  for  most  of  the  simulations.  This 
suggests  the  need  to  improve  the  mesh  density  and  the 
grid  alignment  with  tire  wake  direction  downstream  of 
tire  rotor. 

3.2.4  Radial  plots  for  98%  of  choked  mass  flow 

3.2.4. 1  Stagnation  temperature  plots 

•  At  mid-span,  Kang  and  Hirsch's  results  are  slightly  too 
low,  while  CMOTT,  SWIFT,  CANARI  using  both 
Michel’s  turbulence  model  and  tire  k-e  turbulence 
model,  and  TRACE-S  agree  with  tire  experimental 
values,  and  all  tire  other  simulations  produce  too  high  a 
level  of  temperature. 
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NEAR  PEAK  EFFICIENCY  (98%  Choke)  -  70%  SPAN 


Percent  Axial  Chord 


Fig  3.17  Rotor  wake  decay  (McNulty) 


Fig  3.18  Variation  of  performance  with  mass  flow  (Hah  and  Loellbach,  1997) 


•  Near  the  casing,  all  simulations  produce  a  strong 
increase  of  T0  after  90%  of  the  span,  except  TASCflow1 
which  shows  a  strong  decrease. 

•  Near  the  hub  wall,  in  between  0-8%  of  the  span,  all  the 
simulations  produce  an  increase  of  temperature,  except 
for  TASCflow  that  shows  a  strong  decrease,  and 
TRANSCode-Baldwin-Lomax  which  produces  a 
constant  value. 

3. 2.4. 2  Stagnation  pressure  plots 

It  has  already  been  mentioned  that  ojily  a  few  sets  of 
results  produced  an  overall  pressure  ratio  lower  than 


1  As  the  report  was  going  to  press,  Dr  Hutchinson  informed 
the  WG  that  a  revised  treatment  of  the  energy  equation  in  the 
near-wall  region  had  recently  been  incorporated  in 
TASCflow,  now  called  CFX-TASCflow,  winch  had  greatly 
improved  tire  predictions  of  overall  performance  for  Rotor  37. 


the  experimental  values.  At  mid-span,  Kang  and  Hirsch 
and  CANARI  using  Michel’s  turbulence  model  produce 
a  low  P0,  two  sets  (TRACE-S,  CANARI  using  the  k-e 
turbulence  model)  are  very  close  to  the  experiments, 
while  all  the  others  produce  a  higher  P0. 

•  Near  the  casing,  the  slope  of  the  decrease  towards  the 
wall  is  greater  than  the  experiments  for  CANARI  using 
Michel’s  turbulence  model,  and  the  Kang  and  Hirsch 
and  SWIFT  predictions. 

•  Near  the  hub  wall,  only  Hall's  simulations  give  the 
decrease  of  P0  observed  under  30%  of  the  blade  height 
for  the  98%  of  the  choked  mass  flow.  However  some  of 
them  show  a  sort  of  plateau  in  this  region  (Kang  and 
Hirsch,  CANARI  using  the  k-e  turbulence  model).  Most 
give  a  supplementary  increase  in  between  0-8%  of  the 
blade  height.  Hah  and  Loellbach  (1997)  report  also  a 
much  better  agreement  with  the  experimental  results 
for  P0  and  T0  if  the  mass  flow  is  increased  from  98%  to 
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NASA  Rotor  37  Cavity  Analysis 


Predicted  Rotor  Exit  Profiles  (No  Cavity,  Fine  Mesh) 
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Fig  3.19  Influence  of  overall  pressure  ratio  on  rotor  exit  flow:  coarse  mesh  and  fine  mesh,  no  cavities  (Hall) 
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99%  (Fig.3.18).  Similar  results  are  reported  by 
Delaney.  An  increase  of  the  mass  flow  by  only  0.3% 
generates  the  decrease  of  P0  at  30%  of  the  blade  height 
(Fig.3.19).  A  sufficient  mesh  density  in  the  radial 
direction  seems  very  important  in  order  to  capture  this 
P0  deficit.  Note  however  that  the  P0  hub  deficit  exists  in 
the  experiment  for  all  the  mass  flows  (Strazisar). 

•  The  existence  of  a  cavity  between  the  rotating  rotor  hub 
and  the  upstream  fixed  hub  strongly  influences  the  P0 
deficit  (Fig. 3. 20),  (Delaney,  Shabbir  et  al,  1997). 

3.2.4.3  Adiabatic  efficiency  plots 

•  Most  of  tlie  predictions  give  a  good  distribution  at  mid¬ 
span,  except  SWIFT  which  is  too  high,  and  Kang  and 
Hirsch  and  TASCflow  which  are  slightly  too  low. 

•  Near  the  casing,  most  of  the  simulations  produce  a 
lower  efficiency  than  expected,  except  CMOTT  and 
TASCflow.  Near  the  hub  wall,  some  of  the  simulations 
give  an  increase  of  the  efficiency  towards  the  wall 
(CMOTT,  Shabbir,  CAN  API  using  Michel’s  turbulence 
model,  TRACE-S). 

3.2.5  Hypotheses  for  the  origins  of  the  differences 

between  results  of  simulations  and  experiment 
Unfortunately,  most  of  the  radial  plots  are  provided  at  station 
4  only.  Station  4  is  also  located  near  or  even  at  the 
downstream  boundary  of  the  mesh.  The  numerical  dissipation 
is  believed  to  be  too  high  between  the  rotor  trailing  edge  and 
station  4,  owing  to  the  size  of  the  cells.  As  a  consequence,  it 
is  very  difficult  to  derive  firm  conclusions. 

a)  The  radial  stagnation  temperature  distribution  gives  a 
good  estimate  of  tire  level  of  tire  axial  velocity 
component  Vz.  To  explain  this  point,  assume  that  the 
relative  flow  angle  p  is  correctly  computed.  This  is  a 
reasonable  hypothesis,  considering  for  instance  tire 
predicted  location  of  tire  wake  in  station  3.  Some 
authors  provide  P  radial  distributions  (ADPAC- 
McNulty)  with  a  good  agreement  between  the  results  of 
simulation  and  experiment.  An  overprediction  of  tire 
work  is  then  generated  by  a  low  value  of  Vz,  and 
conversely.  A  low/high  mass  flow  is  also  consistent 
with  a  high/low  pressure. 

b)  Near  the  tip  wall,  the  overestimate  of  tire  work  may  be 
linked  with  either  a  low  axial  velocity,  associated  with 
high  losses,  and/or  a  decrease  of  the  relative  flow  angle 
p.  Note,  however,  that  tire  flow  Mach  number  is  very 
sensitive  to  small  variations  of  tire  blade-to-blade 
relative  flow  direction  as  the  exit  Mach  number  is 
nearly  sonic.  According  to  tire  efficiency  distribution, 
high  losses  occur  near  tire  tip  wall  for  all  the 
simulations,  except  for  CMOTT  and  TASCflow.  For 
CMOTT,  there  should  be  a  low  absolute  value  of  the 
relative  flow  angle  P  near  tire  wall.  For  TASCflow, 
there  is  actually  air  underprediction  of  tire  work  at  tire 
wall,  which  is  then  compatible  with  low  losses  and  a 
high  velocity.  By  comparison  with  tire  measurements, 
tire  flow  blockage  generated  by  the  tip  wall  layer  should 
then  create  air  increased  axial  velocity  in  tire  mid-span 
region,  with  a  corresponding  reduction  of  the  work. 
This  is  not  always  observed  in  practice,  (except  for 
Kang  and  Hirsclr's  results)  perhaps  because  of  tire 
influence  of  the  error  in  tire  choked  mass  flow,  as 
shown  by  Clrinra. 

c)  For  most  simulations,  near  tire  hub  wall,  a  high  value  of 
tire  work  is  observed  with  a  high  efficiency,  and  a  high 
stagnation  pressure.  This  is  compatible  with  a  small 
thickness  of  tire  wnll  shear  layer,  and  then  a  high  value 


of  the  axial  velocity  and  a  high  flow  deflection  in  tire 
relative  frame  of  reference  (passage  vortex  effect).  This 
is  particularly  true  for  CMOTT,  Shabbir,  CANARI 
using  Michel’s  turbulence  model,  and  TRACE-S.  The 
flow  in  tire  hub  region  is  probably  dominated  in  tire 
experiment  by  a  strong  comer  stall  near  the  suction  side 
(Hair  and  Loellbach,  1997,  Delaney,  Kang  and  Hirsch). 
This  comer  stall  is  extremely  sensitive  to  small 
modifications  of  tire  boundary  conditions,  such  as  tire 
global  mass  flow  variations  that  act  upon  tire  Mach 
number  ahead  of  tire  leading  edge  at  tire  hub,  tire 
upstream  cavity,  or  tire  radial  mesh  distribution. 

It  is  important  to  notice  that  tire  flow  at  nrid  span  is  also 
strongly  dependent  upon  tire  flow  near  tire  hub  and  tip  walls. 
This  is  typical  for  this  sort  of  supersonic  compressor.  Tire 
large  loss  that  should  occur  in  the  hub  region  acts  as  a  strong 
aerodynamic  blockage.  It  pushes  the  mass  flow  towards  mid- 
span,  thereby  reducing  the  blade  work  at  mid  span  (Kang  and 
Hirsch).  An  exception  is  observed  for  Hah  and  Loellbaclr's 
results:  although  they  obtain  a  large  comer  stall  at  the  hub, 
they  overpredict  the  pressure  ratio  at  mid  span.  Owing  to 
their  parallel  overestimation  of  tire  efficency  at  mid  span,  it 
may  be  suspected  that  the  shock  boundary  layer  loss 
generation  is  minimised  in  their  simulation. 

3.2.6  Plots  of  turbulent  quantities 

Detailed  plots  of  turbulent  viscosity  have  been  provided  by 
Calvert  with  TRANSCode,  McNulty  and  Heidegger  with 
ADPAC,  Hildebrandt  with  TRACE-S,  Couaillier  with 
CANARI,  and  Clrima  with  SWIFT.  Although  very  different 
turbulent  models  have  been  used,  all  tire  results  present 
similar  trends. 

At  35%  of  the  blade  height  near  the  trailing  edge  at  station  3, 
tire  flow  behaves  as  a  quasi-3D  flow  outside  the  wake  or 
viscous  layers.  Tire  Mach  number  is  also  ahnost  constant  over 
a  large  part  of  the  blade  span,  outside  tire  wake  region.  The 
turbulent  viscosity  has  a  significant  value  only  in  tire  wake 
area,  where  a  strong  radial  movement  is  also  observed.  (Fig. 
3.21). 

By  contrast,  at  90%  blade  height  (Fig.3.22),  tire  flow  feels 
the  influence  of  tire  tip  wall.  Tire  increase  of  the  turbulent 
viscosity  on  tire  pressure  side  near  the  trailing  edge  may  be 
tire  result  of  tire  accumulation  of  tire  tip  leakage  flow  that 
crosses  tire  blade  passage  from  the  suction  side  (Clrinra, 
1996b).  Simultaneously,  tire  contour  maps  of  Mrd  in  tire 
blade-to-blade  section  show  a  strong  reduction  of  tire  Mach 
number  hr  the  middle  of  tire  blade  passage. 

This  interaction  of  the  primary  flow  with  tire  leakage  flow 
may  be  then  dominated  by  a  separation  of  tire  wall  layer,  as 
shown  by  Clrinra,  that  could  create  the  very  high  values  of 
turbulent  viscosity,  particularly  for  tire  low  mass  flow. 

3.2.7  Assessment 

1 )  Most  of  tire  simulations  generated  a  low  axial  velocity 
at  mid-span.  This  could  be  induced  by  tlrree 
phenomena.  Firstly,  most  of  tire  simulations  predict  a 
slightly  low  value  of  tire  choked  mass  flow.  It  can  then 
be  an  error  to  compare  simulations  at  a  fixed 
percentage  of  tire  simulated  choked  mass  flow. 
Secondly,  most  of  lire  simulations  do  not  capture  tire 
hub  comer  stall;  this  is  probably  tire  most  important 
effect.  As  a  consequence,  a  higher  proportion  of  tire 
mass  flow  passes  through  the  hub  region,  thereby 
decreasing  the  axial  velocity  at  mid-span.  Thirdly,  tire 
shock/boundary  layer  interaction  may  not  be  predicted 
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Fig  3  23  Viscosity  ratio  at  station  3  (Calvert) 


correctly.  A  more  detailed  study  of  the  results  is 
suggested,  to  establish  which  explanation  is  correct  in 
this  case. 

2)  The  tip  region  generates  too  much  loss  in  most  of  the 
simulations. 

3)  In  the  hub  region,  the  predicted  pressure  ratio  and 
temperature  ratio  are  both  too  high,  because  the  comer 
stall  is  either  not  predicted  or  is  minimized.  The  comer 
stall  has  been  linked  by  Shabbir  et  al  (1997)  to  the 
existence  in  the  experiment  of  an  axial  gap  between  the 
upstream  fixed  hub  and  the  rotor.  The  predicted  comer 
stall  is  also  very  sensitive  to  the  radial  mesh  density, 
the  global  mass  flow,  and  the  turbulence  model.  A  high 
flow  deflection  of  the  hub  boundary  layer  is  a  possible 
explanation  that  would  be  consistent  with  an 
underestimation  of  the  comer  stall. 

4)  All  the  turbulence  models  are  unable  to  predict  the 
strong  flow  deceleration  in  the  outer  wall  boundary 
layer.  Tire  high  loss  seems  to  appear  as  a  consequence 
of  the  interaction  between  the  passage  shock  and  tire 
leakage  flow  “vortex”.  The  flow  reaction  is  probably 
too  abrupt  across  tire  shock,  which  is  a  common  feature 
of  all  equilibrium  turbulence  models,  such  as  the 
mixing  length  or  linear  k-e  models.  It  is  important  to 
notice  tire  low  value  of  tire  turbulent  viscosity  near  tire 
tip  wall,  compared  with  tire  value  obtained  near  70- 
80%  of  tire  span  at  tire  trailing  edge  (Fig. 3. 23). 

5)  Tire  mixing  length  model  generates  more  loss  in  tire  tip 
wall  region  than  tire  k-E  model.  Tire  interaction  between 
tire  leakage  flow  and  tire  primary  flow  is  certaiirly  far 
more  complex  than  a  simple  boundary  layer  situation 


lor  which  tire  mixing  length  model  was  developed. 
Various  flow  phenomena  have  been  observed:  tire  shear 
layer  linked  with  the  suction  side  comer  at  the  blade 
tip,  tire  3D  separation  line  followed  by  tire  so-called 
leakage  “vortex”,  and  tire  radial  transport  of  tire  suction 
side  boundary  layer  followed  by  its  abrupt  transfer 
toward  tire  pressure  side  under  tire  influence  of  tire 
leakage  flow.  All  these  phenomena  have  to  be  treated 
along  the  tip  wall,  for  which  the  turbulence  model  only 
“sees”  a  single  length  and  a  single  velocity  scale,  while 
tire  tip  wall  flow  is  dominated  by  several  different 
lengths  and  velocity  scales. 

Solutions  usiirg  Baldwin-Lomax  turbulence  model  are 
mostly  similar,  despite  using  a  wide  range  of  different 
codes,  grid  topologies  and  number  of  points,  as  long  as 
at  least  300K  points  are  used. 


3.3  THE  RESULTS  FOR  THE  DLR  CASCADE 
3.3.1  Computations  requested 

Tire  contributors  were  requested  to  compute  tire  flow  field 
through  the  annular  turbine  cascade  corresponding  to  the 
operating  point  defined  in  Section  2.2.5.  Measurements  were 
taken  by  means  of  both  L2F  anemonretry  and  five-hole 
pressure  probes  in  a  plane  (MP3)  located  at  40%  of  an  axial 
chord  downstream  of  tire  nrid-span  trailing  edge.  For  the 
purpose  of  comparison  tire  pressure  probes  data  have  been 
used  for  pressures,  but  the  flow  angles  and  Mach  numbers 
measured  by  the  anemometer  have  been  preferred  to  those 
measured  by  tire  pressure  probes.  Tire  main  requested  plots 
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Fig  3.24  DLR  cascade:  Mach  number 
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(a)  algebraic  /  mixing  length  turbulence  models 
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(b)  turbulent  transport  models 


Fig  3  25  DLR  cascade  total  pressure  ratio 
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(a)  algebraic  /  mixing  length  turbulence  models 
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(b)  turbulent  transport  models 


Fig  3.26  DLR  cascade:  static  pressure  /  reference  total  pressure 
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(a)  algebraic  /  mixing  length  turbulence  models 


(b)  turbulent  transport  models 


Fig  3.27  DLR  cascade:  circumferential  flow  angle  (measured  from  the  tangential  direction) 
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were  the  spanwise  distributions  of  pitchwise-averaged  static 
pressure,  total  pressure  and  flow  angles,  and  the  contours  of 
the  same  quantities  in  the  measuring  plane.  In  addition,  the 
isolines  of  Mach  number  and  total  pressure  were  requested  in 
Hie  blade-to-blade  surfaces  at  10%,  50%  and  90%  of  the 
span.  In  order  to  assess  tire  behaviour  of  the  turbulence 
models,  contributors  had  to  provide  isolines  of  the  ratio  of 
turbulent  to  laminar  viscosity  |Vp  in  the  above  mentioned 
surfaces  as  well  as  a  pitchwise  Mach  number  traverse  at 
midspan  across  the  wake. 

When  averaging  the  computational  results,  pressures  and 
velocity  components  were  mass-averaged,  and  the  angles 
derived  from  the  averaged  velocities.  A  specific  check 
showed  that  when  the  measured  angles  were  averaged  in  the 
same  way  (instead  of  area-averaged)  the  changes  from  the 
reported  results  were  less  than  0.1°. 

Figs  3.24-27  compare  the  results  of  all  the  computations  with 
the  pitchwise-mean  measurements. 

3.3.2  Flow  description 

The  following  flow  features  can  be  deduced  from  the 
experimental  data  in  the  measuring  plane  (Figs  2.10-2.32): 

1)  The  radial  distribution  of  static  pressure  is  in  agreement 
with  the  tangential  flow  angles. 

2)  Moderate  negative  values  of  radial  flow  angle  (i.e.  flow 
directed  toward  the  hub)  can  be  found  in  the  wake  region,  as 
expected  in  an  annular  stator  cascade. 

3)  A  core  of  total  pressure  loss,  associated  with  the  secondary 
flow  vortices,  is  clearly  visible  near  the  hub.  In  the  tip  region 
the  loss  core  is  shifted  toward  the  midspan  and  there  is 
evidently  a  marked  distortion  of  the  wake. 

3.3.3  Overall  performance 

Few  contributors  specified  the  computed  values  of  mass  flow 
rate  and  overall  total  pressure  loss:  the  flow’  agrees  with 
orifice-measured  mass  flow  rate  within  1%,  but  the  loss 
varies  very  widely.  The  calculated  spanwise  distribution  of 
total  pressure  loss  is  reported  by  all  authors.  From  Fig  3.25,  it 
can  be  argued  that  the  overall  loss  prediction  disagrees  in 
some  cases  up  to  40%  with  the  measured  value.  It  can  also  be 
seen  that  the  local  discrepancies  at  each  spanwise  position 
are  on  average  much  higher. 

3.3.4  Radial  plots 

3.3.4. 1  Static  pressure  (Fig  3.26) 

The  slope  predicted  by  the  various  calculations  is  in  good 
agreement  with  the  experiments.  This  means  that  the 
condition  on  the  radial  distribution  of  pressure  imposed  at  the 
outlet  boundary  (usually  a  simplified  form  of  radial 
equilibrium)  does  not  impair  the  pressure  distribution  at  tire 
measuring  station.  Differences  among  the  values  at  fixed 
radial  positions  are  mainly  due  to  the  different  values  of 
pressure  assumed  by  the  contributors  at  the  downstream 
reference  radius.  More  surprising  is  the  influence  of 
turbulence  model  and  of  grid  topology  on  the  pressure 
distribution  computed  with  the  same  flow  solver  (TRACE-S). 

3. 3. 4. 2  Total  pressure  (Fig  3.25) 

Most  of  the  codes  display  total  pressure  radial  distributions 
grouped  within  a  reasonable  band  around  the  experimental 
values,  but  not  one  of  them  proved  able  to  capture  fairly  well 
all  of  the  significant  features  of  the  measurements.  Generally 
the  loss  core  at  the  hub  is  better  captured  than  the  one  at  the 
tip.  Looking  at  the  simulations  characterized  by  a  small 
number  of  points  in  the  radial  direction,  it  seems  that 
predictions  improve  with  a  radial  refinement  of  tire  grid. 


Fig  3.28  Static  pressure  at  MP3  /  reference  total 
pressure  (Bassi/Savini) 


However,  one  of  the  best  results  (TRANSCode)  does  not 
employ  a  great  number  of  points.  Once  again  the  influence  of 
turbulence  model  and  grid  topology  on  the  total  pressure 
distribution  predicted  by  TRACE-S  is  remarkable. 

3. 3.4. 3  Mach  number  (Fig  3.24) 

These  graphs  are  only  a  combination  of  the  pressure  and  total 
pressure  ones  and  their  general  trend  is  satisfactory, 
differences  worth  mentioning  occurring  only  in  the  near  wall 
regions. 

3.3.4.4  Circumferential  flow  angle  (Fig  3.27) 

The  scatter  in  the  computed  solutions  is  1-2°  and  the 
differences  between  measured  and  computed  tangential  flow 
angles  near  the  endwalls  are  significant.  All  codes  fail  to 
reproduce  the  correct  amount  of  overturning  at  the  tip. 
Computations  using  the  greatest  number  of  points  in  the 
spanwise  direction  exhibit  the  best  agreement  and  this 
observation  suggests  that  it  is  necessary  to  use  a  large  number 
of  points  to  obtain  a  correct  resolution  of  the  vortical 
structures  near  the  endwalls.  Maybe  an  even  greater  number 
of  points  is  needed  to  reach  a  truly  grid-independent  solution. 

3.3.5  Contour  plots  in  the  measuring  plane 

3.3.5. 1  Static  pressure 

Every  submitted  plot,  for  example  Fig  3.28,  shows  the  same 
tangential  wave-like  oscillations  of  pressure  displayed  by  the 
post-processing  of  the  measured  data.  These  are  an  outcome 
of  the  real  3D  blade  configuration  and  they  are  gradually 
damped  moving  downstream.  Notice  that  the  assumption  of 
pitchwise  uniform  static  pressure  at  the  outlet  boundary  made 
by  some  contributors  did  not  spoil  the  solution  at  the 
measuring  plane. 

3. 3. 5. 2  Total  pressure  (Fig  3.29) 

The  relevant  feature  of  the  comparison  among  the  various 
authors  and  the  experiments  is  tire  more  or  less  pronounced 
distortion  of  the  wake.  Concerning  this  item,  DLR’s  TRACE- 
S  computations  (performed  with  the  greatest  number  of 
spanwise  points)  gave  the  best  results  (b).  Nevertheless, 
approaching  the  tip,  the  computed  results  show  a  radial  band 


(e)  TRANSCode  using  B-L  turbulence  model  (Stapleton) 


(g)  using  k-ro  turbulence  model  (Bassi/Savini) 


Fig  3.29  Contours  of  total 


i)  TRACE-S  using  k-ro  turbulence  model  (Lisiewicz) 


I)  HAH3D  using  k-e  turbulence  model  (Hah) 
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(f)  TIP3D  using  mixing  length  model  (Denton) 


(h)  C  ANAR  I  using  k-i:  turbulence  model  (Couaillier) 


total  pressure  at  MP3 


Fig  3.30  Contours  of  circumferential  flow  angle  at  MP3 


(c)  TRACE-S  using  extended  k-e  model  (Lisiewicz) 


(a)  TRACE-S  using  k-o)  turbulence  model  (Lisiewicz) 


(c)  using  k-w  turbulence  model  (Bassi/Savini) 


(b)  TRACE-S  using  extended  k-e  model  (Lisiewicz) 


(d)  CANARI  using  k-e  turbulence  model  (Couaillier) 


(e)  CANARI  using  Michel  turbulence  model  (Couaillier) 


Fig  3.31  Contours  of  eddy  viscosity  ratio  at  MP3 
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with  almost  no  losses.  Good  agreement  with  the  experiments 
were  also  obtained  by  HAH3D  and  DRA  codes  (d,  e)  with  a 
somewhat  medium  radial  number  of  points.  No  grid 
refinement  studies  were  reported  for  TRACE-S,  so  it  is  not 
clear  whether  117  points  in  the  radial  direction  are  either 
necessary  or  sufficient.  Grid  refinement  studies  by  Stapleton 
did  not  show  significant  changes  for  the  DLR  turbine  for 
>330K  points  (49  radial).  McNulty  tried  grids  with  25  and  49 
radial  points,  and  his  results  clearly  show  that  25  radial 
points  are  insufficient,  but  not  how  many  more  are  needed. 
Examining  the  influence  of  turbulence  models,  the  TRACE-S 
results  using  the  k-co  and  the  extended  k-s  models  are  nearly 
equal  (b,  c).  Moreover,  Denton  (mixing  length),  Bassi-Savini 
(k-to)  and  Onera’s  CANARI  (k-e)  (f,  g,  h)  showed  quite 
similar  results,  thereby  suggesting  that  tire  turbulence  model 
influences  tire  mean  total  pressure  field  less  than  tire  grid 
radial  spacing;  but  Stapleton  found  tire  reverse:  a  change  in 
the  turbulence  model  from  Baldwin-Lomax  to  Spalart- 
AlLnraras  made  a  large  difference. 

3. 3. 5. 3  Radial  flow  angle 

No  great  differences  were  displayed  and  all  codes  predicted 
negative  flow  angles  in  tire  wake  region  and  positive  flow 
angles  in  tire  mainstream. 

3. 3. 5. 4  Circumferential  flow  angle  (Fig  3.30) 

The  plots  submitted  by  tire  various  contributors  support  tire 
observation  that,  in  general,  fine  grids  are  necessary  to 
resolve  the  details  of  tire  vortical  structure  near  tire  hub.  Tire 
TRACE-S  computation  performs  well  hr  resolving  such 
details,  but  even  finer  grids  would  be  advisable.  Tire  flow 
angle  distribution  hr  tire  relatively  larger  vortical  structure 
near  tire  tip  region  is  better  described,  especially  by  tire  fine 
grid  computations.  Tire  TRACE-S  computations  on  tire  same 
grid  with  the  extended  k-e  and  tire  k-co  turbulence  models  do 
not  show  significant  differences  in  tire  circumferential  flow 
angle. 

3. 3.5.5  Fig  3.31) 

These  contour  plots  reveal  many  surprising  facts.  Tire 
TRACE-S  turbulent  viscosity  computed  with  different  two- 
equation  models  (a,  b)  differs  by  more  than  one  order  of 
magnitude  (even  hr  tire  mainstream  region)  and  tire  same 
observation  applies  to  tire  TRACE-S  and  Bassi-Savini  results 
obtained  with  the  same  model  (k-co)  (a,  c).  Also  CANARTs 
and  TRACE-S's  k-s  results  are  not  similar  (b,  d),  leading  the 
fomrer  to  much  stronger  gradients  of  turbulent  viscosity  hr 
the  wake.  Finally,  tire  use  of  the  Michel  algebraic  model 
(CANARI)  led  to  a  complete  disappearance  of  turbulent 
viscosity  in  tire  measuring  plane  apart  from  tire  endwalls  (e). 
From  these  results  it  is  hard  to  extract  any  meaningful 
comment  on  tire  turbulence  models  and  it  is  felt  that  such 
strange  behaviour  is  mainly  due  to  numerical  aspects  of  tire 
codes  employed.  Further  studies  of  tire  implementation  of 
turbulence  models  hr  RANS  codes  seem  to  be  needed. 

3.3.6  Contour  plots  in  the  blade-to-blade  planes 
Contours  of  tire  TRACE-S  solution  on  tire  O-H-Fl-H  grid 
show  unphysical  behaviour  at  tire  downstream  junctions  of 
grid  subdomains  and  approaching  tire  outflow  boundary;  this 
is  dramatic  for  tire  case  of  Pt,  shown  hr  Fig  3.32,  and  pt/p 
isolhres.  This  happens  because  tire  grid  is  very  coarse  and 
diffusive  in  that  region.  It  illustrates  tire  difficulties  inherent 
in  implementing  complex  grid  schemes,  and  why  some  code 
writers  avoid  them. 

All  other  solutions  show  results  qualitatively  similar  to  each 


Fig  3.32  Total  pressure  ratio  contours  at  mid-span: 
TRACE-S  with  O-H-H-H  grid  and  k-e  model  (Lisiewicz) 


Fig  3.33  Viscosity  ratio  at  mid-span:  CANARI  using  two 
different  turbulence  models  (Couaillier) 


Mach  number 
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Fig  3.34  Mid-span  circumferential  Mach  number 
traverses 

other.  No  flow  separation  is  evident  on  the  blade  rear  suction 
surface  and  the  highest  losses  occur  in  the  near  hub  sections. 
This  is  obviously  to  be  ascribed  to  the  highest  freestream 
velocities.  Results  concerning  turbulence  characteristics  are 
similar  to  those  discussed  in  tire  previous  subsection.  It  is 
only  worth  noting  how  tire  turbulence  model  affects  tire 
turbulent  viscosity  in  the  blade  channel  obtained  with 
CANARI  code.  Using  tire  Michel  model,  (Vp.  values  are 
everywhere  low  (not  exceeding  100  even  close  to  the  profile) 
(Fig  3.33a)  whilst,  using  the  k-e  model,  values  are  greater 
(Fig  3.33b)  and  there  is  also  a  marked  turbulence  production 
in  the  suction  side  acceleration  region. 

Tire  only  two  tangential  midspan  Mach  traverses  submitted 
present  opposite  characteristics:  the  extent  of  the  wake,  that 
can  be  associated  with  the  region  of  deficit  in  the  Mach 
number  profile,  covers  the  whole  pitch  in  the  TRANSCode 
solution  and  about- 40%  of  tire  pitch  in  the  Bassi-Savini 
solution  (Fig  3.34);  the  peak  value  of  deficit  in  tire  latter  is 
twice  as  big  as  that  in  the  former.  (TRANSCode  used  49 
points  across  the  pitch  and  Bassi-Savini  37  points.) 

3.3.7  Assessment 

1)  The  computations  performed  clearly  pointed  out  tire 
importance  of  an  adequate  radial  resolution  of  tire  mesh.  This 
is  of  great  concern  particularly  outside  the  endwall  boundary 
layers  as  they  are  normally  well  enough  resolved.  Moreover, 
tire  computations  stressed  the  usefulness  of  a  fully-3D 
numerical  approach  when  simulating  such  a  flow,  it  is  felt 
that  most  of  the  relevant  3D  characteristics  cannot  be 
simulated  with  any  type  of  Quasi-3D  model. 

2)  Global  features  are  much  better  solved  than  local  features, 
especially  in  the  near  endwall  regions  where  tire  3D  effects 
induced  by  the  secondary  flows  are  significant.  One  way  to 
reduce  this  shortcoming  could  be  tire  use  of  solution-adaptive 
embedding.  Anyway  these  computations,  representative  of  tire 
state  of  the  art,  suggest  that  nearly  500,000  grid  points, 
adequately  distributed,  should  be  sufficient  to  solve  the 
steady  flow  field  in  an  annular  turbine  cascade. 

3)  Turbulence  models  are  often  considered  to  be  the  main 
source  of  discrepancy  in  the  validation  of  CFD  codes.  This  is 
not  always  true,  as  in  tire  present  case,  when  grid  resolution 


and  boundary  conditions  also  seem  to  play  a  major  role.  Big 
differences  are  still  hidden  in  the  numerics  rather  than  in  the 
physical  models  employed;  the  differences  between  tire 
TRACE-S  solutions  with  CH  and  OHHH  grids  suggest  that 
difficulties  with  multiblock  grids  can  occur. 

4)  The  huge  scatter  in  the  predicted  behaviour  of  turbulence 
quantities  is  surprising.  This  is  believed  to  be  actually  more 
related  to  the  practical  implementation  of  the  turbulence 
models  (e.g.  evaluation  of  source  terms,  clipping  of  maxima 
and  minima,  near-wall  treatment,  boundary  conditions,  etc.) 
than  to  the  physical  and  mathematical  aspects  of  tire  models 
themselves. 

Notwithstanding  these  uncertainities,  turbulence  models 
seem  to  have,  in  this  test  case,  a  moderate  influence  on  the 
aerodynamic  mean  flow  quantities;  however  it  can  be  argued 
that  much  stronger  effects  should  be  visible  analysing  the 
mixing  process  and  the  wall  heat  transfer  coefficients. 

3.4  ANALYSIS  OF  THE  COMPUTATIONAL  GRIDS 

3.4.1  Introduction 

Turbomachinery  CFD  code  developers  have  found  that  grid 
characteristics  can  have  significant  effects  on  the  accuracy  of 
a  computation.  Specifically,  those  grid  characteristics  known 
to  have  an  impact  on  the  solution  are: 

•  grid  type 

•  grid  size 

•  near-wall  characteristics,  including  normal  spacing  and 
cell  aspect  ratio 

•  grid  distortion  parameters,  including  stretching  and 
skewness 

•  tip  clearance  treatment. 

Discussions  of  these  parameters,  and  observations  of  their 
impact  on  the  accuracy  of  the  simulations  are  provided  in  the 
sections  that  follow.  Reference  is  made  to  Table  3.2  which 
summarizes  the  grid  information  for  each  of  the  simulations 
submitted  by  the  WG  26  contributors  for  the  Rotor  37  and 
DLR  cascade  test  cases.  The  table  presents  the  grid  type,  grid 
size,  y+  normal  to  the  wall  spacing,  minimum  and  maximum 
spacings  in  each  of  the  coordinate  directions,  near-surface 
maximum  cell-to-cell  spacing  ratios  (clustering),  numbers  of 
leading  and  trailing  edge  points,  points  in  the  spanwise 
direction  in  the  tip  gap  region,  maximum  cell  aspect  ratios, 
maximum  skew,  number  of  blocks  and  minimum  and 
maximum  values  of  the  axial  coordinate  for  the  solution 
domain. 

3.4.2  Grid  Type 

Brief  descriptions  of  tire  grid  systems  used  by  the  WG  26 
contributors  are  given  in  the  ensuing  subsections.  In  all  cases, 
the  grids  were  of  the  non-adaptive  structured  type.  They  are 
denoted  by  type  as  H,  I,  C  and  O  grids.  Many  contributors 
used  composites  of  these  grids  which  were  coupled  along 
contiguous  boundaries  and  others  used  overset  grid 
combinations.  All  of  these  3-D  grid  systems  were  constructed 
by  stacking  2-D  grids  generated  on  blade-to-blade  surfaces  of 
revolution  with  tire  hub  and  tip  grids  conforming  to  the 
flowpatlr  boundaries.  It  should  be  noted  that  composite  grids 
generally  reduce  the  undesirable  features  associated  with 
grids  of  a  single  grid  type  but  that  they  introduce  additional 
complexity  in  tire  flow  solver  along  tire  component  grid 
boundaries  and  can  lead  to  non-physical  features  in  the 
solutions  in  these  areas. 
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Fig  3.36  Rotor  37  H-l-H  grid  (Hildebrandt) 


3.4.2. 1  HGrids 

An  H  grid  is  the  most  widely  used  grid  type  for 
turbomachinery  flow  calculations  and  it  is  preferred  by  most 
researchers  because  of  ease  of  construction,  hi  its  simplest 
form,  it  is  generated  with  “straight”  line  segments  along  the 
x,  r  and  6  coordinate  lines  using  algebraic  techniques.  To 
resolve  high  gradient  regions  near  the  airfoil  surfaces,  hub 
and  casing,  and  airfoil  leading  and  trailing  edges,  clustering 
laws  in  the  form  of  exponential  functions  are  employed. 
Generalizations  of  the'  H  grid  with  curved  grid  lines  which 
reduce  shear  and  increase  airfoil  leading  and  trailing  edge 
resolution  have  been  used  by  some  researchers.  The  primary 
drawback  with  H  grids  is  the  high  grid  shear  inherent  with 
compressor  airfoils  at  high  setting  angles  and  turbine  airfoils 
with  high  turning.  Another  disadvantage  is  that  grid 
independent  solutions  are  difficult  to  achieve  because  as  the 
grid  is  refined,  the  near-leading-edge  and  trailing-edge  grid 
cells  become  more  highly  sheared  which  introduces  large 
metric  gradients  and  additional  discretization  errors  into  the 
solution.  As  a  result  of  this  leading-  and  trailing-edge 
behavior,  highly  accurate  H  grid  simulations  are  generally 
only  achieved  by  very  experienced  users  who  have  extensive 
background  in  generating  acceptable  H  grids.  The  H  grid  is 
generally  better  suited  for  compressor  computations  than 
turbine  computations  because  the  compressor  airfoils  are 
normally  tliinner  and  produce  less  shear  around  lire  airfoil 
leading  and  trailing  edges.  Calvert,  Denton,  Shabbir  and 
Hutchinson  used  this  grid  type  for  Rotor  37  simulations  and 
Stapleton  and  Denton  used  it  as  well  for  the  DLR  cascade  test 
case.  The  Rotor  37  H  grid  employed  by  Calvert  is  shown  in 
Fig  3.35. 

3.4. 2. 2  I  Grids 

I  grids  were  introduced  to  reduce  the  undesirable  grid  shear 
inherent  with  H  grids.  They  are  normally  constructed  with 
algebraic  techniques  which  align  the  blade-to-blade  grid  lines 
in  file  direction  normal  to  the  stream  wise  grid  lines.  I  grids 


Fig  3.37  DLR  cascade  H-C  grid  (Lisiewicz) 

substantially  reduce  grid  shear  but  introduce  complications 
for  the  flow  solver  along  the  grid  periodic  boundaries  where 
the  grid  lines  are  generally  not  continuous  and  not  periodic. 
An  additional  difficulty  is  introduced  at  the  upstream  and 
downstream  boundaries  in  the  specification  of  boundary 
conditions  because  those  boundaries  are  skewed  relative  to 
tlie  axial  planes  utilized  with  H-type  grids.  Hah  employed  I 
grids  for  his  DLR  cascade  and  Rotor  37  calculations. 
Hildebrandt  devised  a  novel  periodic  composite  H-I-H  grid 
for  Rotor  37,  shown  in  Fig  3.36,  which  had  none  of  the 
undesirable  features  normally  associated  with  I  grids,  being 
periodic,  continuous  along  the  periodic  boundaries  and 
having  axial  grid  planes  at  the  upstream  and  downstream 
boundaries  for  ease  of  specification  of  inflow  and  outflow 
boundary  conditions. 

3.4.2. 3  C  Grids 

C  grids  wrap  around  the  airfoil  leading  edge  in  a  “C”  pattern 
and  have  the  desirable  feature  of  providing  a  blade 
conforming  orthogonal  grid  in  that  region  for  improved 
resolution  of  high  leading  edge  flow  gradients.  Similar  to  H 
grids  downstream,  they  align  with  file  flow  direction  in  that 
region  for  good  wake  definition.  Disadvantages  of  C  grids 
are:  (1)  they  are  limited  in  file  distance  which  they  can  be 
extended  upstream  of  fire  airfoil  row,  and  (2)  they  become 
more  sparse  with  distance  upstream  -  which  can  be  a  concern 
for  compressors  operating  at  supersonic  inlet  conditions  with 
respect  to  resolving  shock  strong  waves  propagating 
upstream.  C-type  grids  are  normally  constructed  with  elliptic 
grid  generation  techniques.  Periodic  C  grids  were  employed 
by  Bassi/Savini,  Dadone/DePalma,  and  Chima  for  the  DLR 
cascade.  Non-periodic  C  grids  with  non-contiguous  points 
along  file  wake  cut  were  used  by  McNulty  for  the  DLR 
cascade  and  by  Amone  for  Rotor  37.  The  block  composite 
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Degenerate  0-Grid 


Fig  3.38  Rotor  37  H-O-H  grid  (McNulty) 


periodic  H-C  grid  shown  in  Fig  3.37  was  used  by  Lisievvicz 
for  the  DLR  cascade  with  an  H  grid  extension  upstream. 

3.4.2. 4  0  Grids 

O  grids  wrap  around  the  airfoil  and  have  one  family  of  grid 
lines  conforming  to  the  airfoil  surface.  With  sufficient  grid 
density,  they  can  be  used  to  provide  high  resolution  of  both 
leading-  and  trailing-edge  flows.  Typically,  O  grids  are 
orthogonal  or  nearly  orthogonal  along  the  airfoil  surface, 
thereby  minimizing  tire  undesirable  grid  shear  inherent  with 
H  type  grids.  They  share  the  same  disadvantages  of  C  grids  in 
that  they  are  limited  in  the  distance  they  can  be  extended 
upstream  and  downstream  and  they  become  sparse  with 
distance  upstream  and  downstream.  There  is  an  additional 
concern  with  the  latter  feature  downstream  as  it  is  desirable 
to  maintain  tight  clustering  in  the  wake  region.  Some 
researchers  have  tried  to  circumvent  this  problem  by  moving 
the  O  grid  downstream  boundary  closer  to  the  airfoil  to 
maintain  a  high  density  mesh  locally  in  the  trailing  edge 
region  and  adding  a  high  density  H  grid  downstream.  Others 
have  used  O  grids  only  in  the  local  region  of  the  airfoil  and 
coupled  separate  H  grids  upstream,  downstream  and  within 
the  airfoil  passage.  Like  C  grids,  O  grids  are  normally 
constructed  with  elliptic  grid  generation  schemes.  It  is 
generally  accepted  that  O  grids  are  best  suited  for  turbine 


airfoils  because  they  can  be  used  to  define  large  trailing 
edges  accurately  as  is  needed  to  predict  the  airfoil  base 
pressure.  Couaillier  used  non-periodic  H-O-H  grids  for  the 
Rotor  37  and  DLR  test  cases  and  for  Rotor  37  he  also 
investigated  a  local  O  grid  coupled  with  a  background  H  grid 
along  a  contiguous  patched  grid  boundary.  McNulty 
employed  the  H-O-H  grid  and  the  degenerate  O  grid  in  the  tip 
clearance  shown  in  Fig  3.38  for  Rotor  37.  Lisiewicz 
performed  a  calculation  of  the  DLR  cascade  with  a  novel 
periodic  O-H-H-H  grid  consisting  of  4  separate  grids  coupled 
along  patch  boundaries.  Weber  used  an  overset  O-H  grid, 
shown  in  Fig  3.39,  which  did  not  suffer  from  the  undesirable 
O-grid  “comers”  in  patched  O-H  grid  systems.  With  this 
system,  the  grid  cell  sizes  in  the  O  and  H  grids  in  the  overset 
region  are  maintained  approximately  the  same  to  ensure 
accurate  transfer  of  boundary  data  between  the  grids  by 
interpolation. 

3. 4. 2. 5  Conclusions  on  Grid  Type 

hi  tlie  computations  for  Rotor  37  and  the  DLR  cascade 
undertaken  for  WG26,  no  single  grid  type  stood  out  as  being 
superior  to  the  others.  In  general,  the  O  and  C  grids  provided 
the  best  definition  of  the  airfoil  leading  edge  and  gave  the 
best  resolution  of  the  leading  edge  flow  gradients  and  the 
bow  shock  structure  for  Rotor  37.  One  exception  to  this 
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Fig  3.39  Rotor  37  overset  O-H  grid  (Weber) 


generalization  was  the  Rotor  37  simulation  of  Hildebrandt 
performed  on  the  H-I-H  grid  shown  in  Fig  3.36  with  only  4 
points  defining  the  leading  edge.  In  this  case,  excellent  shock 
resolution,  shown  in  Fig  3.40,  was  achieved  by  aligning  the 
grid  lines  in  the  blade-to-blade  surfaces  along  and  normal  to 
the  mean  flow  directions,  thereby  aligning  the  grid  lines 
parallel  and  normal  to  the  passage/bow  shock  system.  With 
the  exception  of  the  near  hub  region,  Hildebrandt's  prediction 
showed  remarkably  good  agreement  with  the  experimental 
spanwise  profiles  of  total  pressure,  total  temperature  and 
adiabatic  efficiency.  His  overall  total  pressure  predictions 
also  showed  good  agreement  with  the  measurements  over  the 
entire  speed  line.  It  should  be  noted,  however,  that  the 
excellent  results  achieved  in  this  case  were  likely  the  result 
of  aligning  the  grid  lines  with  the  dominant  shock  features 
and,  therefore,  that  it  may  be  difficult  to  achieve  comparable 
levels  of  accuracy  for  all  airfoil  shapes  and  flow  conditions 
with  this  grid  system. 

In  a  related  investigation  on  leading  edge  flow  resolution, 
Amone  performed  a  Q3D  study  of  the  effect  of  leading  edge 
definition  for  Rotor  37  with  a  C  grid  in  which  he  varied  tire 
number  of  points  around  the  leading-edge  circle  from  4  to  40. 
He  determined  that  the  predicted  total  pressure  rise 
characteristic  changed  monotonically  with  increased 
resolution  and  converged  to  a  single  characteristic  with  20 
points  around  the  leading  edge.  However,  the  necessary 
number  of  points  must  surely  depend  on  the  type  of  leading 
edge  flow,  the  transition  assumption  and  the  turbulence 
model  used. 

Unfortunately,  neither  the  Rotor  37  nor  the  DLR  cascade  test 
cases  provided  airfoil  surface  data  needed  to  determine  the 
effect  of  leading  edge  definition  on  tire  airfoil  surface  flow 
conditions  properly.  It  is  well  known  that  any  errors  incurred 
at  the  leading  edge  show  up  as  entropy  increases  which 
convect  back  onto  the  airfoil  surfaces  and  adversely  affect  tire 
quality  of  the  solution  over  the  entire  surface.  Additionally, 
any  errors  on  the  pressure  and  suction  surfaces  generated  at 
the  leading  edge  combine  downstream  to  affect  the  wake 
prediction.  Because  of  the  potential  significance  of  the 
leading  edge  flow  modeling  on  tire  accuracy  of  the  cascade 


Fig  3.40  TRACE-S  Rotor  37  relative  Mach  contours 
at  90%  span  and  98%  choke  flow  (Hildebrandt) 


flow  simulations,  further  experimental/computational  studies 
are  recommended.  For  these  investigations,  2-D  cascade 
experiments  and  computations  are  recommended  in  order  to 
eliminate  any  3-D  effects  which  could  mask  leading  edge 
effects. 

Weber's  remarkably  good  results  for  Rotor  37  may  be 
attributable,  at  least  in  part,  to  tire  use  of  the  overset  O-H  grid 
system  shown  in  Fig  3.39  .  This  grid  does  not  suffer  from  the 
undesirable  sharp  O-grid  “comers”  inherent  with  patched  O- 
H  grid  systems  where  metrics  gradients  are  large  and  the 
grids  are  sparse.  With  the  overset  O-H  grid  system,  tire  grid 
cell  sizes  in  tire  O  and  H  grids  in  tire  overset  region  are 
maintained  at  nearly  the  same  levels  for  good  interpolation 
data  transfer  between  grids  and  the  metric  variations  in  both 
grids  are  smooth  throughout.  Also  on  tire  endwall,  the  use  of 
the  “background”  H  grid  in  the  overset  O-H  grid  system 
ensures  smooth  development  of  the  endwall  boundary  layers 
entering  the  bladed  region,  which  may  be  a  significant  factor 
with  respect  to  predicting  the  hub  suction  surface  separation 
observed  with  Rotor  37.  The  primary  drawback  for  the  use  of 
tire  overset  O-H  system  is  that  file  flow  solver  is  non¬ 
conservative  as  a  result  of  the  required  interpolation 
procedures  used  to  transfer  data  between  file  grids.  For  the 
Rotor  37  test  case,  however,  these  procedures  appeared  to 
result  in  no  appreciable  negative  impact  on  file  solution. 

Regarding  the  airfoil  trailing  edge  flow  and  the  wake 
structure,  the  C,  I  and  H  grids  that  align  with  the  wake 
generally  provided  the  best  wake  definition.  These  grids  have 
one  family  of  grid  lines  aligned  with  file  wake  and  have  a 
large  number  of  points  across  file  wake.  Hie  O  grid  solutions, 
on  file  other  hand,  provide  high  resolution  of  the  local  trailing 
edge  wake  structure  but  rapidly  diffuse  the  gradients  as  file 
grid  opens  up  downstream  of  the  trailing  edge.  For  turbines, 
the  “best”  grid  for  trailing  edge  flow  and  wake  resolution  is 
considered  by  some  WG  26  members  to  be  the  H-O-H  grid 
which  can  be  used  to  provide  sufficient  local  trailing  edge 
definition  to  predict  file  airfoil  base  pressure  and  H-grid 
alignment  with  file  wake  downstream  to  resolve  file  wake 
structure,  though  that  conclusion  is  not  based  on  the  results  of 
file  DLR  cascade  case. 
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3.4.3  Grid  Size 

A  number  of  investigators  studied  tire  effect  of  grid  size  or 
grid  fineness  on  solution  accuracy.  In  general,  these  studies 
showed  that  a  moderate  grid  of  200,000  points  was  needed  to 
capture  overall  performance  characteristics  and  that  much 
finer  grids  with  as  many  as  1,000,000  may  be  needed  to 
isolate  the  detailed  flow  features  such  as  endwall  secondary 
flows  and  tip  clearance  flows.  If  wall  functions  are  used,  less 
grid  points  are  needed  near  the  walls,  so  less  points  are 
needed  overall. 

McNulty  studied  the  effect  of  grid  size  for  the  DLR  cascade 
test  case.  He  employed  three  C  grids  with  approximately 
160,000,  300,000,  and  1,200,000  points.  The  two  variations 
from  tire  baseline  160,000  point  grid:  (1)  doubled  tire  number 
of  radial  grid  points,  and  then  (2)  doubled  tire  number  of 
normals  along  the  airfoil  surface  and  C  contour  grid  lines, 
respectively.  McNulty  determined  that  at  least  49  grid  points, 
which  he  used  in  tire  radial  direction  with  tire  finer  meshes, 
were  needed  to  represent  the  endwall  secondary  flows  near 
tire  hub  and  casing.  This  study  also  showed  that  even  more 
points  were  needed  in  tire  radial  direction  to  fully  capture  tire 
endwall  secondary  flow  vortices.  The  DLR  cascade  solutions 
by  Lisiewicz  on  C-H  and  O-H-H-H  grids,  both  with  117 
points  in  tire  radial  direction,  confirmed  this  finding  and 
suggested  that  even  finer  grids  in  the  radial  direction  may  be 
needed  to  match  the  data. 

Stapleton  studied  tire  effect  of  grid  size  for  tire  DLR  cascade 
with  H  grids.  He  used  a  baseline  135x49x49  grid  with 
approximately  325,000  points  and  variants  with:  (1)  more 
than  twice  tire  number  of  points  (1 17)  in  the  radial  direction, 
and  (2)  24%  more  in  the  axial  and  49%  more  in  tire 
circumferential  and  radial  directions.  He  found  that  all  grids 
produced  nearly  the  same  pitchwise-average  performance  and 
that  tire  total  pressure  contours  for  all  of  tire  solutions 
qualitatively  matched  the  measured  contours,  but  that  the  fine 
radial  grid  solutions  gave  slightly  better  definition  of  the  near 
endwall  secondary  flow  structure. 

McNulty  also  looked  at  tire  effect  of  grid  size  on  solution 
accuracy  for  Rotor  37  using  H-O-H  grids  with  degenerate  O 
grids  in  the  tip  clearance  region.  Three  grids  with 
approximately  200,000,  400,000  and  350,000  points  were 
studied  to  look  at  the  effect  of  first  marginally  increasing 
(20%-40%)  tire  grid  size  in  each  direction  and  then  doubling 
tire  number  of  points  only  hr  the  radial  direction  from  the 
baseline  200,000  point  grid.  In  each  case,  nine  points  were 
used  in  tire  radial  direction  hr  tire  tip  clearance  region.  All  of 
these  solutions  gave  remarkably  similar  pitchwise-average 
performance  predictions. 

3.4.4  Near-Wall  Characteristics 

Those  grid  properties  locally  affecting  solution  accuracy  hr 
tire  region  of  tire  wall  shear  layers  are  tire  spacing  normal  to 
tire  wall  and  tire  cell  aspect  ratio.  A  general  nomral-to-the- 
wall  grid  spacing  guidelhre  for  calculations  that  solve  down 
to  tire  wall  and  enforce  no-slip  wall  boundary  conditions  is  to 
maintain  a  normal  spacing,  y+,  of  less  than  1 .0,  which  ensures 
that  the  first  pohrt  off  tire  wall  lies  within  tire  viscous 
sublayer  of  tire  boundary  layer.  However,  it  depends  on  tire 
turbulence  model  used.  None  of  the  contributors  who  solved 
down  to  tire  wall  reported  y+  values  less  than  1  throughout, 
and  only  a  few  maintained  y+  levels  less  than  10  at  every  wall 
point. 

For  those  computations  employing  wall  functions  without 
solvhrg  down  to  tire  wall,  tire  accepted  guidelhre  is  to  place 
tire  first  point  off  the  wall  within  tire  fully  turbulent  inner 


region  of  tire  boundary  layer  with  y*  values  of  the  order  of  50. 
Most  contributors  employing  wall  functions  adhered  to  this 
guideline  and  reported  y+  values  between  20  and  60. 

Cell  aspect  ratios  near  tire  airfoil  and  endwall  surfaces  can 
become  extremely  large  as  points  are  clustered  in  those 
regions.  This  is  especially  true  with  those  solvers  that  do  not 
employ  wall  functions,  but  rather  solve  down  to  the  wall. 
Typical  maximum  cell  aspect  ratios  reported  from  those 
contributors  not  employing  wall  functions  ranged  from 
approxhrrately  200  to  800.  In  general,  for  those  contributors 
employing  wall  functions,  tire  reported  maximum  cell  aspect 
ratios  were  less  than  100. 

3.4.5  Grid  Distortion  Parameters 

Grid  distortion,  which  can  have  an  adverse  effect  on  solution 
accuracy,  can  be  characterized  by  the  degree  of  grid 
stretching  in  any  coordinate  direction  and  the  degree  of  grid 
skewness  or  shear  between  any  two  coordinates. 

Normally,  turbomachinery  grids  are  clustered  toward  the 
walls  to  resolve  the  endwall  and  airfoil  surface  boundary 
layers  without  incurring  the  large  CPU  times  required  for 
calculating  the  flow  on  a  uniformly  distributed  dense  mesh. 
All  of  the  WG  26  contributors  utilized  grid  stretching  to 
cluster  points  near  Arose  surfaces.  This  approach  is 
acceptable  for  small  mesh  stretching  ratios  (cell-to-cell  size 
ratio);  however,  it  does  introduce  error  into  the  solution  that 
is  proportional  to  the  degree  of  stretching  for  a  given  code,  hi 
general,  it  is  accepted  that  the  cell-to-cell  size  ratio  should 
not  exceed  a  value  of  1.3.  Most  contributors  who  reported 
grid  clustering  information  stayed  within  this  guideline. 

Grid  skewness  is  known  to  have  adverse  effects  on  the 
accuracy  and  stability  of  a  numerical  solution.  That  is,  as  the 
skewness  or  shear  increases,  the  accuracy  decreases  and  the 
solution  becomes  less  stable.  For  the  purpose  of  this  study, 
maximum  skew,  as  reported  in  Table  3.2,  is  defined  as  the 
minimum  angle  between  the  surface  and  the  grid  line 
extending  away  from  the  surface.  In  all  cases,  it  was  assumed 
that  the  maximum  skewness  occurred  on  the  airfoil  surface. 
In  general,  the  H-  and  C-  type  grids  produced  the  most  grid 
shear  with  it  occurring  at  or  near  the  trailing  edge  of  the 
airfoil.  Stapleton  reported  a  minimum  skew  angle  of  about  4° 
for  his  H-grid  DLR  cascade  solutions  and  Calvert  21°  for  his 
H-grid  Rotor  37  solution.  Likewise,  Chima  reported  low 
skew  angles  of  18  and  30°  for  his  C-grid  DLR  cascade  and 
Rotor  37  solutions,  respectively.  I  grids  employed  for  Rotor 
37  computations  by  Hah  and  Plildebrandt  produced  much  less 
shear  with  reported  skew  angles  of  about  60°  for  each.  The  O 
grid  employed  by  Lisiewicz  for  the  DLR  cascade  had  no  shear 
with  a  reported  minimum  skew  angle  of  90°. 

Grid  distortion  tolerance  probably  depends  partly  on  the  code; 
for  example  cell  vertex  methods  may  be  more  tolerant  than 
cell  centred  schemes. 

3.4.6  Tip  Clearance  Treatment 

A  number  of  different  treatments  have  been  employed  to 
model  the  flow  through  the  tip  clearance  region  in  rotors. 
These  range  from  the  simplest  so-called  “pinched-tip”  model 
to  more  elaborate  fully  gridded  models.  The  pinched  grid 
model  treatments  reduce  the  airfoil  tip  to  zero  thickness  in 
one  mesh  spacing  and  employ  just  a  few  points  in  the 
clearance  region.  The  rationale  for  this  very  simple  model  is 
the  fact  that  the  tip  clearance  flow  is  dominated  by  inviscid 
effects  and  can  be  approximately  computed  knowing  the 
pressure  difference  across  the  blade  tip.  These  models  are 
easily  implemented  with  H  grids  where  the  grid  lines  on  tire 
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pressure  and  suction  surfaces  are  aligned  in  the 
circumferential  direction.  Normally  the  computational  tip  gap 
is  not  the  same  as  the  actual  gap  size  but  is  determined  by 
numerical  experimentation  in  order  to  allow  for  the  vena 
contracla.  Rotor  37  computations  by  Calvert,  Denton  and 
Shabbir  employed  the  pinched-tip  model.  Hutchinson 
performed  two  calculations  using  the  pinched-tip  treatment 
with  4  and  6  nodes  in  the  tip  gap. 

The  fully  gridded  tip  treatments  attempt  to  model  tire  actual 
airfoil  tip  shape  and  tire  clearance  region.  McNulty  and 
Couaillier  employed  degenerate  O  grids  in  the  tip  clearance 
for  Rotor  37  that  were  patched  along  tire  airfoil  surface  to  the 
passage  O  grids  along  contiguous  patch  grid  interfaces.  In 
these  cases,  about  10  grid  points  were  used  in  tire  radial 
direction.  Hildebrandt  studied  tire  effect  of  number  of  points 
in  a  fully  gridded  H  grid  treatment  in  tire  tip  gap  for  Rotor  37 
and  performed  calculations  with  4  and  7  points  in  the  radial 
direction. 

A  study  by  Clrirrra  (1996b)  showed  relatively  little  difference 
between  simple  and  sophisticated  approaches  to  tire  grid  in 
tire  tip  clearance  region. 

Because  of  tire  lack  of  detailed  flow  measurements  and 
computational  results  in  tire  tip  clearance  region  for  Rotor  37, 
it  was  not  possible  to  determine  the  precise  impact  of  the  tip 
clearance  modeling  on  the  flow  solutions.  It  was  generally 
felt  by  tire  Working  Group  members  that  tire  fully  gridded 
models  provided  a  better  basis  for  predicting  the  tip  clearance 
flow  than  tire  “pinched  tip”  models,  especially  for  predicting 
the  tip  clearance  losses.  However,  in  some  cases,  tire 
predictions  with  tire  “pinched  tip”  models  produced  better 
agreement  with  tire  experimental  spanwise  profiles  of  total 
pressure  and  efficiency  than  the  fully  gridded  models. 
Additional  experiments  and  supporting  computations  are 
needed  to  determine  the  impact  of  tire  tip  clearance  model  on 
the  global  simulation  more  clearly. 


3.5  ANALYSIS  OF  THE  ALGORITHMS 

In  Chapter  1 ,  tire  various  algorithms  used  for  turbonrachinery 
computations  were  reviewed.  They  fall  into  two  classes:  time 
marching  and  pressure  correction.  Within  tire  time  marching 
class,  used  by  most  contributors,  there  are  many  variations. 
Tire  turbonrachinery  CFD  specialists  have  generally  adopted 
algorithms  proposed  by  external  aerodynamics  research 
workers,  and  unlike  the  recent  Fluid  Dynamics  Panel  meeting 
(AGARD,  1995)  turbonrachinery  CFD  meetings  have  rarely 
discussed  this  aspect  of  their  work. 

Initially,  tire  time  marching  algorithms  were  only  applicable 
to  compressible  flows  with  Mach  number  above  about  0.3, 
and  that  limitation  applies  to  many  of  tire  codes  employed  for 
the  WG,  while  the  pressure  correction  codes  were  limited  to 
subsonic  flow.  More  recently,  these  limitations  have  been 
overcome,  but  that  has  not  made  specialists  switch  their 
chosen  approach. 

Tire  choice  of  algorithm  is  made  to  achieve  fast  and  rapidly 
convergent  computations,  and  tire  present  WG  has  focused  on 
accuracy  rather  than  speed  or  stability.  It  is  generally  believed 
that  tire  choice  of  algorithm  has  little  influence  on  accuracy, 
all  solutions  of  the  same  equations  converging  -  if  they 
converge  at  all  -  to  tire  same  results,  within  engineering 
accuracy. 

The  question  of  computational  speed  is  difficult  to  judge  in  a 
study  of  this  kind,  because  of  tire  wide  variety  of  computer 


hardware  used.  Tire  codes  were  often  run  on  a  network,  so 
that  tire  true  running  time  is  not  known.  Tire  WG  has  made 
no  attempt,  therefore,  to  rate  relative  speeds. 

From  tire  contributions  made  to  WG  26,  there  is  no  basis  for 
assessing  tire  relative  merits  of  tire  schemes  used  in  tire 
various  codes. 

3.6  ANALYSIS  OF  THE  TURBULENCE  MODELS 
It  is  quite  difficult  to  compare  tire  different  turbulence  models 
employed,  since  some  of  the  differences  between  the  various 
computations  of  tire  DLR  cascade  are  believed  to  depend  on 
tire  grids  and  flow  solvers  used  rather  than  on  the  turbulence 
models.  So  tire  turbulence  models  are  compared  for  that  case 
only  under  tire  same  conditions,  i.e.  when  using  the  same 
code  and  grid. 


3.6.1  Mixing  length  models 
Kang  and  Hirsch 


Code 

EURANUS/TURBO 

Test  case 

Rotor  37 

Model  description 

The  Baldwin-Lomax  (1978) 
algebraic  turbulence  model  is 
used.  All  the  constants  of  the 
model  are  tire  standard  ones 
except  Cwk=l  .0. 

Model  Implementation 

In  EURANUS/TURBO,  ym„x  is 

7  J  max 

the  value  of  y  corresponding  to 
tire' maximum  of  F. 

Transition  model 

No  transition  model.  Fully 
turbulent  flow  is  assumed. 

The  pitchwise  average  total  pressure  and  temperature  and 
flow  angle  exhibit  some  of  the  experimental  radial  variations, 
particularly  in  tire  hub  region.  Other  authors,  using  the  same 
Baldwin-Lomax  model,  with  similar  quantity  of  nresh  points 
have  produced  results  with  a  linear  variation  in  this  hub 
region.  According  to  tire  authors,  these  better  simulations 
may  be  associated  with  tire  model  implementation,  and  the 
yn,ax  definition.  They  suggest  that  other  similar  codes  take  tire 
ymax  corresponding  to  tire  F  maxima  furthest  away  from  a 
wall.  The  treatment  of  separated  wall  layers  could  be 
improved  for  this  reason. 

The  radial  distributions  agree  with  the  observations  presented 
in  paragraph  3.2.3.  Tire  ratio  between  the  computed  choked 
mass  flow  and  the  measured  one  is  0.993.  The  increase  of  tire 
computed  axial  velocity  at  mid-span  (R=0.22nr)  at  station  4 
can  be  estimated  with  the  following  observations: 

•  Tire  stagnation  temperature  is  correct 

•  Tire  absolute  flow  angle  is  underestimated  by  3° 

The  ratio  of  axial  velocity  is  then  VzSin,uiA,zeXpe=  1.113.  This 
means  that  a  large  amount  of  flow  blockage  is  generated  at 
the  tip  and  hub  walls  as  seen  in  Kang’s  plots  of  tire  efficiency 
and  entropy.  Note  also  that  the  mesh  is  coarse  for  station  4, 
so  that  some  numerical  dissipation  exists  between  tire  blade 
trailing  edge  and  station  4. 

A  large  amount  of  entropy  is  generated  along  tire  line  of 
interaction  between  tire  leakage  flow  and  the  primary  flow.  A 
second  zone  of  high  entropy  is  also  observed  on  tire  suction 
side,  after  the  passage  shock.  This  is  connected  with  a  strong 
local  aerodynamic  blockage  for  Vz  It  is  remarkable  that  tire 
turbulent  viscosity  has  a  medium  value  after  tire  impact 
between  the  3D  leakage  “vortex”  and  tire  passage  shock  (pt 
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/(j.=l 50-250).  Tlic  largest  value  of  the  turbulent  viscosity  is 
observed  further  from  the  tip  wall,  near  tire  suction  side,  as  a 
probable  consequence  of  the  separated  layer,  and  the  strong 
radial  migration  towards  the  tip  wall  (|it  /p-700-940). 

Note  also  that  preliminary  results  of  the  authors  with 
Cwk=0.25  (the  originally  proposed  value  of  Baldwin-Lomax) 
generate  a  large  separation  on  the  blade  suction  side. 

Chima 


Code 

SWIFT  -  multiblock 

Test  case 

Rotor  37 

Model  description 

The  simulation  neglects  the  viscous 
terms  in  the  streamwise  direction, 
according  to  a  thin  layer 
approximation.  The  3D  results  have 
been  obtained  with  the  Baldwin- 
Lomax  model  (1978).  All  the 
constants  of  tire  model  are  the 
standard  ones  except  C  =1.216  and 
=  0.646  which  were  shown  to 
give  a  better  agreement  with  the 
Cebeci-Smith  model  (Chima,  Giel, 
and  Boyle,  1993).  Tire  adapted 
model  predicts  the  stall  point  better 
than  tire  original  implementation 
(Chima,  1996a). 

2D  results  at  70%h  have  been 
obtained  also  with  the  k-co  Wilcox's 
turbulence  model  (1994). 

Model 

Tire  turbulence  model  is  adapted  in 

Implementation 

tire  clearance  region.  Tire  inner 
fonnulation  is  used  near  the  blade 
tip  and  tire  casing.  A  constant  outer 
turbulent  viscosity  is  used  across  the 
rest  of  the  gap.  For  the  outer 
formulation,  Fmax  is  taken  as  tire 
maximum  of  the  function  F  across 
tire  entire  gap,  and  Ymax  is  taken  as 
the  distance  to  the  nearest  of  tire  tip 
or  casing  walls.  In  comers,  tire 
distance  is  computed  with  Buleev's 
formulation  for  tire  inner  layer.  For 
tire  outer  layer,  the  classical  distance 
is  used. 

Boundary  conditions 

Integration  down  to  the  wall. 

Transition  model 

Tire  hub  and  casing  walls  are 
assumed  fully  turbulent.  Tire  blade 
boundary  layers  are  allowed  to 
undergo  transition  using  tire  model 
proposed  by  Baldwin-Lomax.  The 
turbulence  is  activated  if  |Vn>14. 

Special  wake 

Tire  wake  model  of  Baldwin-Lonrax 

treatment 

is  used. 

At  mid-span,  for  98%  of  the  choked  mass-flow,  the  radial 
plots  in  station  4  show  an  overestimate  of  the  stagnation 
temperature,  while  the  absolute  flow  angle  is  correct.  This 
suggests  that  the  predicted  axial  flow  velocity  is  too  high. 
This  is  consistent  with  a  large  amount  of  loss  near  the  tip 
wall,  that  produces  a  low  efficiency. 

The  gradients  of  stagnation  pressure  are  smoothed  out  in  the 
lower  50%  of  the  span.  The  good  overall  efficiency  results 
from  a  balance  between  an  overestimate  by  2%  at  mid-span, 
while  the  losses  are  too  high  near  the  tip  wall.  The  increase 


of  mesh  density  in  the  radial  direction  from  63  to  96  lines 
improves  the  stagnation  pressure  and  temperature  gradients 
in  the  40-85%  region. 

The  blade-to-blade  plot  of  Mrei  shows  a  passage  shock  which 
is  too  strong. 

The  2D  simulations  of  the  blade-to-blade  section  at  70% 
height  with  the  k-co  and  Baldwin-Lomax  turbulence  models 
do  not  show  significant  differences  in  the  wake  profiles. 
However,  the  pt  contours  obtained  with  Wilcox's  model  are 
smoother  than  with  Baldwin-Lomax's  model;  the  non-zero  pt 
values  are  also  concentrated  ill  the  wakes. 

Chima 


Code 

RVC3CD 

Test  case 

DLR  cascade 

Model  description 

Baldwin-Lonrax  model  with  Ccp  = 

1.216,  Cueb  =  0.646 

Transition  model 

as  for  SWIFT 

Arnone  and  Marconcini 


Code 

TRAF3D 

Test  case 

Rotor  37 

Model  description 

Baldwin-Lonrax's  model,  with 
Cwk=l  -0 

Chima,  Giel,  Boyle's  (1993) 
model, 

Cebeci-Snrith's  model 

Transition  model 

No  transition  model 

The  three  turbulence  models  (Baldwin-Lomax,  Chima-Giel- 
Boyle,  Cebeci-Smith  )  give  the  same  overall  pressure  ratio. 
The  absolute  flow  angle  is  correct  from  10  to  80%  of  the 
blade  span.  In  this  region,  the  stagnation  pressure  and 
temperature  are  too  high  in  the  prediction.  This  suggests  an 
excess  of  the  axial  velocity  component,  induced  by  too  many 
losses  near  the  tip  wall,  as  shown  by  the  efficiency 
distribution. 

The  wake  deficit  is  also  too  high  at  mid-span  in  station  3. 
However,  because  of  the  excess  of  work  input,  the  efficiency 
is  correctly  predicted  in  this  region. 

Dadone  &  De  Palma 


Code 

Test  case 

DLR  cascade 

Model  description 

Baldwin-Lomax's  model  + 
Degani-Schiff  modification 

McNulty 


Code 

ADPAC 

Test  case 

Rotor  37  &  DLR  cascade 

Model  description 

Baldwin-Lomax's  model  +  wall 
functions 

Transition  model 

No  transition  model 

For  Rotor  37,  the  radial  plots  at  98%  of  the  choked  mass  flow 
include  the  relative  velocity  W  and  the  relative  flow  angle  (3. 
There  is  a  low  value  of  W  over  the  whole  span,  and 
simultaneously  a  low  p  value  except  at  mid-span  where  the 
prediction  is  good.  The  blade-to-blade  plots  of  relative  Mach 
number  show  low  predicted  values  downstream  of  the  shock. 
This  induces  a  high  temperature,  particularly  near  the  end 
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walls.  The  axial  velocity  component  must  then  be  low  all 
over  tire  span.  This  reduction  of  the  axial  velocity  is 
estimated  at  -3.8%.  This  is  a  surprisingly  low  value,  that  does 
not  fit  with  tire  mass  flow. 

Tire  radial  plots  at  tire  near  stall  condition  show  a  good 
agreement  with  tire  measured  values,  except  for  P0  under 
30%  of  the  blade  height.  Tire  blade-to-blade  plots  of  Mrei 
show  a  strong  detached  shock,  behind  which  tire  predicted 
level  of  Mrei  is  correct. 

Hall  has  presented  good  results  in  tire  hub  region,  by 
increasing  tire  number  of  nodes  in  tire  radial  direction  from 
49  to  97.  Tire  effect  of  an  upstream  cavity  is  observed  only 
with  this  new  fine  grid,  as  a  consequence  of  a  strong  comer 
stall. 


Calvert  and  Stapleton 


Code 

TRANSCode 

Test  case 

Rotor  37  &  DLR  cascade 

Model  description 

Baldwin-Lomax's  model,  with 
Cwk=l  -0 

Model  Implementation 

Calculation  is  based  on  tire 
nearest  wall 

Boundary  conditions 
and  implementation 

Log-law  is  applied  to  the  viscous 
forces  on  tire  wall  cells  for  y+  >5 

Transition  model 

Tire  flow  is  assumed  fully 
turbulent 

Special  wake  treatment 

The  turbulent  viscosity  is 
assumed  constant  in  the  wake 

A  full  description  of  tire  BL  implementation  is  given  in 
Chapter  1;  in  addition,  the  search  for  Fmax  is  limited  in  order 
to  avoid  spurious  values.  Tire  parameter  Cwk  was  found  to 
have  a  large  effect  on  tire  solutions.  For  example,  for  a 
quasi-3D  solution  of  the  mid-span  section  of  Rotor  37,  tire 
variation  of  CWk  from  0.25  to  1 .0  produces  a  strong  decrease 
of  the  exit  Mach  number  and  hence  air  increase  in  pressure 
ratio  from  1.763  to  2.125,  with  a  corresponding  increase  in 
efficiency  from  82.3%  to  86.3%.  The  value  CWk  =1.0  was 
used  for  the  3D  solutions  and  this  produced  some  over¬ 
prediction  of  temperature  ratio  as  well  as  pressure  ratio. 
Owing  to  the  high  loss  near  tire  tip  wall,  the  axial  velocity  is 
perhaps  overestimated  at  nrid-span. 

The  contour  plots  of  the  turbulent  viscosity  show  a  strong 
influence  of  tire  tip  wall  at  90%  of  the  blade  height,  that  is 
amplified  at  the  near  stall  condition. 

For  both  cases,  calculations  were  also  performed  with  the 
Spalart-Allmaras  (SA)  model  (see  Section  3.6.2).  Looking  at 
the  radial  profile  of  total  pressure  it  seems  that  in  the  tip 
region  the  BL  model  does  a  better  job,  but  the  contours  of 
total  pressure  in  tire  measuring  plane  predicted  by  the  SA 
model  are  closer  to  the  experimental  ones,  with  a  more 
pronounced  distortion  of  tire  wake.  Moreover,  tire  SA  radial 
distribution  of  the  tangential  flow  angle  is  in  fairly  good 
agreement  with  the  experiment,  even  in  the  tip  region.  It  is 
worth  mentioning  that,  with  the  SA  model,  the  code 
converged  better  (even  if  more  slowly)  without  needing  the 
damping  in  the  pt  changes  sometimes  required  by  tire  BL 
model. 


Shabbir 


Code 

VSTAGE 

Test  case 

Rotor  37 

Model  description 

Baldwin-Lomax  algebraic 
turbulence  model 

Transition  model 

transition  at  10%  of  suction 
surface;  fully  turbulent  pressure 
surface 

A  full  set  of  results  is  given  in  Shabbir,  Zhu,  and  Celestina 
(1996).  There  is  a  good  prediction  of  the  radial  variation  of 
tire  absolute  flow  angle  a,  while  the  temperature  is  slightly 
overestimated.  This  is  consistent  with  a  high  axial  velocity  at 
nrid-span,  probably  created  by  a  flow  blockage  near  the  tip 
wall.  See  Section  3.2.3  for  a  comparison  with  two  k-s 
models. 

Couaillicr 


Code 

CANARI 

Test  case 

Rotor  37  &  DLR  cascade 

Model  description 

Michel's  model 

Model  Implementation 

Definition  of  tire  boundary  layer 
thickness  in  term  of  a  specific 
amount  of  tire  vorticity.  Tire 
distance  to  the  wall  is 
determined  from  Butaev's 
formulation. 

Transition  model 

fully  turbulent 

Special  wake  treatment 

no  wake  model 

For  Rotor  37,  the  turbulent  viscosity  disappears  in  tire  wake 
for  the  nresh  points  far  from  tire  casing.  The  losses  generated 
near  the  casing  are  very  large  compared  with  other 
simulations. 

For  comments  on  tire  DLR  cascade  results,  see  Section 
3.6.3. 1. 


3.6.2  One  equation  models 
Calvert  and  Stapleton 


Code 

DRA  TRANSCode 

Test  case 

Rotor  37  &  DLR  cascade 

Model  description 

Spalart-Allmaras'  model 

The  same  boundary  conditions  were  used  as  for  the  Baldwin- 
Lomax  model.  The  kinematic  eddy  viscosity  is  prescribed  at 
the  inlet  according  to  a  typical  turbulent  wall  boundary  layer, 
is  extrapolated  at  the  outlet,  and  tire  profile  boundary  layer  is 
tripped  to  turbulence  near  the  leading  edge.  Tire  Rotor  37 
results  are  changed  significantly  compared  with  the  Baldwin- 
Lomax  solution,  with  tire  pressure  ratio  now  being 
underestimated,  in  a  similar  way  to  the  CANARI  results  with 
this  model.  The  reason  seenrs  to  be  that  the  interaction  of  tire 
passage  shock  with  the  suction  surface  boundary  layer  is 
overestimated  over  tire  outer  half  of  the  span.  The  Spalart- 
Allmaras  model  has  been  used  with  success  on  a  fundamental 
test  case  for  3D  shock-boundary  layer  interaction  (Delery), 
and  so  it  is  possible  that  the  predicted  pre-shock  boundary 
layer  is  incorrect.  Tire  model  does  give  some  indication  of  the 
hole  in  tire  total  pressure  profile  near  tire  hub,  suggesting  a 
better  prediction  of  tire  sensitivity  of  tire  flow  in  this  region. 

For  discussion  of  tire  results  on  the  DLR  cascade,  see  Section 
3.6.1 
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Couaillier 


Code 

CANARI 

Test  case 

Rotor  37 

Model  description 

Spalart-Allmaras'  model 

The  results  have  been  obtained  with  the  same  grid  as  for  the 
mixing  length  and  k-e  models.  Spalart's  model  does  not  show 
any  improvement,  compared  with  the  other  models,  and  does 
not  perform  better  than  the  mixing  length  model. 

Weber 


Code 

OVERFLOW 

Test  case 

Rotor  37 

Model  description 

Spalart-Allmaras'  model  with  a 
modification  for  the  strong  local 
rotation  from  Dacles  and  Mariani 
(1995) 

The  turbulence  production  term  was  calculated  using 
P  =  S-MAX(O.W-S),  where  S  =  strain  rate  tensor,  W= 
vorticity.  Although  S  and  W  are  very  close  in  thin  shear 
layers  (close  to  any  solid  wall,  in  general),  using  S  in  the 
production  term  produced  a  different  solution  from  using  W 
or  this  combination  of  the  two. 

Using  an  overset  O-H  grid,  remarkably  good  predictions  of 
the  spanwise  distributions  of  pressure  and  efficiency  were 
obtained  except  near  the  tip,  without  modelling  the  hub 
leakage  flow. 

3.6.3  Two-equation  models 
3.6.3. 1  Low  Reynolds  number  models 

Couaillier 


Code 

CANARI 

Test  case 

Rotor  37,  DLR  cascade 

Model  description 

k-s  of  Jones-Launder 

Model 

Integration  to  the  wall,  cut-off  for 

Implementation 

k  and  e  (zero  machine)  and 
limitation  of  |it  /  p  >  5000 

Boundary  conditions 

k=0,  e=0  at  the  wall, 
k  and  e  prescribed  upstream, 
k  and  s  extrapolated  downstream 

Initialisation 

From  tire  algebraic  results,  k  and  e 
are  deduced  from  p,  and 
Bradshaw's  relations  with  a  local 
equilibrium  assumption. 

Solution  strategy 

Tire  k-E  equations  are  decoupled 
from  the  mean  flow  equations 

Transition  model 

No 

Special  wake 

Included  in  the  turbulent 

treatment 

equations 

of  the  losses  is  observed  with  the  use  of  the  k-e  model 
compared  with  Michel's  model. 

The  turbulent  viscosity  ratio  reaches  a  high  value  in  the  wake 
at  mid-span  (1200),  and  a  lower  value  near  the  casing  on  the 
pressure  side  (700).  A  small  pocket  of  high  turbulent 
viscosity  occurs  in  the  interaction  of  the  leakage  “vortex”  and 
the  passage  shock. 

Interestingly  enough  Hie  k-e  field  was  decoupled  from  the 
mean  flow  field,  computing  it  with  the  updated  values  of  the 
mean  flow  at  each  time  step.  The  same  strategy  was  applied 
in  the  TRACE-S  code  and  is  probably  done  in  order  to  relieve 
stability  problems. 

For  the  DLR  cascade,  tire  radial  profile  of  circumferential 
flow  angle  seems  better  captured  by  tire  k-e  model;  on  tire 
other  hand  it  predicts  a  less  satisfactory  Pt  profile, 
overestimating  tire  losses.  Peaks  in  tire  isolines  are  better 
defined  with  the  algebraic  model,  particularly  the  loss  core  in 
tire  hub  region;  on  the  contrary,  the  shape  of  the  wake  is  less 
satisfactory.  Also  in  this  case  the  huge  difference  in  tire 
turbulent  viscosity  (using  tire  two-equation  model  it  is  much 
greater)  does  not  produce  a  clear  trend;  for  example  using 
Michel's  model  the  wake  is  a  little  thicker,  but  peaks  in  the 
contour  lines  are  sharper. 

Hah 


Code 

HAH  3D 

Test  case 

Rotor  37  &  DLR  cascade 

Model  description 

k-s  with  Chien's  low  Reynolds 
number  modification  ( 1 982) 

Model  Implementation 

Integration  to  tire  wall  when  y+ 
is  less  than  1 1 .  Wall  functions 
are  used  when  y+  is  greater  than 

1 1 .  No  limitation  or  artificial 
cut-off  for  k  and  s. 

Boundary  conditions 

k=0,  s=0  at  tlie  wall 
k  and  e  are  prescribed  upstream 
using  data  or  local  equilibrium 
conditions 

k  and  e  are  extrapolated 
downstream 

Initialisation 

Uniform  distributions  of  k  and  £ 

Solution  strategy 

Fully  coupled 

Transition  model 

Criteria  based  on  low  Reynolds 
number  modification 

Special  wake  treatment 

None 

For  Rotor  37,  the  turbulence  model  used  with  tire  current 
numerical  procedure  calculates  the  hub-comer  flow 
separation  both  at  98%  and  99%  choked  mass  flow.  The 
numerical  solution  at  99%  mass  flow  agrees  very  well  with 
tlie  measured  data.  The  code  strongly  underpredicts  tire  losses 
at  mid  span,  and  produces  as  a  consequence  too  high  a  level 
of  efficency  and  pressure  ratio. 


For  Rotor  37,  the  k-e  simulation  provides  better  radial 
profiles  of  stagnation  pressure  and  temperature  than  tire 
Michel  model  (see  Section  3.6.1,  Couaillier).  The  level  of 
stagnation  temperature  is  however  slightly  too  high, 
particularly  above  60%  of  tire  blade  span.  This  is  also 
observed  for  tire  stagnation  pressure  which  is  slightly  above 
tire  experimental  values  while  tire  agreement  is  veiy  good  in 
between  30  and  60%  of  tire  blade  height.  Tire  losses  are 
however  too  high  near  tire  tip  wall,  although  a  large  reduction 
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3. 6. 3. 2  Models  with  wall  functions 
Lisiewicz 


Code 

TRACE-S 

Test  case 

DLR  cascade 

Model  description 

Two  versions  of  k-e  model  and 
one  k-co  model 

Model  Implementation 

see  below 

Boundary  conditions 

see  below 

Transition  model 

assumed  fully  turbulent 

Special  wake  treatment 

None 

Three  two-equation  turbulence  models  were  tested:  the 
standard  k-co  (Wilcox)  and  k-e  (Launder  and  Spalding) 
models,  and  the  Kato-Launder  version  of  tire  latter.  Both  the 
k-e  and  k-co  models  use  Spalding's  law  of  tire  wall  in  the  first 
computational  cell;  at  the  inflow  k  is  prescribed  from 
turbulence  intensity  whilst  e  and  co  are  set  specifying  a 
characteristic  turbulence  length  scale;  at  the  outflow 
everything  is  extrapolated.  Hie  extended  k-e  model  and  the 
k-co  model  led  to  practically  equal  radial  profiles. 
Nothwithstanding  this  fact,  the  contours  of  Pt  and  tangential 
angle  (the  most  significant  quantities)  are  similar  but  not 
equal;  the  characteristic  feature  is  that  peaks  are  more 
smeared  using  the  k-co  model,  as  a  consequence  of  the  much 
higher  levels  of  turbulent  viscosity.  The  equal  values  of  the 
pitchwise  integrals  seem  thus  to  indicate  that  there  is  a 
balance  in  the  two  models  between  profile  losses  and  mixing 
losses,  making  the  influence  of  fit  values  somewhat  small. 


Bassi  and  Savini 


Code 

Test  case 

DLR  cascade 

Model  description 

Wilcox  standard  k-co  model 

Martelli 


Code 

FLOS3D 

Test  case 

DLR  cascade 

Model  description 

Wilcox  standard  k-co  model 

Shabbir,  Zhu,  and  Celestina 


Code 

VSTAGE 

Test  case 

Rotor  37 

Model  description 

Two  turbulence  models  are 
used:  SKE  (k-e  of  Launder- 
Spalding,  1974),  CKE  (  k-s 
CMOTT  of  Shih  et  al.,  1995) 
that  avoids  a  local  equilibrium 
assumption 

Boundary  conditions  and 
implementation 

Use  of  wall  functions  of 
Launder-Spalding  (1974)  if 
y+>25  for  flic  first  mesh  point. 

If  y+<ll,  computation  of  tire 
wall  shear  stress  directly  from 
tire  velocity  profile,  and  use  of 
wall  functions  of  Shih  and 

Lumley  ( 1 993). 

Transition  model 

None 

Special  wake  treatment 

None 

The  radial  distributions  of  P0,  T0  and  a  are  greatly  improved 
with  the  CKE  model  compared  with  the  SKE  model.  The 
agreement  with  the  experimental  values  is  good  between  10 
and  70%  of  tire  blade  span  for  tire  CKE  model.  The  k-e  model 
produces  an  improved  description  of  the  wakes  compared 
with  Baldwin-Lomax's  model.  The  two  k-s  models  give 
similar  relative  Mach  number  distributions  in  tire  wake  for 
stations  3  and  4  at  tire  three  radial  positions  examined.  The 
comparisons  with  tire  measurements  are  correct  for  station  4, 
but  show  tire  classical  overprediction  of  tire  wake  deficit  in 
station  3.  As  the  CKE  model  avoids  the  use  of  the  local 
turbulent  equilibrium,  it  could  be  assumed  that  it  has  some 
favourable  effect  on  the  prediction  of  tire  shock-boundary 
layer  interaction.  However,  this  modification  has  almost  no 
effect  on  improving  the  tip  leakage  description.  The  losses 
are  even  higher  in  this  region  for  the  CKE  model. 

The  two  sets  of  boundary  conditions  do  not  show  any 
modification  of  the  results. 

The  SKE  model  is  very  sensitive  to  the  upstream  condition 
for  ,  while  tire  CKE  model  shows  ahnost  no  sensitivity  to 
this  parameter. 

Shabbir  et  al.  (1997)  show  a  strong  influence  of  the  hub 
cavity  on  the  blade  suction  side  comer  stall. 

Ivanovic  and  Hutchinson 


Code 

TASCflow 

Test  case 

Rotor  37 

Model  description 

standard  k-s,  with  wall 
functions. 

Transition  model 

None 

The  stagnation  temperature  shows  a  large  decrease  near  tire 
tip  wall,  with  an  increase  of  the  efficiency.  This  behaviour  is 
not  shown  by  other  contributors* . 

Hildebrandt  and  Vogel 


Code 

TRACE-S,  TASCflow 

Test  case 

Rotor  37 

Model  description 

standard  k-e,  with  wall  functions 

Boundary  conditions 
and  implementation 

20<  y+  <60 

Solution  strategy 

k-s  equations  are  uncoupled 
from  the  mean  flow  equations 

Transition  model 

none 

The  higher  wake  dissipation  observed  in  the  simulation 
between  stations  2  and  3  is  certainly  comiected  with  the 
coarse  mesh  used  in  the  downstream  H  zone. 

From  the  blade-to-blade  turbulent  viscosity  maps,  high  values 
of  |ii/p  are  seen  on  the  suction  side  wake  («1200)  at  mid¬ 
span.  At  90%  span,  lower  values  are  seen  on  the  suction  side 
(«700),  but  they  increase  on  the  pressure  side  near  the 
trailing  edge  (»300).  At  mid-height  near  tire  trailing  edge  in 
station  3,  the  turbulent  viscosity  has  a  significant  value  only 
in  tire  wake  area.  By  contrast,  at  90%  span  tire  flow  feels  tire 
influence  of  tire  tip  wall.  The  increase  of  the  turbulent 
viscosity  on  the  pressure  side  near  the  trailing  edge  may  be 
the  result  of  the  accumulation  of  the  tip  leakage  flow  that 
crosses  tire  blade  passage  from  tire  suction  side. 


1  but  see  footnote  to  Section  3.2.4. 1 
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converged  operating  point  as  an  initial  guess  for  the  next 
operating  point.  However,  the  approach  largely  eliminates 
dependence  on  the  definition  adopted  for  the  residuals:  for 
example,  the  effects  of  features  such  as  Hie  magnitude  of  the 
steps  between  iterations,  residual  smoothing  and  multigrids 
cancel  out,  though  care  is  necessary  if  the  solution  is  started 
on  a  coarse  grid.  Because  of  the  differences  between  the 
definitions,  it  is  not  generally  possible  to  compare  absolute 
levels  of  convergence  from  different  codes. 

Targets  for  Hie  drop  in  residuals  are  generally  between  3  and 
5  decades  relative  to  Hie  initial  values  for  an  arbitrary  first 
guess,  hi  addition,  Hie  iteration  history  of  overall  performance 
parameters  such  as  inlet  and  exit,  mass  flows  and  pressure 
ratio  are  usually  monitored  to  ensure  that  the  solution  has 
reached  a  steady  state,  at  least  in  engineering  terms.  This  can 
be  particularly  important  for  compressor  cases  near  stall, 
where  Hie  solution  may  be  changing  slowly  but  steadily. 

3.7.1 .1  Convergence  Information  For  Contributed  Solutions 
ADPAC  (McNulty) 

(a)  Targets:  Convergence  is  assessed  by  examining 
maximum  and  mis  values  of  an  unspecified  residual. 
No  target  value  has  been  quoted.  In  addition,  Hie 
changes  in  inlet  and  exit  mass  flow  and  in  either 
pressure  ratio  and  efficiency  (Rotor  37)  or  lift  and  loss 
coefficients  and  exit  flow  angle  (DLR  cascade)  are 
monitored. 

(b)  Rotor  37:  The  maximum  and  mis  residuals  dropped  by 
about  27%  and  4  decades  respectively  after  400 
iterations  at  Hie  maximum  efficiency  condition  for  a 
solution  started  on  a  4h  grid  (where  h  is  Hie  cell 
dimension  for  Hie  finest  grid).  Plots  of  Hie  iteration 
histories  indicated  that  inlet  and  exit  mass  flows  were 
constant  and  equal  to  plotting  accuracy  (about  ±0.01%), 
and  that  pressure  ratio  and  efficiency  had  changed  by 
less  than  0.002  and  ±0.02%  respectively  over  the  last 
50  iterations.  No  convergence  data  were  submitted  for 
the  near  stall  point. 

(c)  DLR  cascade:  The  maximum  and  mis  residuals  dropped 
by  about  27%  and  4  decades  respectively  after  350 
iterations  for  a  solution  started  on  a  4h  grid.  Plots  of  the 
iteration  histories  indicated  that  inlet  and  exit  mass 
flows  were  equal  to  within  about  0.1%;  lift  coefficient 
was  constant  to  plotting  accuracy;  and  loss  coefficient 
and  exit  flow  angle  were  dropping  slowly  at  the  rates  of 
about  2%  (of  Hie  current  value)  and  0.1°  for  100 
iterations. 

Bassi/Savini 

(a)  Targets:  A  “stopping  criterion”  of  5  orders  of  magnitude 
decay  for  the  mis  of  Hie  time  derivative  of  density  is 
quoted,  but  Bassi  and  Savini  (1992)  suggest  that  three 
to  four  orders  of  magnitude  is  more  typical. 

(b)  DLR  cascade:  The  convergence  history  was  not 
recorded.  It  is  stated  that  (lie  first  3  orders  of  magnitude 
of  residual  decay  was  obtained  using  multigrid 
technique  on  Hie  mean  flow  variables,  then  Hie  code 
was  run  in  single-grid  mode. 

CANARI  (Couaillier) 

(a)  Targets:  Convergence  is  assessed  by  examining 
maximum  and  mis  residuals  of  each  conservative 
variable.  No  target  value  has  been  quoted,  hi  addition, 
Hie  evolution  of  Hie  mass  flow  rate  at  different  axial 
sections  of  Hie  computational  domain,  including  Hie 
inlet  and  exit  boundaries,  is  monitored. 


(b)  Rotor  37:  For  Hie  solution  at  Hie  maximum  efficiency 
condition  using  the  mixing  length  turbulence  model.  Hie 
mis  residuals  for  all  five  Navier-Stokes  equations 
dropped  by  between  1.5  and  2  decades  after  1500 
iterations  and  then  remained  constant.  Inlet  and  exit 
mass  flow  levels  were  constant  after  about  1000 
iterations,  but  there  was  an  increase  in  mass  flow  of 
about  0.5%  between  inlet  and  exit.  No  convergence 
data  were  submitted  for  Hie  near  stall  point. 

(c)  DLR  cascade:  For  the  solution  using  the  mixing  length 
turbulence  model  the  rms  residuals  for  all  five  Navier- 
Stokes  equations  dropped  by  between  2  and  2.5  decades 
after  1500  iterations  and  then  remained  constant.  Inlet 
and  exit  mass  flow  levels  were  constant  after  about 
1000  iterations,  but  there  was  an  increase  in  mass  flow 
of  about  0.3%  between  inlet  and  exit. 

Dadone/De  Palma 

(a)  Targets:  Convergence  is  assessed  by  examining  the 
residual  of  Hie  continuity  equation.  No  target  value  has 
been  quoted. 

(b)  DLR  cascade:  It  is  stated  that  the  residual  of  the 
continuity  equation  drops  about  2  orders  of  magnitude 
in  12000  multigrid  cycles.  No  plots  of  convergence 
history  have  been  submitted. 

E  URANUS/TURBO  (Kang  and  Hirsch) 

(a)  Targets:  Convergence  state  is  assessed  by  examining 
both  density  residuals  and  mass  flow  convergence 
histories.  No  target  values  are  quoted. 

(b)  Rotor  37:  For  the  solution  at  Hie  maximum  efficiency 
condition,  the  density  residual  dropped  by  about  2 
decades,  and  Hie  inlet  and  exit  mass  flows  were 
constant  to  within  about  0.5%  during  the  last  quarter  of 
the  run.  For  the  final  solution  there  was  a  drop  in  mass 
flow  of  0.39%  between  inlet  and  exit.  No  convergence 
data  were  submitted  for  the  near  stall  point. 

FLOS3D  (Martelli) 

(a)  Targets:  The  residuals  of  all  the  equations,  including 
the  turbulence  model  if  appropriate,  are  considered. 
Convergence  is  assumed  when  a  drop  in  the  residuals  of 
4  orders  of  magnitude  is  obtained  relative  to  an 
arbitrary  first  guess  with  a  linear  pressure  distribution 
between  inlet  and  exit. 

(b)  DLR  Turbine:  The  residuals  for  all  five  flow  equations 
dropped  by  between  27%  and  3  orders  of  magnitude 
within  1500  iterations  and  then  remained  constant. 

HAH3D  (Hah) 

(a)  Targets:  Hie  absolute  values  of  the  residuals  of  each 
finite  difference  equation  are  integrated  over  the  entire 
computational  domain.  Hie  solution  is  considered  to  be 
converged  when  the  total  integrated  residuals  of  all  Hie 
equations  have  been  reduced  by  four  orders  of 
magnitude  from  the  initial  values,  and  the  error  in  mass 
flow  between  inlet  and  exit  is  less  than  0.01%. 

(b)  Rotor  37:  It  is  stated  that  the  momentum  equations  and 
the  turbulence  transport  equations  reached  Hie 
convergence  criterion.  No  plots  of  convergence  history 
have  been  submitted. 

(c)  DLR  cascade:  It  is  stated  that  Hie  momentum  equations 
and  Hie  turbulence  transport  equations  reached  Hie 
convergence  criterion.  No  plots  of  convergence  history 
have  been  submitted. 


70 


OVERFLOW  (Weber) 

(a)  Targets:  Convergence  is  assessed  by  considering  the 
residuals  for  each  grid,  and  also  the  histories  of  changes 
in  mass  flow,  total  pressure,  total  temperature  and 
efficiency.  A  target  of  3  decades  drop  is  quoted  for  the 
residuals. 

(b)  Rotor  37:  The  drop  in  residuals  for  the  background  H- 
grid  did  not  quite  reach  3  decades  for  the  HO  grid.  This 
was  attributed  to  the  high  number  of  interpolations  for 
the  hole  boundary  points  around  the  embedded  rotor 
blade  and  blade  grid.  A  drop  of  about  3  decades  was 
achieved  for  all  grids  in  the  H-O-H  grid  system. 

SWIFT/RVC3D  (Chima) 

(a)  Targets:  Convergence  is  assessed  by  examining 
maximum  and  nns  values  of  tire  change  in  the  variable 
for  the  continuity  equation.  A  target  of  three  decades 
drop  is  set.  hi  addition,  any  changes  in  exit  mass  flow 
and  total  pressure  and  temperature  should  be  only  in 
the  4th  digit,  and  the  exit  flow  profiles  should  be 
converged  to  plotting  accuracy. 

(b)  Rotor  37:  The  maximum  and  rms  residuals  dropped  by 
just  under  2  and  2'A  decades  respectively  at  both  the 
maximum  efficiency  and  stall  conditions.  About  2000 
iterations  were  needed  for  most  operating  points,  but 
tlus  increased  to  over  4000  near  stall.  Plots  of  the 
iteration  histories  of  mean  exit  total  pressure  and 
temperature  indicate  constant  conditions  to  plotting 
accuracy  (about  ±0.0005  of  the  inlet  values). 

(c)  DLR  cascade:  The  maximum  and  rms  residuals  dropped 
by  about  3  and  3'A  decades  respectively  after  2000 
iterations.  Plots  of  the  iteration  histories  indicated  that 
mass  flow  and  total  temperature  were  constant;  mean 
exit  total  pressure  was  dropping  slowly  and  reduced 
from  0.9792  to  0.9780  over  the  last  500  time  steps  (an 
increase  of  5%  in  loss). 

TASCfiow  (Hutchinson  and  Hildebrandt) 

(a)  Targets:  The  time  histories  of  the  nns  residuals  of  the 
U,  V  and  P  equations  are  considered.  No  target  value 
has  been  quoted. 

(b)  Rotor  37:  For  the  solution  with  the  250K  grid  at  95%  of 
the  choke  flow,  the  U  residual  drops  by  2.2  decades,  the 
V  residual  by  3.3  decades  and  the  P  residual  by  3.7 
decades.  No  histories  of  overall  performance 
parameters  have  been  submitted.  No  convergence  data 
were  submitted  for  the  near  stall  point. 

(c)  DLR  cascade:  No  data  currently  submitted. 

TIP3D  (Denton) 

(a)  Targets:  Tire  time  histories  of  the  nns  residuals  for  tire 
meridional  momentum  equation,  the  difference  between 
inlet  and  exit  mass  flow,  and  the  rms  value  of  velocity 
are  considered. 

(b)  Rotor  37:  For  a  solution  at  choked  flow,  the  nns 
residuals  drop  by  about  2.5  decades  after  10000  steps 
and  the  maximum  error  in  continuity  is  0.5%. 

(c)  DLR  cascade:  The  mis  residuals  drop  by  about  2' A 
decades  and  the  error  in  continuity  is  less  than  0.1%. 

TRACE-S  (Lisiewicz,  Vogel  and  Hildebrandt) 

(a)  Targets:  Convergence  is  assessed  by  examining  the 
residual  of  tire  density  (ie  the  change  in  the  variable) 
and,  for  turbulent  calculations,  the  residual  of  the 
square  root  of  both  turbulent  quantities.  The  target  for 
convergence  is  usually  machine  accuracy  (ie  maximum 
residual  in  density  converged  up  to  6  magnitudes, 


averaged  residual  converged  up  to  10  magnitudes,  both 
relative  to  the  maximum  residual  on  time  step  zero). 

(b)  Rotor  37:  For  the  grid  with  500K  points,  the  maximum 
and  nns  density  residuals  dropped  by  about  5  decades 
in  3000  steps  for  the  maximum  flow  point,  when 
started  from  an  initial  guess  with  zero  flow.  The  nns 
residuals  for  the  k  and  e  equations  in  the  turbulence 
model  dropped  by  about  8  decades.  Inlet  and  exit  mass 
flows  were  essentially  constant  after  2000  steps  and 
agreed  to  within  0.01%.  No  convergence  data  were 
submitted  for  the  near  stall  point. 

(c)  DLR  cascade:  The  density  residual  appears  to  reach  a 
limit  of  4  decades  drop  after  about  6500  time  steps. 
Inlet  and  exit  mass  flows  are  constant  and  equal  to 
within  the  plotting  accuracy  (±0.02%)  over  the  last 
2500  time  steps. 

TRAF3D  (Amone  and  Marconcini) 

(a)  Targets:  Convergence  is  assessed  by  examining  a 
residual  defined  as  the  vector  sum  of  the  changes  in  the 
five  conservation  variables,  together  with  parameters 
such  as  inlet  and  exit  absolute  Mach  number.  For 
multigrid  solutions  the  residual  is  based  on  the  changes 
for  the  finest  grid  only.  The  target  for  convergence  is 
stated  as  half  an  order  above  single  precision  machine 
accuracy. 

(b)  Rotor  37:  The  rms  residual  dropped  by  just  over  4 
decades  relative  to  the  level  after  the  first  iteration  in 
300  multigrid  cycles  (including  drops  of  about  0.5 
decades  when  the  grid  was  refined  after  50  and  100 
cycles).  The  maximum  residual  dropped  by  3  decades, 
and  the  inlet  and  exit  Mach  numbers  were  constant  to 
within  about  0.001  over  the  last  100  cycles. 

TRANSCode  (Calvert  and  Stapleton) 

(a)  Targets:  Convergence  is  assessed  by  examining 
maximum  and  nns  values  of  the  imbalance  in  the  axial 
momentum  equation  on  the  fine  grid.  A  drop  of  two 
(and  preferably  three)  decades  is  considered  highly 
desirable.  In  addition,  global  flow  quantities  such  as 
inlet  mass  flow  and  overall  pressure  ratio  or  loss  should 
have  reached  steady  values,  and  mass  flow  should  be 
conserved  to  better  than  0.5%  at  all  quasi-orthogonal 
planes.  Overall  mass  flow  conservation  (inlet  to  exit) 
should  be  better  than  0.05%. 

(b)  Rotor  37:  At  maximum  efficiency,  the  maximum  and 
rms  residuals  dropped  by  3  and  3  'A  decades 
respectively,  including  a  drop  of  about  1  decade  due  to 
applying  heavy  damping  to  tire  changes  in  turbulent 
viscosity  and  time  step  on  the  last  200  steps.  Mass  flow 
was  conserved  to  0.01%  between  inlet  and  exit,  with  a 
maximum  error  at  any  plane  of  0.2%,  and  it  was 
constant  to  within  0.005%  over  the  last  200  steps.  At 
the  near  stall  operating  condition,  the  drops  in  residuals 
were  about  2' A  decades,  with  similar  values  for  the 
other  parameters.  No  damping  of  turbulent  viscosity 
was  needed  with  the  Spalart  and  Allmaras  turbulence 
model. 

(c)  DLR  cascade:  The  maximum  and  rms  residuals  dropped 
by  about  4  decades,  including  a  drop  of  about  1  decade 
due  to  applying  heavy  damping  to  the  changes  in 
turbulent  viscosity  and  time  step  on  the  last  500  steps. 
Mass  flow  was  conserved  to  0.002%  between  inlet  and 
exit,  with  a  maximum  error  at  any  plane  of  0.03%,  and 
it  was  constant  to  within  0.01%  over  the  last  600  steps. 
No  damping  of  turbulent  viscosity  was  needed  with  the 
Spalart  and  Allmaras  turbulence  model. 
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3.6.4  Assessment  of  the  models 

The  following  conclusions  may  be  derived  about  the 

influence  of  the  turbulence  model  on  the  predicted  results  for 

the  NASA  37  rotor: 

a)  Most  of  the  simulations  that  use  Baldwin-Lomax's 
model  produce  very  similar  results  (see  Amone's 
results),  provided  that  enough  nodes  are  used  in  the 
mesh. 

b)  The  mixing-length  models  can  produce  good  results  at 
mid-span.  For  Baldwin-Lomax's  model,  it  is  important 
to  correct  the  original  value  of  tire  constant  CWk  from 
0.25  to  a  value  close  to  1.0.  This  reduces  the  growth  of 
an  important  separation  zone  on  the  suction  side  after 
tire  shock  interaction. 

c)  The  mixing-length  models  give  incorrect  results  for  the 
tip  leakage  flow.  This  class  of  model  has  been  derived 
for  the  description  of  wall  boundary  layers.  The  tip 
leakage  flow  is  a  configuration  with  many  scales  that 
are  not  directly  linked  to  the  local  wall.  The  tip  leakage 
flow  is  dominated  by  an  almost  inviscid  supersonic  jet, 
that  meets  the  supersonic  primary  flow.  At  the 
boundary  of  the  leakage  jet,  there  is  a  shear  layer  that 
comes  from  the  big  difference  of  flow  direction,  starting 
at  the  suction  side  comer  at  tire  blade  extremity.  Tire 
scales  that  are  associated  with  tire  leakage  vortex  ought 
to  be  linked  with  that  shear  layer. 

d)  The  turbulence  models  with  two  transport  equations 
give  improved  results  compared  with  the  mixing  length 
model.  The  improvement  seems,  however,  to  be  rather 
marginal  in  a  two-dimensional  configuration,  as  shown 
by  Chima  with  the  k-co  model  of  Wilcox.  This  suggests 
that  models  with  two  equations  are  better  able  to  deal 
with  strongly  three-dimensional  turbulent  wall  flow. 

e)  All  simulations  are  unable  to  give  a  good  prediction  of 
the  tip  leakage  flow  area,  although  the  k-e  model 
performs  a  little  better.  It  seems  that  the  generation  of 
loss  is  too  high  in  this  region,  where  there  is  a  flow 
separation  from  the  tip  wall  downstream  of  the  passage 
shock.  The  value  of  the  turbulent  viscosity  is  also 
smaller  than  on  tire  suction  side  at  a  lower  radial 
position  along  tire  span.  It  may  be  that  the  interaction 
between  the  leakage  vortex  and  tire  passage  shock 
produces  an  incorrect  reaction  of  tire  turbulence  model. 
This  is  not  surprising  as  most  of  the  models  have  been 
calibrated  for  boundary  layer  flow’s,  without  strong 
three-dimensional  effects  of  tire  leakage  flow  type. 

(f)  The  non-equilibrium  k-s  model  (CMOTT),  with  wall 
law  functions,  and  the  low-Reynolds  number  k-s  model 
give  better  results  than  the  classical  equilibrium  model 
with  w’all  law  functions. 

g)  It  is  difficult  to  deduce  a  firm  conclusion  from  tire 
results  presented  with  the  one-equation  turbulence 
model  (Spalart-Allmaras).  Calvert  and  Couaillier 
obtained  similar  results,  with  the  shock/boundary  layer 
interaction  being  significantly  over-estimated  over  the 
outer  half  of  the  span,  leading  to  under-estimation  of 
tire  pressure  ratio  and  efficiency.  Weber,  with  tire 
model  modified  according  to  Dacles  and  Mariani 
(1995),  achieved  much  better  results  using  two 
different  grid  schemes.  Both  solutions  match  the  test 
results  closely  between  40  and  90%  span,  and  tire 
solution  on  the  overset  O-H  grid  also  matches  the 
measurements  nearer  tire  hub. 


For  the  DLR  cascade,  tire  conclusions  are: 

a)  Moving  from  algebraic  to  one-equation  and  two- 
equation  models,  tire  details  of  tire  flow  field  seem 
generally  better  captured,  as  they  should  be  since  more 
physics  are  included  in  the  turbulence  model,  but 
sometimes  the  integral  values  computed  with  low  level 
approximations  are  closer  to  tire  experimental  data.  As 
normally  differential  turbulence  models  requires 
substantially  more  mass  storage  and  CPU  time  than 
algebraic  ones  and  lead  to  a  more  stiff  system  of 
equations,  tire  question  arises  whether  there  is  a  pay-off 
for  this  increased  complexity.  Tire  results  presented 
here  are  far  from  giving  a  clear  answer,  but  there  is  an 
impression  that  the  good  behaviour  of  a  simple 
turbulence  model  in  complex  flows  relies  on  tire 
somehow  lucky  cancellation  of  opposing  errors  when  an 
integration  is  performed. 

b)  Tire  surprising  differences  showed  between 
computational  results  obtained  with  tire  same 
turbulence  model  can  be  reasonably  ascribed  mainly  to 
their  implementation  and  to  the  artificial  smoothing 
properties  of  each  solver  associated  with  the  grid 
construction.  Tire  conclusions  that  naturally  stem  are 
that  the  best  way  to  analyze  the  relative  performance  of 
each  model  is  to  carry  out  extensive  testing  with  a 
single  type  of  code  and/or  computational  mesh,  and  that 
detailed  comparisons  between  the  results  from  similar 
codes  are  an  important  aid  to  code  developers. 


3.7  CONVERGENCE  CRITERIA 

3.7.1  Convergence  levels  achieved 

Convergence  of  a  numerical  scheme  can  mean  either 

(i)  that  the  solution  to  the  finite  difference  equations 
approaches  the  true  solution  to  tire  partial  differential 
equations  as  the  mesh  is  refined;  or 

(ii)  that  the  iterative  process  has  been  repeated  until  the 
magnitude  of  the  difference  between  tire  function  at  tire 
n+1  and  n  iteration  levels  is  as  small  as  required  at 
every  grid  point. 

Definition  (i)  requires  that  the  scheme  is  stable  and 
consistent.  The  second  definition  is  generally  adopted  when 
considering  engineering  applications  of  CFD  and  it  is  this 
“iteration  convergence”  which  is  considered  here. 

Ideally  the  criterion  for  convergence  would  be  that  tire 
differences  between  successive  iteration  levels  (usually 
termed  residuals)  should  be  reduced  to  the  level  implied  by 
double  precision  machine  accuracy,  and  it  is  highly  desirable 
that  codes  are  demonstrated  to  be  capable  of  this  for  simple 
flows.  However,  this  level  of  convergence  is  often  not 
possible  for  more  complex  flows  because  of  some  local 
instability,  which  may  be  linked  to  a  real  physical 
phenomenon  such  as  vortex  shedding  or  due  to  a  numerical 
feature.  Tire  “correct”  response  to  such  instabilities  is  to 
investigate  them  using  a  time-accurate  method,  but  this  is  not 
practical  with  present  computing  resources. 

Most  practitioners  therefore  adopt  more  pragmatic  measures 
to  assess  whether  a  solution  is  sufficiently  converged  to  be 
useful.  Typically,  convergence  is  assumed  to  have  occurred 
when  the  residuals  have  dropped  to  a  given  proportion  of  tire 
levels  at  tire  start  of  tire  solution.  This  indicator  has  the 
disadvantage  that  it  is  dependent  on  tire  initial  conditions 
assumed:  an  improvement  to  tire  accuracy  of  tire  first  guess 
gives  an  apparent  drop  in  tire  convergence  level  attained  and 
this  must  be  taken  into  account  when  using  air  already 
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VSTAGE  (Shabbir/Zhu/Celestina) 

(a)  Targets:  Convergence  is  determined  by  monitoring  the 
maximum  and  mis  residuals  of  density.  Typically  a 
drop  of  two  decades  is  needed.  The  global  parameters, 
such  as  the  mass  flow  at  the  inlet  and  exit  planes, 
pressure  ratio  and  efficiency,  are  also  monitored  to 
ensure  that  these  have  reached  steady  state  values.  The 
mass  flow  at  the  inlet  and  exit  should  match  within 
0.02%,  and  they  should  both  remain  constant  to  within 
0.005%  for  the  last  few  hundred  iterations.  Mass 
conservation  at  all  planes  is  typically  better  than  0.5%. 

(b)  Rotor  37:  It  is  stated  that  the  mass  flow  at  the  inlet  and 
exit  matched  to  within  0.017%  and  that  both  were 
constant  to  within  0.004%  for  the  last  200  iterations  at 
the  high  flow  point  (m/mchoke  =  0.98).  No  plots  of 
convergence  history  have  been  submitted. 

3.7.2  Assessment  of  convergence  parameters 
The  convergence  levels  actually  achieved  for  the  WG  26  test 
cases  did  not  reach  the  nominal  targets  for  many  of  the 
solutions  submitted  (see  Table  3.1).  The  levels  achieved 
mainly  ranged  from  2  to  4  decades  for  both  Rotor  37  and  the 
DLR  cascade.  Little  information  was  submitted  for  Rotor  37 
at  the  near  stall  condition,  but  this  tended  to  indicate  slightly 
poorer  levels  of  convergence  than  at  points  nearer  choke.  Tire 
iteration  histories  indicated  that  the  solutions  had,  in  most 
cases,  reached  sensibly  constant  conditions  for  the  main 
overall  performance  parameters.  A  few  of  the  contributors 
stated  that  they  carried  out  additional  investigations  into 
whether  more  detailed  aspects  of  the  solution  (such  as  the 
exit  radial  pressure  and  temperature  profiles)  were  also 
constant,  but  no  examples  were  submitted. 

The  best  level  of  convergence  submitted  to  the  Working 
Group  was  for  the  TRACE-S  solution  of  Rotor  37  at 
maximum  flow  (see  Fig  3.41).  This  achieved  drops  of  about  5 
decades  in  the  density  residuals  and  of  about  8  decades  in  the 
turbulence  model  equations  for  k  and  e.  A  more  typical  result 
is  the  TRANSCode  solution  for  the  DLR  turbine  cascade  (Fig 
3.42)  where  there  is  a  drop  of  only  about  2.5  decades,  with  a 
further  improvement  of  one  decade  if  heavy  damping  is 
applied  to  the  values  of  turbulent  viscosity  calculated  by  the 
Baldwin-Lomax  turbulence  model.  However,  the  overall 
performance  parameters  are  completely  steady  and  the 
maximum  error  in  mass  flow  conservation  at  any  plane  is  less 
than  0.03%. 

To  summarise,  the  convergence  levels  achieved  for  the  WG 
26  test  cases  generally  fell  short  of  the  ideal.  No  tune- 
accurate  solutions  were  presented  and  it  is  not  clear  whether 
the  limits  reached  were  due  to  physical  phenomena  such  as 
vortex  shedding  or  to  numerical  features.  However,  most  of 
the  solutions  were  sensibly  converged  in  terms  of  overall 
performance  parameters,  and  they  are  probably  typical  of 
those  generally  produced  for  high-speed  turbomachinery 
blade  rows.  Given  this  situation  it  is  suggested  that 
experienced  operators  are  necessary  for  applying  current  CFD 
codes  to  turbomachinery  and  for  assessing  tire  adequacy  of 
the  results. 
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Fig  3.41  Convergence  history  for  TRACE-S  solution  for  Rotor  37  (Hildebrandt) 
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Fig  3.42  Convergence  history  forTRANSCode  solution  for  DLR  cascade  (Stapleton) 
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4.1  CONCLUSIONS 

Turbomachinery  flows  embody  many  complex  physical 
phenomena,  which  Reynolds-averaged  Navier-Stokes  codes 
are  able  to  predict  to  some  extent.  The  key  to  improving  the 
predictions  is  to  understand  the  physical  phenomena  and  the 
features  of  the  codes  which  model  them.  Different  aspects  of 
the  code  are  important,  depending  on  the  physical  nature  of 
the  flows  being  predicted;  it  follows  that  what  is  best  for 
predicting  external  flows  around  aircraft  or  for 
meteorological  predictions  may  not  be  best  for 
turbomachines.  The  two  test  cases  studied  in  detail  enabled 
the  WG  to  assess  the  ability  of  CFD  to  predict  several  of  the 
aerodynamic  phenomena  which  dominate  the  performance  of 
turbomachines.  The  Rotor  37  case  included  shock  waves, 
comer  stall,  and  tip  clearance  effects  in  a  compressor  rotor 
with  supersonic  inlet  relative  Mach  number  across  the  whole 
span.  The  DLR  annular  cascade  showed  the  effect  of 
spanwise  pressure  gradients  on  the  complex  secondary  flow 
field.  The  test  cases  did  not,  on  the  other  hand,  include 
important  transition  effects,  heat  transfer,  or  high  incidence 
conditions,  nor  of  course  the  three-dimensional  separations 
found  within  a  centrifugal  impeller.  More  importantly,  the 
WG  did  not  attempt  to  cover  the  difficult  area  of  row-to-row 
interference  due  to  unsteady  flow. 

In  both  cases,  it  proved  difficult  to  clarify  all  the  details  of 
the  flow  from  the  measurements  alone,  but  with  the  help  of 
the  CFD  results  a  good  understanding  of  all  the  flow  fields 
was  reached.  The  use  of  computer  graphics  to  display 
streamlines  in  complex  three  dimensional  flows  was 
particularly  helpful.  The  secondary  flow  in  a  turbine  cascade 
has  been  well  known  for  many  years,  but  the  complexities  of 
the  Rotor  37  flow  have  only  now  been  illuminated  by 
experiments  and  CFD  studies  conducted  by  NASA  Lewis. 
This  research  was  not  specifically  done  for  the  WG,  but  the 
widespread  international  interest  created  by  the  ASME 
Turbomachinery  Committee’s  exercise  and  now  by  the  WG 
have  spurred  on  NASA’s  research. 

The  overall  performance  of  both  test  cases  predicted  by  a 
wide  range  of  CFD  codes  fell  short  of  the  accuracy  engine 
designers  need  as  a  basis  for  decision-making.  It  is  well 
known  that  designers  “calibrate”  the  predictions  of  their 
current  codes  against  the  measured  performance  of  their 
company’s  similar  compressors  and  turbines;  but  CFD 
specialists  obviously  need  to  identify  the  reasons  for  the 
inaccuracies  and  so  eliminate  them. 

4.1.1  Rotor  37 

The  NASA  Rotor  37  case  was  chosen  by  ASME  in  1993  as  a 
“blind”  test  case  offered  to  the  CFD  community.  Most  codes 
predicted  correctly  the  shape  of  the  performance  curves,  but 
too  high  a  pressure  ratio,  and  too  low  an  efficiency  (Figs  3.1 
and  3.2).  The  measurements  suggest  there  is  a  comer  stall 
which  prevents  the  rotor  achieving  the  predicted  pressure 
ratio.  The  error  in  efficiency  prediction,  on  the  other  hand, 
arises  largely  because  the  codes  overestimate  the  tip 
clearance  loss. 

There  is  a  shock  wave  at  the  rotor  leading  edge,  and  the 
predicted  Mach  number  in  the  blade  passage  after  it  is 


always  too  low.  The  ability  to  predict  shock/boundary  layer 
interaction  correctly  is  the  issue  here.  It  has  been  shown  that 
this  result  is  sensitive  to  the  local  mass  flow.  Possibly,  the 
local  mass  flow  predicted  by  the  codes  and  deduced  from  the 
experiments  do  not  agree,  either  because  the  overall 
measured  mass  flow  is  wrong,  or  because  inaccuracies  in 
predicting  the  flow  in  the  endwall  regions  has  led  to  errors  in 
the  spanwise  distribution  of  the  flow,  and  hence  an 
underestimate  of  the  axial  velocity  around  mid-span.  A  more 
detailed  study  of  this  question  is  suggested. 

All  the  predictions  supplied  to  ASME,  although  generally 
encouraging,  failed  to  reproduce  the  hub  comer  stall 
suggested  by  the  measurements,  and  nor  did  most  of  the 
more  recent  predictions  supplied  to  the  WG.  This  led  NASA 
to  make  an  important  discovery:  the  flow  in  the  hub  region 
was  seriously  affected  by  the  presence  of  a  small  axial  gap 
(0.75  mm)  in  the  hub  annulus  line  just  upstream  of  the  rotor 
leading  edge  (because  the  hub  wall  ahead  of  that  point  did 
not  rotate).  Shabbir  et  al  (1997)  demonstrated  by  both 
measurements  and  computations  that,  although  the  gap  led 
only  to  a  blind  cavity,  air  was  pumped  in  and  out  of  it  as  the 
rotor  rotated.  The  presence  of  a  shock  wave  even  at  the  hub 
results  in  sufficiently  large  pitchwise  static  pressure  changes 
to  cause  the  pumping.  Shabbir  et  al  (1997)  showed  that  CFD 
predictions  in  which  the  inflow-outflow  was  simulated 
predicted  the  comer  stall,  while  the  same  code  without  the 
inflow-outflow  did  not.  The  important  effect  of  the  gap  on 
the  comer  stall  is  a  major  discovery,  but  it  came  at  a  late 
stage  of  the  WG’s  activities,  and  only  one  other  member  had 
time  to  try  to  model  the  gap-flow  in  his  code.  Only  a  few 
codes  predicted  a  comer  stall  without  modelling  the  flow  in 
and  out  of  the  gap. 

Fig  3.6  shows  the  comer  stall,  as  predicted  by  one  of  the 
codes.  It  is  evident  that  air  is  convected  away  from  the  hub 
in  the  separated  region  on  the  suction  surface.  The  effect  of 
the  comer  stall  on  the  pitchwise-mean  pressure  distribution 
is  a  local  depression  in  total  pressure  around  20%  span, 
visible  in  Fig  3.3,  a  depression  which  most  of  the  codes  do 
not  predict.  Comer  stall  has  been  observed  in  many 
compressors,  and  there  is  no  suggestion  that  it  only  arises 
when  there  is  a  gap  or  when  the  flow  is  transonic,  although 
in  this  particular  case  the  comer  stall  disappears  at  speeds 
below  design  speed  (when  the  inlet  relative  velocity  at  the 
hub  is  subsonic).  It  has  been  postulated  (Povinelli,  1997)  that 
the  interaction  of  the  glancing  shock  wave  from  the  rotor 
pressure  surface  with  the  hub  boundary  layer  has  some 
influence  on  the  generation  of  secondary  flow  and  the  total 
pressure  loss  at  the  hub  surface.  This  influence  disappears  at 
subsonic  operating  conditions. 

One  or  two  of  the  solutions  supplied  to  the  WG  predicted  the 
comer  stall  without  simulating  the  gap.  It  is  probably  very 
sensitive  to  the  mesh  topology  and/or  the  grid  density  as  well 
as  the  turbulence  model. 

The  flow  pattern  in  the  tip  clearance  region  has  been  studied 
by  Chima  (1996b)  and  by  Suder  and  Celestina  (1996).  They 
showed  that  the  interaction  between  the  overtip  flow  and  the 
mainstream  generated  very  high  local  shear,  a  region  in 
which  the  axial  velocity  is  reversed,  and  a  region  in  which 
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the  casing  wall  boundary  layer  is  separated.  This  complex 
flow  pattern  proved  too  difficult  for  most  turbulence  models, 
which  tended  to  predict  too  low  a  turbulent  viscosity  very 
near  the  wall,  giving  rise  to  the  separation,  and  hence  too 
high  pressure  losses  in  the  outermost  5%  of  the  span  outside 
the  casing  boundary  layer. 

The  blade  wakes  even  in  the  quasi-two-dimensional  region 
of  the  flow  near  mid-span  are  inaccurately  predicted  (Fig 
3.15).  Near  the  trailing  edge  they  are  too  deep  and  narrow, 
but  as  they  move  downstream  they  dissipate  too  fast  in  some 
solutions,  so  that  by  the  downstream  traverse  plane  they  are 
too  shallow  and  too  wide.  This  dissipation  results  from  using 
too  coarse  a  grid  in  that  region,  or  a  grid  which  is  not  aligned 
with  the  stream. 

4.1.2  DLR  cascade 

The  overall  pressure  loss  predicted  by  some  of  the  codes  was 
as  much  as  40%  in  error.  The  measurements  of  total  pressure 
show  that  the  secondary  loss  region  near  the  hub  is  confined 
close  to  the  wall  by  the  spanwise  static  pressure  gradient, 
and  that  the  location  of  the  peak  loss  is  generally  well 
predicted,  though  its  magnitude  varies  very  widely  between 
the  codes.  The  secondary  loss  region  at  the  casing,  on  the 
other  hand,  is  convected  well  away  from  the  wall,  the 
measured  loss  peak  being  at  55%  span  at  the  measurement 
plane  (40%  chord  downstream).  Most  codes  were  unable  to 
predict  the  location  and  magnitude  of  this  peak.  So  it  is  clear 
that  the  codes  are  unable  to  predict  correctly  the  highly 
three-dimensional  secondary  flow  and  the  resulting 
migration  of  low  energy  fluid.  Nevertheless,  the  solutions 
implementing  turbulent  transport  models  were  generally 
closer  to  the  experimental  results  than  solutions 
implementing  algebraic  turbulence  models,  provided  fine 
grids  were  used  in  the  radial  direction. 

One  code  (TRACE-S)  was  run  with  the  same  turbulence 
model  (k-e)  but  two  different  grid  topologies.  Both  the 
pressure  losses  and  the  flow  angles  changed  appreciably,  and 
exposed  the  difficulties  introduced  at  the  interface  between 
an  O-grid  and  an  H-grid  just  behind  the  trailing  edge. 

4.1.3  Choice  of  algorithms 

The  algorithms  used  for  the  WG  were  mostly  of  the  time¬ 
marching  type,  but  pressure  correction  methods  were  also 
represented.  While  the  details  of  the  algorithm  must  surely 
control  the  stability  and  convergence  of  the  code  as  well  as 
its  running  time,  the  WG  has  no  evidence  to  suggest  that  it 
has  any  effect  on  the  accuracy  of  the  converged  result. 

4.1.4  Grid  construction 

The  desirable  features  of  a  computational  grid  are  well 
known:  it  should  be  fine,  have  approximately  square  cells, 
and  be  aligned  with  the  stream.  Since  these  properties  are 
impossible  to  achieve  simultaneously  in  a  turbomachine 
context,  the  choice  of  grid  always  represents  a  compromise 
between  the  various  desirable  properties  of  the  grid  and  the 
complications  introduced  into  the  algorithms  when  complex 
mixed  grid  schemes  are  chosen.  It  is  also  known  that  the 
choice  cannot  be  dissociated  from  the  algorithm  and  the 
turbulence  model.  Nevertheless,  some  general  conclusions 
are  possible  in  the  case  of  a  structured  grid.  No  unstructured 
grid  solutions  were  offered  to  the  WG. 

In  the  present  test  cases,  no  single  grid  type  stood  out  as 
being  superior  to  the  others.  In  general,  the  O-  and  C-grids 
proved  better  in  the  leading  edge  region,  although  the 


present  test  cases,  being  at  nearly  zero  incidence,  were 
relatively  insensitive  to  the  leading  edge  region.  In  neither 
test  case  was  it  possible  to  measure  the  surface  boundary 
layers,  so  the  influence  of  the  grid  construction  on  profile 
loss  could  not  be  quantified. 

At  the  trailing  edge,  the  C-,  I-,  and  H-grids  that  aligned  with 
the  wake  provided  the  best  wake  definition.  The  O-grid 
solutions  provided  high  resolution  of  the  flow  near  the 
trailing  edge  but  diffused  the  wake  too  rapidly  as  the  grid 
opened  up  further  downstream. 

The  most  successful  grids  used  by  contributors  were  an  H-I- 
H  grid  (because  it  aligned  with  the  shock  waves)  (Fig  3.36) 
and  an  overset  O-H  grid  (Fig  3.39),  which  avoids  some  of 
the  problems  which  have  been  identified  in  other 
configurations.  Both  were  applied  to  Rotor  37. 

The  grid  lines  must  be  clustered  progressively  near  solid 
surfaces,  and  ought  to  be  clustered  in  any  regions  of  strong 
shear.  All  the  solutions  submitted  used  grids  well  clustered 
near  surfaces,  'flic  total  number  of  cells  varied  widely,  with 
several  contributors  conducting  grid  refinement  studies  to 
establish  how  many  cells  were  needed  to  make  the  solution 
grid-independent  to  engineering  accuracy.  The  minimum 
number  of  cells  must  depend  on  the  flow  being  computed, 
and  the  algorithm  and  turbulence  model,  so  it  is  difficult  to 
generalise.  For  examining  localised  flow  features,  such  as 
leading  edge  separation  bubbles,  the  grid  must  be  locally 
refined.  However,  it  became  clear  that  in  the  cases  used  by 
the  WG  at  least  50  grid  lines  hub-to-tip,  around  50  blade-to- 
blade  (if  wall  functions  are  used),  and  around  300,000  cells 
in  all  are  needed  if  the  pitchwise-mean  performance  is  to  be 
resolved.  If  wall  functions  are  not  used,  a  finer  grid  is 
needed  near  the  walls.  To  capture  the  three-dimensional 
detail  of  the  secondary  flow  vortices  in  the  DLR  cascade, 
more  than  100  grid  lines  hub-to-tip  may  be  needed.  For 
overall  performance  and  for  blade  surface  pressure 
distributions,  on  the  other  hand,  a  grid  of  around  200,000 
cells  may  be  adequate. 

However,  the  WG  believes  that  having  a  large  number  of 
grid  points  (perhaps  one  million  per  row)  is  not  necessarily 
sufficient  to  obtain  an  accurate  solution.  The  turbulence 
model  must  also  be  adequate. 

The  “pinched  tip”  model  of  the  tip  clearance  region  chosen 
by  some  contributors  is  unsatisfactory,  in  that  the  actual 
clearance  is  not  used;  some  empirical  “effective”  clearance 
is  chosen  instead,  which  may  be  dependent  on  the  turbulence 
model.  From  the  fully-gridded  solutions  submitted,  which 
were  in  this  Rotor  37  case  just  as  inaccurate  as  the  pinched 
tip  solutions,  the  WG  was  unable  to  recommend  a  minimum 
number  of  cells  within  the  clearance  region.  Some  WG 
members  suggested  about  ten,  on  the  basis  of  other 
experience. 

4.1.5  Turbulence  model 

Turbulence  models  can  broadly  be  divided  into  mixing 
length  types  (including  the  popular  Baldwin-Lomax  model) 
and  turbulent  convection  types  (including  k-e  models  and 
one-equation  models).  Most  of  the  Rotor  37  solutions 
contributed  to  the  WG  using  the  Baldwin-Lomax  model 
produced  very  similar  results,  and  they  were  generally 
inferior  to  the  predictions  obtained  using  turbulent 
convection  models  in  those  regions  where  the  flow  is 
separated  or  highly  three-dimensional.  It  is  logical  that  any 
turbulence  model  requiring  a  “distance  to  the  nearest  wall” 
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as  a  way  of  defining  the  shape  of  the  whole  boundary  layer, 
and  tuned  to  predict  two-dimensional  turbulent  boundary 
layers,  must  encounter  serious  difficulties  in  such  a  region. 
Some  turbulent  transport  models  also  use  the  “distance  from 
the  wall”,  but  only  to  adjust  the  calculation  within  the 
laminar  sub-layer,  so  that  objection  is  no  longer  valid.  The 
WG  confirmed  that  the  turbulent  convection  models  used  in 
the  present  study  tended  to  give  better  solutions  in  separated 
flow  regions  than  the  mixing  length  models. 

For  flows  which  are  subsonic  and  nearly  two-dimensional, 
and  where  the  viscous  phenomena  are  primarily  of  the  nature 
of  a  boundary  layer,  any  of  the  well-known  turbulence 
models  are  adequate,  since  they  were  set  up  for  boundary 
layers.  But  the  grid  needs  to  be  adequately  fine  near  the 
walls.  It  is  generally  considered  that  for  algebraic  models  a 
y+  value  less  than  5  is  advisable  unless  wall  functions  are 
used;  for  codes  in  which  the  calculation  extends  fully  to  the 
wall  using  a  low  Reynolds  number  turbulence  model  a  y+ 
value  less  than  1  is  needed  (and  as  low  as  0.1  for  heat 
transfer  calculations),  though  no  contributor  used  as  low  a 
value  as  that.  In  codes  using  wall  functions,  y+  =  50  seems 
adequate. 

Massive  differences  were  noted  between  the  values  of  eddy 
viscosity  predicted  by  different  turbulence  models.  These 
differences,  together  with  locally  coarse  grids  in  some 
solutions,  led  to  differences  in  loss  prediction  which  took 
some  of  the  losses  outside  an  acceptable  range  of  accuracy. 
Some  of  the  difference  may  result  not  from  the  modelling 
concept  but  from  the  way  it  is  implemented  within  the 
particular  code.  Although  most  of  the  Baldwin-Lomax 
solutions  for  Rotor  37  were  similar,  this  was  not  so  for  the 
DLR  cascade  solutions;  and  more  generally  solutions  using 
nominally  the  same  turbulence  model  did  not  always  agree 
with  each  other. 

The  WG  was  unable  to  identify  any  one  turbulence  model 
which  always  gave  good  loss  predictions.  It  is  well  known 
that  this  is  an  area  of  continuing  vigorous  research,  and  it 
needs  to  be. 

Transition  predictions  were  not  thought  to  be  important  for 
the  test  cases  chosen  by  the  WG,  but  it  is  well  known  that 
many  current  CFD  codes  cannot  predict  transition  or  re¬ 
attachment  satisfactorily.  Predictions  of  the  flow 
downstream  of  a  leading  edge  separation  bubble  are  in  some 
applications  critical  to  aerodynamic  loss  prediction. 
Transition  prediction  is  also  the  key  to  good  heat  transfer 
prediction. 

4.1.6  Convergence 

Code  developers  generally  aim  to  reduce  chosen  residuals  by 
around  five  decades,  but  in  general  the  solutions  contributed 
to  the  WG  only  achieved  between  two  and  three  decades. 
Users  also  check  that  overall  performance  parameters  have 
stabilised.  It  is  not  clear  whether  the  inability  of  the  codes  to 
converge  better  is  due  to  numerical  problems  (induced  by  the 
mesh  topology  near  to  the  wall  and  in  specific  regions  such 
as  the  leading  and  trailing  edges)  or  whether  it  is  the  result 
of  inadequate  flow  modelling  (such  as  a  turbulence  model  or 
the  forcing  of  a  steady  solution  to  a  flow  field  which  is 
known  to  contain  unsteady  shed  vortices).  So  experienced 
code  users  are  needed  to  assess  the  adequacy  of  results. 

4.1.7  Working  Group  procedures 

The  key  activity  of  the  WG  was  to  analyse  and  interpret  the 
results  of  the  computations.  The  full  result  of  each 


computation  was  several  million  numbers,  far  too  large  a 
data  file  to  be  circulated  as  a  floppy  disc  or  an  e-mail.  The 
WG  therefore  chose  to  start  by  examining  the  pitchwise- 
mean  solutions  at  the  traverse  planes  at  which  the 
measurements  had  been  made.  Those  solutions  were  sent  to 
one  member,  who  plotted  them  together.  However,  the  WG 
members  who  undertook  the  detailed  analyses  found  that 
they  needed  selective  plots  of  parameters,  at  other  planes 
too.  In  retrospect,  the  establishment  of  a  general  data  base  of 
the  selected  results,  which  all  members  could  access  (using 
e-mail  or  perhaps  by  Internet)  would  have  helped  the  WG. 

4.1.8  Closing  remarks 

The  aim  of  the  WG  was  to  understand  why  current  CFD 
codes  are  sometimes  unable  to  predict  the  measurements 
made  even  on  isolated  turbomachine  stages,  and  to  clarify 
the  role  of  the  grids  and  turbulence  models  used  in  achieving 
good  predictions.  These  are  not  simple  questions  to  which  a 
final  answer  could  be  expected,  but  the  Group  did  throw 
light  on  several  aspects  of  turbomachinery  CFD,  which  may 
point  research  workers  in  the  right  direction  for  the  future. 
The  WG  members  most  deeply  involved  in  the  analysis 
agreed  how  valuable  the  study  had  been  to  them.  The 
improvement  in  the  quality  of  the  Rotor  37  predictions 
resulting  from  the  ASME  exercise,  the  subsequent  NASA 
research  and  the  WG  activities  is  obvious.  The  experience  of 
the  WG  has  also  provided  a  timely  reminder  that  good 
quality  detailed  experimental  measurements  are  essential  to 
the  continuing  development  of  CFD. 

4.2  RECOMMENDATIONS  FOR  RESEARCH 

These  recommendations  apply  both  to  code  developers  and 
code  users. 

1.  Turbomachinery  CFD  should  turn  away  from 
mixing  length  turbulence  models  and  develop  the 
use  of  turbulent  transport  models. 

2.  If  losses  are  to  be  predicted,  around  300,000  well- 
chosen  grid  points  are  needed  per  blade  row, 
refined  in  regions  of  high  aerodynamic  shear. 

3.  If  no  wall  function  scheme  is  used,  a  finer  grid  is 
essential  near  the  walls,  and  hence  many  more  grid 
points. 

4.  If  fuller  details  of  the  three-dimensional  flow 
pattern  are  needed,  a  finer  radial  grid  and  hence 
500,000  or  more  points  in  all  are  needed, 
depending  on  the  turbulence  model. 

5.  Some  intensive  research  should  be  focused  on  tip 
clearance  effects. 

6.  Computations  allowing  for  full  geometric  details, 
leakage  flows  and  annulus  wall  gaps  should  be 
more  widely  developed. 

7.  More  use  should  be  made  of  computer  graphics  to 
visualise  complex  three  dimensional  flows. 

8.  There  is  still  a  need  for  detailed  experimental 
measurements  as  a  basis  for  future  developments  in 
turbomachinery  CFD,  especially  on  multiple  blade 
rows. 

4.3  RECOMMENDATIONS  TO  PEP 

1 .  Detailed  experiments  to  provide  test  cases  are  very 
expensive.  PEP  should  encourage  international 
collaboration  in  devising  and  undertaking  suitable 
experiments. 

2.  Similar  studies  should  be  undertaken,  specifically 
on  heat  transfer  predictions  in  turbomachines  and 
later  on  combustor  flows.  Full  use  should  be  made 


of  electronic  data  transfer  facilities  for  conducting 
those  studies. 

In  about  five  years’  time,  a  Working  Group  should 
be  set  up  to  return  to  the  turbomachinery  CFD  area, 
this  time  choosing  test  cases  with  more  than  one 
row,  and  focusing  attention  on  the  way  steady 
codes  can  represent  unsteady  interference  effects. 
By  then,  sufficient  unsteady  computations  will  be 
available. 
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